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FLEXIBLE ROBOT CONTROL: MODELING AND EXPERIMENTS 

By 

Irving J. Oppenheim 
Camegle-Mellon University 
Pittsburgh. Pennsylvania 

and 

Isao Shimoyama 
University of Tokyo 
Tokyo. JAPAN 


ABSTRACT 


A dynamic model fills several roles In the development of flexible manipulators and their 
control structures. A proper dynamic model permits Identification of the proper state 
variables for control, completes the mathematical model used in design studies and In 
simulation, and provides the forward transform needed In model-based control. While 
there exist many proven analytical approaches, and although numerous models have been 
constructed and tested, there remains a need for simple models which capture all the 
important behavior while otherwise suppressing modeling complexities and 
computational demands. Such simple models are necessary for online applications 
because of their computational compactness, and are advantageous for design and 
simulation studies because of their accessibility by users. For manipulator control 
applications, an Ideal (simple) model might contain Independent variables no greater in 
number than the state variables required for acceptable control. This paper describes such 
a model and Its use In experimental studies of flexible manipulators. 

The analytical model developed In this research uses the equivalent of Rayleigh’s method 
lo approximate (he displaced shape of a flexible link as the static elastic displacement 
which would occur under end rotations as applied at the Joints. The generalized 
coordinates are thereby expressly compatible with joint motions and rotations In serial 
link manipulators, because the amplitude variables are simply the end rotations between 
the flexible link and the chord connecting the end points. The equations for the system 
dynamics are quite simple and can readily be formulated for the multi-link, three- 
dimensional case. When the flexible links possess mass and (polar moment of) Inertia 
which are small compared to the concentrated mass and Inertia at the Joints, the 
analytical model is exact and displays the additional advantage of reduction In system 
dimension for the governing equations. 

Four series of pilot tests have been completed. Studies on a planar single-link system were 
conducted at Camegle-Mellon University, and tests conducted at Toshiba Corporation on 
a planar two-link system were then Incorporated into the study. A single link system 
under three-dimensional motion, displaying biaxial flexure, was then tested at Camegle- 
Mellon. The most recent tests, also conducted at Camegle-Mellon, studied a three- 
dimensional system In which coupled (biaxial) flexural-torsional vibrations were present. 
In every test series effective control of the flexible system was accomplished; performance 
of the proposed model was studied and confirmed. 
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FLEXIBLE ROBOT CONTROL: 
MODELLING AND EXPERIMENTS 


Irving J. Oppenheim, Carnegie-Mellon University 
Isao Shimoyama, University of Tokyo 


Describing a simple dynamic model: 

• Useful for rapid prototyping and control 
system development 

• Useful during manipulator design 

• Applicable for real-time computation 


Describing experimental results: 

• Single-link 2-D 

• Single-link 3-D 

• Two-link 2-D 

• Two-link 3-D 
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Modelling Link Flexibility Effects 


Problems: 

• Manipulators are non-linear by their 
configuration 

• All models for flexible dynamics must 
approximate the solutions to PDE’s 

• Generalized co-ordinates (mode shapes) are 
often utilized 

• Truncated mode shape models: OK, but not 
fully consistent with manipulator control 


Demands: 

• Generalizable to M-DOF manipulators 

• Simple to formulate and use in simulation 

• Computable in real-time 
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Intended Users 


Laboratory research in flexible manipulator control 

These restrictions are common: 

• Single-link distributed mass systems 

• Direct drive motors 

• Planar systems 

• Modelling based on truncated mode shapes 

Our experimental target: 

• General multi-link, 3-D system 

• Mechanical actuation, with friction, backlash, 
etc. 

• Possibly joint-dominated in mass 
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Experimental Apparatus at Carnegie-Mellon 


Reconfigurable manipulator; modular design, up to 6DOF 

• Mechanical 

• Each joint: DC-motor/harmonic drive/ 
potentiometer 

• Reconfiguration and link changeout 
using tubing and NPT hardware 

• Computational: 

•Motorola VMEbus System 1000: 
68010, VersaDOS, Pascal 

• Smalltalk-80 (One concurrent 
implementation under VersaDOS, a 
second Unix implementation on Sun-3) 

• Sensing: 

• Position (rotation) sensing on joints 

• Strain sensing on links 

• Vision end-point sensing (not used for 
control) 
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MANIPULATOR IN A 3-DOF TWO LINK CONFIGURATION 



MANIPULATOR IN A 4-DOF TWO LINK CONFIGURATION FOR FORCE COGNITIVE EXCAVATION 
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A Simple Model for a General Flexible Link 


Starting points to consider: 

• Link motion results from concatenation to 
other links; the non-linear configuration 
problem, present in rigid manipulators as well. 

• The link itself deforms as a result of the end- 
forces and the inertial forces acting on it. 

• Which quantities can be observed or sensed? 

• Which quantities can be controlled? 

• How is the (approximate) solution to the PDE 
for link deformation to be contained within the 
dynamics equation? (What are the amplitude 
variables for the generalized co-ordinates 
chosen?) 


First step in the approach: View first the motion of the 
chord connecting the end points, and then refer the 
(elastic) deformations to that chord. 



Kinematics/Mechanics of the Simple Model 

Chord motion: 

• Denote rotation by 0, equivalent to a rigid-link 
formulation. 

• Include dynamic effects of concentrated 
masses and inertias at joints. 

• Assume that inertial effects of the link are 
modelled (from 0 and co) by the translation 
and rotation of the (c.m. of) the chord. 


Deformations (displacements) of the flexible link: 

• Displacements y(x) are referenced to the 
chord. 

• Assume that the displacements equal those 
resulting from static application of end- 
rotations <|> and \y. 

• Displacement and potential energy: 

y(x) = xCx-l>[CfU0x - 

U = * t\ j/£ 
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A SINGLE LINK SYSTEM; EQUATIONS OF MOTION FROM APPLICATION OF THE SIMPLE MODEL 
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COMPARISON OF TRANSFER FUNCTIONS, SIMPLt 
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Properties of the Simple Model 


The model is equivalent to Rayleigh’s method, using 
assumed shape with two amplitude variables, <)> and \|/. 

• Inertial effects of the joints are properly 
modelled, and are consistent with the 
mathematical formulation. (Functions in 0) 

• The model is also a lumped mass assumption 
of m acting on the chord. (If this is the 
dominant link inertial effect, then the error is 
small.) 

• The assumed shape has only 2 "dof," and can 
only approximate the real shape. 

• Some higher order effects are plainly 
"missed," as they would be for a truncated 
mode solution. 

• The formulation would be useful for control, 
because the variables 0, <|> and \y can be 
measured and actuated. 

• For joint-dominant systems the model should 
be very accurate, and if joint inertias are small 
the equations reduce in order. 
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Applications to Manipulator Control 


Equations of motion can be used as follows: 

• To confirm the number and the identity of 
state variables for control. 

• To perform simulation studies. 

• To set gains from classical control theory. (In 
principle) 

• To compute variable gains for a non-linear 
system. (In principle) 

• To accomplish model-based (shaped) control. 

• To accomplish model-based feedforward 
control; requires real-time performance. 



Single Link Systems 


1. Planar (2-D) motion 

• Actuator was a direct drive DC motor. 

• The simple model produces 3x3 equations of 
motion. 

• Tip has mass but low inertia; system order 
reduces to 2; state variables are identified as 
0 and <|>. 

• Sensing of rotation (0+<|>) and strain (~<)>). 

• Perform experiments; set gains by trial and 
error. 

• Discussion of friction effects. 


2. Spherical (3-D) motion 

• Actuation using two joints of the modular 
manipulator. 

• See videotaped results. 

• Controllable despite friction and backlash. 
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Planar Two Link System 


• Experiments performed at Toshiba. 

• Air table, 2-D manipulator. 

• Four state variables: 0’s and <t>’s. 

• Compare experimental and simulation results. 

• Friction in actuators causes vibration. 

• Feedforward control is attempted, inclusive of 
friction effects. 

• Model based feedforward control limits 
vibration. 
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EXPERIMENTAL RESULTS, FEEDBACK CONTROL: NOTE VIBRATION EFFECT, FROM FRICTION 
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EXPERIMENTAL RESULTS, FEEDFORWARD CONTROL: NOTE ELIMINATION OF VIBRATION 
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Combined Flexural-Torsional (3-D) Motion 


Experiments performed April 1988: 

• Three actuated DOF (yaw, pitch, roll). 

• Two links; one flexible, one rigid. 

• Linear feedback control; gains by trial and 
error. 

• Coupled flexural and torsional vibrations. 

• See videotape; see experimental results. 

• Actuator properties by system identification 
(in process). 

• Analtyical model used in simulation studies (in 
process). 

• Next phase: distal link made flexible, 4 (or 5) 
actuated DOF. 
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Roll strain Pitch strain Yaw strain Roll angle, radians Pitch angle, radians Yaw angle, radians 






COUPLED FLEXURAL-TORSIONAL MOTION ROTATIONS AND STRAINS 
FEEDBACK CONTROL ON POSITION; NO FEEDBACK ON STRAIN 
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COUPLED FLEXURAL-TORSIONAL MOTION: STRAINS 
FEEDBACK CONTROL ON POSITION; NO FEEDBACK ON STRAIN 
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COUPLED FLEXURAL-TORSIONAL MOTION: STRAINS 
FEEDBACK CONTROL ON POSITION AND ON STRAIN 
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Discussion and Conclusions 


• The model may be well suited for serial link 
manipulators, including joint-dominated 
systems. 

• Accuracy for MDOF systems, non-linear in 
configuration, remains to be examined. 

• Control experiments must be extended 
beyond linearized regions. 

• The major application is model-based control, 
still to be studied in depth. 

• Effects of friction, backlash, deadband should 
be included. 

• Friction or torque ripple can excite higher 
modes. 

• Frequencies of unmodelled (higher) modes 
can be adjusted by inserting "redundant" 
actuators. 
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MINIMUM-VARIANCE REDUCED-ORDER ESTIMATION ALGORITHMS PROM 
PONTRYGDTS MINIMUM PRINCIPLE 

By 

Yaghoob S. Ebrahlml 
Boeing Commercial Airplanes 
Seattle, Washington 


ABSTRACT 


It has become apparent with the introduction of modem control and estimation theory 
(hat the entire knowledge of a system cannot be included in the system design for most 
practical applications. Such mechanization of a total system usually results In a model 
exceeding the capacity of a real-time processor, thus requiring a reduction in the state 
size. In addition, linear Altering and smoothing problems have been extensively 
Investigated for a case wherein the Alter or smoother is of the same state dimension as 
(he dimension of the "best state" model available. 

This paper presents a uniform derivation of minimum-variance reduced-order (MVRO) 
filter-smoother algorithms from Pontrygln's Minimum Principle. An appropriate 
performance Index for a general class of reduced order estimation problem Is 
formulated herein to yield optimal results over the entire time interval of estimation. 
These results provide quantitative criteria for measuring the performance of certain 
classes of heurlstically designed, suboptimal reduced-order estimators as well as 
explicit guidance to the suboptimal Alter design process with both continuous and 
discrete Alter-smoolher algorithms being considered. 

By the duality principle, the algorithms of reduced-order estimation can be easily 
extended to the deterministic problems of optimal control (l.e.. the regulator and linear 
tracking problem). 
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MODIFYING HIGH-ORDER AEROELASTIC MATH MODEL OF A JET TRANSPORT 
USING MAXIMUM LIKELIHOOD ESTIMATION 


By 

Amir A. Anissipour and Russell A. Benson 
The Boeing Company 
Seattle, Washington 


ABSTRACT 


The design of control laws to damp flexible structural modes requires accurate math 
models. Unlike the design of control laws for rigid body motion (e.g., where robust 
control is used to compensate for modeling inaccuracies), structural mode damping 
usually employs narrow band notch Alters. In order to obtain the required accuracy In 
the math model, maximum likelihood estimation technique Is employed to improve 
the accuracy of the math model using flight data. This paper presents all phases of this 
methodology: (1) pre-flight analysis (i.e.. optimal Input signal design for flight test, 
sensor location determination, model reduction technique, etc.). (2) data collection and 
preprocessing, and (3) post-flight analysis (i.e., estimation technique and model 
verification). In addition, a discussion Is presented herein of the software tools used for 
this study and the need for future study In this field. 
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Modifying High-Order Aeroelastic Math Model of a Jet Transport 
Using Maximum Likelihood Estimation 

| 

Amir A. Anissipouf 
Russell A. Benson 


The Boeing Company 
Boeing Commercial Airplanes 
P.O.Box 3707 M/S9W-38 
Seattle ,Wa. 98124 


ABSTRACT 

The design of control laws to damp flexible structure modes requires accurate math 
models of the dynamic system. To obtain the required accuracy of a math model, 
the parameter estimation technique using maximum likelihood estimation is 
employed to improve the accuracy of the model based on flight data. This paper 
presents all phases of this methodology: pre-flight analysis (i.e., optimal input signal 
design for flight test, sensor location determination, model reduction technique, 
etc.), data collection and preprocessing, and post-flight analysis (i.e., estimation 
technique and model verification). The results of this study indicate that the 
parameter estimation technique (i.e., maximum likelihood estimation) is an 
effective and powerful technique in modifing high-order aeroelastic aircraft models. 
However, the accuracy of the results depends upon the fidelity of the theoretical 
model with regards to the correct number of dominant modes for the desired 
frequency bandwith in the model (i.e., model order). If the number of modes in the 
model are not representative, then an identification problem can occure in the 
parameter estimation technique. Nevertheless, this problem can be overcome using 
the system identification technique. 
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INTRODUCTION 


Having an accurate mathematical representation is fundamental to any aircraft 
control system design. In general, aircraft models are developed from a theoretical 
basis and modified by analyzing the experimental data (i.e., wind-tunnel data for 
aerodynamic models or ground shake test data for structural models). Although 
present techniques provide very good dynamic models for the design stages of an 
aircraft, often these models do not match the actual dynamic flight response. This 
problem has generated a need for advanced system identifiaction and parameter 
estimation techniques in upgrading dynamic models of an aircraft based on flight 
test data. This modeling problem is more apparent with high-order aeroelastic 
models with which our experince with modeling techniques is limited. 

Low-frequency structural modes are easily excited for a jet transport with a long 
fuselage. This excitation causes a lateral ride discomfort in certain flight conditions. 
In order to design a yaw damper to dampen Dutch roll response and suppress the 
undesirable low-frequency structure modes by means of active control, an accurate 
aeroelastic model of the aircraft must be available. In this study, parameter 
estimation technique is applied to upgrade the high-order aeroelstic math model of 
a jet transport. The following is a summary of the parameter estimation technique 
using maximum likelihood estimation. 


Maximum Likelihhood Estimation 

Suppose the actual system is described by (Reference 1): 

x(t) = A x (t) + B u(t) + S s(t)+F n (t) 
z(t. ) = C x(t.) + Du(t.)+Hs(t.) + G m(t.) 


where 



x (t) 

state vector 

u(t) 

control vector 

z(t,) 

measurement vector 

s(t) 

bias vector 

n(t) 

process noise 

m(t,) 

measurement noise 

ti 

time sample 

A ,B ,C ,D,S ,H,F ,G 

system matrices with unknown parameters 


n(t)and m(t) are zero mean .Gaussian and independent noise 

Assume k is the vector of unknowns that contains elements of the system matrices 
A, B, C, D, S, H, F and G. The objective is to maximize the probability distribution of 
unknowns (i.e., k) when the measurements z are available. Therefore, maximizing 
P(k/z), where P is the probability distribution function of k given z. 


( 2 ) 


(3) 


By Bayes’ rule: 

P(k/z) P(z) = P(k, z) = P (z / k) P(k) 


or 


P(k/z) = P (z / k) 


POO 

P(z) 


Since in these equations z is given, so P(z) becomes a constant. Assume there is no a 
priori preference for k, so P(k) becomes a constant. Therefore, P(z/k) differs from 
P(k/z) only by a constant. In other words equation (3) becomes: 

P(k/z) = P(z/k) -constant (4) 


Equation (4) indicates that P(z/k) may be maximized instead of P(k/ z). Therefore, 
using Gaussian assumption, the likelihood ratio may be written as: 


-In 


N 


P(z / k) = [(27i) L |GG*| ] expj- -j £ [z k (t 1 )-z(t l )J* (GG*)” 1 [ z k (t 4 ) -z(t t )] 
where 




z k (ti) predicted estimate at time ti 

GG* measurement noise covariance matrix 

L number of measurements 


( 5 ) 
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If the logarithm of equation (5) is taken, the consatnt terms are eliminated by the 
maximization, and the equation is multiplied by -1 to do minimization rather than 
maximization, then equation (6) will be obtained as: 

J(k)=i l{p k (t | )-z(t j )r(GGr , p k (t.)-z(t ) J]} 

i = 1 (6) 


where J(k) is the cost function to be minimized. Two steps are taken to obtain Zk(tj). 
Prediction step: 

x k< t i + i> = 0 x k<ti> + ^ + 

~ X k(*i + l) + ^ U (*i + l) 

where 

-t 

O = e A4t and *¥ = J e A * ds 

J o 

and the correction step: 

x k^i + i^ = x k(*i + i) + ^ [z (t i+1 ) - z k (t i + 1 )] 

K in equation (8) is the Kalman filter gain matrix given by: 

K = P C * (GG* )~ 1 


where P is the solution to the discrete time Riccati equation: 


AP+PA* --^PC*(GG*) -1 CP + FF* = 0 


( 10 ) 


After obtaining the cost function J(k), the Newton-Raphson algorithm is used 
iteratively to minimize the cost function by revising the unknowns parameters. 


k| + i * k. - {V k J (k,)} - ' {V k I(k,)} 


( 11 ) 
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This algorithm requires an intial estimate for the vector of unknowns (ko). A priori 
estimate is available for each unknown parameter through the analytical model. 

The MMLE software tool developed by NASA Dryden is a parameter estimation 
program supporting this estimation technique. This software has been modified by 
Boeing to accept and handle higher order models. A comprehensive description of 
this software tool is described in Reference 1. 


ERE -FLIGHT ANALYSIS 


Math Model 

A sixtieth order linear aeroelastic math model for a flight condition of Mach .6 
speed, 15000 foot altitude, and no turbulance, and cruise configuration of a jet 
transport was provided in the form of: 

Mq+Cq+Kq=u 


where 


M mass matrix 

C damping matrix 

K stiffness matrix 

q generalized coordinate 

u control inputs 


The model is defined in the inertial axis system, and the dynamics (q), consist of 
rigid body and flexible modes. The model is timed using data from ground shake 
testing. The system of equations (12) was transformed into state-space form using 
the following transformation: 
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therefore the system equation (12) becomes: 


x = A m x m +B m u 
y = C m x m +D m u 


where 



- M" 1 C 


This transformation always exists because the mass matrix is positive definite. 
Although this is a well-posed theoretical problem, it is not trivial. The flexible 
model is usually on the order of one hundred states, thus causing numerical 
inaccuracies in the inversion of the mass matrix. In our analysis the software 
package MPAC was used to perform the transformation. (MPAC is a numerically 
robust modem control and analysis software tool developed by the Boeing 
Company.) 

For the identification process, the system equation (13) was transformed into the 
conjugate modal form using the following transformation: 

m = T“ 1 x m 


Equation (13) becomes: 

m = Am + Bu 
y =Cm + Du 

where 
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A = dia (X t ) Xj = i* eigenvalue 

B = T” 1 B controllability matrix 

C =CT observability matrix 

The advantage of using the modalized form given by equation (14) is that all the 
modes through A matrix, along with the controllability and observability matrices 
are readily available for an analyst to quickly locate uncontrollable and unobservable 
modes. In addition, the modes in the A matrix are decoupled and may be 
partitioned into rigid model and elaastic model. 

The order of the model was reduced to nineteen by deleting the modes above 6 Hz. 
Since this model will eventually be used for ride quality study and modal 
supperasion design, only those modes less than 6 Hz were retained. 

The reduced order, modal model (19th order) is represented by: 

m r = A r m r + B r u 

y =C r m r + D r u (15) 

This model contains one state for heading, one for the spiral mode, two for the 
Dutch roll mode, one for roll mode, eight for low-damped elastic modes, and six for 
high-damped elastic modes. 

To support this study, a special set of sensors were installed on the aircraft to 
measure the dynamic response of the jet transport. The locations of these sensors 
were based on the mode shapes of the aircraft determined by the math model and 
physical constraints (Table I). (A complete discussion on sensor selection and 
location placement on the aircraft is omitted herein for proprietary reasons.) 
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TABLE I: Sensor Type and Locations for High-order 
Aeroelastic Modeling 


SENSOR TYPE 

SENSOR LOCATION 

Position Transducer 

On all control surfaces 

Yaw Rate Gyro 

Pilot seat, IRU (a station between CG and 
cockpit below the cabin floor), CG station 

Lateral Accelerometer 

1 Pilot seat, 1 Cockpit ceiling, 8 on the 
passanger cabin floor from the cockpit to the 
aft galley, 1 on the aft galley ceiling, 3 on 

vertical tail (tip and mid section, front 
and rear spar), three on each nacelle, 

1 IRU station 

Vertical Accelerometer 

1 on the pilot seat, 1 IRU, 1 aft galley, 

8 on each wing, 3 on each horizontal tail, 

2 on each nacelle 

Roll rate, Yaw rate, 
Bank angle, Heading, 

IRU and CG stations 


The sensors selected for the analysis were: body roll angle (0), heading angle OP), roll 
rate (p) and yaw rate (r) at the ERU; body yaw rate at the pilot seat; 9 lateral 
accelerometers along the fuselage; 2 lateral accelerometers on the nacelle number 2; 
and 3 lateral accelerometers on the vertical tail. 


Input Signal Desig n 

The flight test input-signal design analysis for high-order aeroelastic modeling was 
performed using the reduced order analytical model (equation 15). Although a 
number of "optimum" input signals have been proposed for flight testing in 
conjunction with parameter estimation, none have been found to be appropriate for 
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high-order aeroelastic modeling. Essentially, all the analytical techniques proposed 
in designing the optimum input signals are based on the analytical model. This 
model is the subject of improvement by the identification and estimation 
techniques. Hence, no "optimum" input signal exists. 

A number of different input signals were evaluated for this study. After a 
comprehensive simulation study, it was determined that a frequency sweep of a 
linear sine-wave with adequate energy to excite all the modes (rigid and elastic) 
yeilds the best results. In addition, the linear sine-wave frequency sweep optimizes 
the most commonly used criterion for input signal design: 

91 = -log (detM) (i6) 

where M is the Fisher information matrix (or sensitivity matrix) defined by: 

M = V’ k J(k> (17) 

J is the cost function defined in equation (6). The criterion 91 defined in equation 
(16) is related to the volume of highest probability density region for the parameters 
k. An interesting property of the determinant criterion is that it is independent of 
scaling parameters (Refemce 2). 

Fifteen tests were designed for the same flight condition. Five frequency sweeps 
were designed for each control surface. Each test was repeated for rudder, aileron, 
and both surfaces in phase. The first frequency sweep covered 0 to 6 Hz to excite all 
the modes in one test. The other four tests were then designed to excite specifically 
high-damped modes by sweeping from .25 Hz below to .25 Hz above the frequency of 
the mode. 

The amplitude of the input signals were designed to be constant for practical 
purposes (i.e., rate limits). The designed input signals were tested in the lab to 
confirm that the signals did not saturate the servos and actuators of the control 
surfaces. However, the output of the actuators during flight test generated signals 
with decaying amplitutes. These decaying amplitudes reduced the energy level 
initially designed for the test. Figures 1 and 2 show the actual control surface 
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Figure 1. Actual Control Surface Deflections for Flight Condition 21 
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Figure 2. Actual Control Surface Deflections for Flight Condition 41 




deflections for rudder sweep alone, and for rudder and aileron surfaces 
simultaneously in phase. 


Sampling Frequency 

To record the data in flight test, a simulation study was conducted to determine the 
required sampling frequency. The analytical model (i.e., system equations 15) was 
assumed to be the true model, and simulated using the designed input signal. A 
considerable amount of noise was added to the simulation data, and then that data 
was treated as pseudo-flight data. The acutal model was used for parameter 
estimation to determine the required sampling frequency. Sampling frequencies of 
20, 25, 50, 100, 200 Hz were considered for this study. One mode or group of modes 
at a time were selected for the estimation process of each sampling frequency. The 
results indicated that 100 Hz is the best sampling frequency for this study. Figure 3 
shows the typical results for identified parameters when different sampling 
frequencies were used. 



Figure 3. Typical Results from Estimation with Different 
Sampling Frequency 
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FLIGHT TEST 


The flight test was performed using designed linear sine-wave frequency sweeps for 
rudder and aileron. The test conditions were conducted at a speed of Mach .6, an 
altitude of 15,000 feet, and minimal turbulence. A preprogrammed frequency 
function generator was used to apply the linear sinusoidal frequency sweeps (0-6 Hz) 
to the aileron and rudder (through the autopilot servo). 

The flight test data were recorded with 100 sample per second, and then filtered 
using a Graham low-pass filter with the cutoff frequency of 10 Hz and rolloff 
frequency of 15 Hz. Prior to estimation analysis, the data were cleaned up by 
removing all the sensor biases and data dropouts. 


POST FLIGHT ANALYSIS 

The analytical model (system equations 15) was simulated using actual control 
surface defelection during flight as input signals. The comparison of flight data with 
the response of the analytical model for flight condition 41, where both rudder and 
aileron frequency sweeps are used, is presented in the Figures 4-11. 

The maximum likelihood estimation software tool (MMLE) developed by NASA 
Dryden was used to minimize the residuals between flight data and response of the 
analytical model in Figures 4-11. At the time of analysis, MMLE was hosted on the 
Cyber mainfram. Due to Cyber having a memory limit, the capability of using 
process noise was not available for analysis. Hence the results obtained herein, are 
preliminary results which do not include the effect of process noise. The final 
results of this study will be reported at the 1989 AIAA Guidance, Navigation and 
Control conference. 

The high-order model was partitioned into two sections: rigid model and elastic 
model. For rigid model identification, 15 seconds of data were used. First the rigid 
portion of the control and measurement matrices were upgraded. Then, the A 
matrix was upgraded. Finally, all the parameters in the rigid section of the A /B and 
C matrices were simultaneously estimated. 
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Figure 4. Time Response of Body Roll at IRU Comparing Flight Data with Math Model 
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Flight condition 41, Rudder and Allwon Input 
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Figure 5. Time Response of Body Yaw at IRU Comparing Flight Data with Math Model 
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Two different approches were taken for the elastic model identification. In the first 
approach, the 19 th order model was used for the analysis with all the elements of B 
and C being estimated. All 70 seconds of data were used for this estimation 
approach, . In this process, those parameters in the B andC that did not contribute 
to the residuals were identified and kept constant for the remainder of the analysis. 
Then, the elements of A were added to the estimation process while keeping some 
of the elements of B andC constant. The results of this estimation approach are 
show in Figures 12-19. 

The second approach was to add one elastic mode at a time to the rigid model. For 
this approach, the first elastic mode was added with 28 seconds of data used for the 
analysis. The corresponding parameters in the B and C matrices were estimated 
every time a mode was added to the model. The result of this approach was not 
satisfactory because several times the algorithem diverged and the residuals were 
big. 

Figures 20-26 show the PSD plots obtained from the analytical model. Figures 27-34 
show the PSD plots obtained from the estimated model. The PSD plots obtained 
from the estimated model, clearly show that the estimation analysis improved the 
accuracy of the model in terms of its modal representation. However, the estimated 
parameters in the B and C matrices are biased. Since an accurate representation of 
the transfer functions was desired for this study rather than true values of the B 
and C matrices, the biased estimates in theB and C matrices did not create any 
problem. 

Figures 16, 17 and 19 indicate that another mode is present in the flight data which is 
not modeled in the analytical or estimated model. This problem can not be solved 
via parameter estimation technique which assumes the structure of the model (i.e., 
the order of the model) is correct. Hence, it is suggested that the system 
identification technique developed by V. klein and J. Batterson of NASA LaRC be 
used to overcome this problem. 
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Figure 12. Time Response of Body Roll at IRU Comparing Flight Data with MLE Model 
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Figure 13. Time Response of Body Yaw at IRU Comparing Flight Data with MLE Model 



Z.#T 


\M 




. 5 - 

Roll Rate 
IRU 

( deg / sec ) 

• . 0 " 


o 

CX> -i.i- 


- 1 . 5 - 


-2.0- 


Fi 




609 


Flight condition 41, Rudder and Aileron Input 
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Figure 15. Time Response of Yaw Rate at IRU Comparing Flight Data with MLE Model 
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Figure 20. Power Spectral Density Response of Body Roll at IRU Comparing Flight Data 
with Math Model 
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Figure 23. 


Frequency, HZ 

Power Spectral Density Response of Yaw Rate at IRU Comparing Flight Data 
with Math Model 
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igure 25. Power Spectral Density Response of Lateral Acceleration at IRU Comparing 
Flight Data with Math Model 
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Power Spectral Density Response of Roll Rate at IRU Comparing Flight Data 
with MLE Model 
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Figure 32. Power Spectral Density Response of Lateral Acceleration at IRU 
Comparing Flight Data with MLE Model 
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Figure 33. Power Spectral Density Response of Lateral Acceleration at AFT Body 
Comparing Flight Data with MLE Model 
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Comparing Flight Data with MLE Model 
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AUTOMATED MODEL FORMULATION FOR TIME-VARYING FLEXIBLE STRUCTURES 

By 

B. J. Glass* and S. Hanagud 
Georgia Institute of Technology 
Atlanta, Georgia 

ABSTRACT 


The control of many types of flexible structures, such as robotic manipulators or large 
space structures, usually requires an accurate analytical model. Once obtained, these 
models are currently compared with observations of the behavior of the structure and 
incremental changes can be made by using Kalman filtering techniques or other 
parameter identification techniques. For time-varying flexible structures, however, 
such changes may occur In sudden changes to boundary conditions or to the form of the 
model differential equations. Some of the primary causes for such changes are growth, 
reconfiguration or damage. This class of changes often requires a reformulation of the 
analytical model. This paper presents an identification technique that uses the sensor 
information to choose a new model out of a finite set of discrete model space, in order to 
follow the observed changes to the given time varying flexible structure. Boundary 
condition sets or other information on model variations are used to organize the set of 
possible models laterally into a search tree with levels of abstraction used to order the 
models vertically within branches. An object-oriented programming approach is used 
to represent the model set In the search tree. A modified A* best first search algorithm 
finds the model where the model response best matches the current observations. 
Several extensions to this methodology will be discussed. Methods of possible 
Integration of rules with the current search algorithm will be considered to give weight 
to Interpreted trends that may be found In a series of observations. This capability 
might lead, for instance, to identifying a model that Incorporates a progressive damage 
rather than with Incorrect parameters such as added mass. Another new direction Is to 
consider the use of noisy time domain sensor feedback rather than frequency domain 
information in the search algorithm to improve the real-time capability of the 
developed procedure. The next logical step will be to automatically expand the model 
space by adding subsets of recognized possible models. This can be accomplished by 
using developed methods of machine learning. Finally, testing of the currently 
developed approach with model spaces that are derived from more complex structures 
will be discussed. 

* Dr. Glass is currently with the System Autonomy Demonstration Project Office of the 
NASA Ames Research Center, Moffett Field, California 94035 
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NUMERICALLY EFFICIENT ALGORITHM FOR MODEL DEVELOPMENT 
OF HIGH-ORDER SYSTEMS 

By 

L. O. Parada 

Calspan Advanced Technology Center 
Buffalo, New York 


ABSTRACT 


Frequency domain parameter identification techniques provide a straightforward 
approach to transfer function estimation. However, for high-order systems, numerical 
difficulties may be encountered during the estimation process. Inaccuracies may result 
because of the large variation of the transfer function polynomial coefficients for high- 
order systems. The lack of numerical precision to represent this variation may cause 
the estimation process to break down. 

This paper presents a technique for estimating transfer functions in partial fraction 
expansion form from frequency response data for a high-order system. The problem 
formulation avoids many of the numerical difficulties associated with high-order 
polynomials and has the advantage of having the option to fix the damping and 
frequency of a mode. If known, during the estimation process. The resulting transfer 
function(s) may be converted to Jordan-Form time domain equations directly. 

During the implementation of this technique, a frequency and amplitude normalizing 
window was developed that maximized the efficiency of the optimization algorithm. 
The combination of estimating the transfer function in factored form, the ability to fix 
previously determined parameters and the effectiveness of the normalizing window led 
to a progressive approach to synthesizing transfer functions from frequency response 
data for high-order systems. 
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Abstract 


Numerically Efficient Algorithm for Model Development of High Order Systems 


L. 0. Parada 

Cal span Advanced Technology Center 
P.0. Box 400 
Buffalo, NY 14225 
(716) 632-7500 


for presentation at the 

NASA Langley Research Center Workshop on 
Computational Aspects in the Control of Flexible Structures 


Frequency domain parameter identification techniques provide a straightforward 
approach to transfer function estimation. However, for high order systems, 
numerical difficulties may be encountered during the estimation process. 
Inaccuracies may result because of the large variation of the transfer function 
polynomial coefficients for high order systems. The lack of numerical precision 
to represent this variation may cause the estimation process to break down. 

This paper presents a technique for estimating transfer functions in partial 
fraction expansion form from frequency response data for a high order system. 
The problem formulation avoids many of the numerical difficulties associated 
with high order polynomials and has the advantage of having the option 
to fix the damping and frequency of a mode, if known, during the estimation 
process. The resulting transfer function(s) may be converted to Jordan-Form 
time domain equations directly. 

During the implementation of this technique, a frequency and ampl itude normal izing 
window was developed that maximized the efficiency of the optimization algorithm. 
The combination of estimating the transfer function in factored form, the 
ability to fix previously determined parameters and the effectiveness of 
the normalizing window led to a progressive approach to synthesizing transfer 
functions from frequency response data for high order systems. 
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NUMERICALLY EFFICIENT ALGORITHM FOR 
MODEL DEVELOPMENT OF HIGH ORDER SYSTEMS 

Statement of Problem 
Development of Mathematical Models: 

Time Domain - Difficult to implement 

- Instrumentation Complement 

- Input Design 

- Noise 

- Computational Load 
Freq Domain - Simplified Implementation 

- Fewer parameters per computation 
cycle 

- Statistical methods applicable 
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PREVIOUS WORK 


Frequency domain parameter identification requires 

Determination of characteristic equation 
(nonlinear or iterative techniques) 

Estimation of numerator polynomials 

Factor characteristic equation 

Estimate zeros or residues 


Inaccuracies (for high order systems) due to: 

Variation of transfer function polynomial coefficients 
Transformation errors 

Sensitivity of polynomial roots to variations in polynomial coefficients 



638 


SUMMARY OF CONTRIBUTIONS 


Development of technique to estimate transfer functions in partial fraction 
expansion form from frequency response (amplitude and phase) data 

Elimination of numerical difficulties associated with high order polynomials 

Incorporation of a priori knowledge of system modes (frequency and 
damping) directly into the estimation process 

Development of frequency and amplitude normalizing window that maxi- 
mizes effectiveness of the optimization algorithm and eliminates the initial 
guess problem 

Stepwise approach for synthesizing transfer functions where order of sys- 
tem is high and unknown 
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FACTORED FORM ESTIMATION 


Classical Nonlinear Regression Problem 

Estimate parameters from measured amplitude and phase data 

Error Function: 

Square of distance between measured and estimated frequency responses summed over all discrete frequency 
points 

M 2 

e - E [FflUj) - G (l“i ) ] 
i = 1 


where: M = # frequency points 

F(j co) = measured frequency response 
GQ w) = estimated frequency response 
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where: N = order of system 

Q = # of second order terms 
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Express measured and estimated frequency responses in terms of real and imaginary 
components 
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Solve for unknown parameters 

Set partial derivatives equal to zero 

Solve using nonlinear optimization technique 


Fifth Order Single Precision Example 

Simulate parameter identification of high order system 
Modes distributed over wide frequency range 

Single precision: Scale down problem 

Reduce number of variables 

5th Order Transfer Function: 

Cascade form 

(S + 5 x 10~ 3 ) (s + 5 x IQ -1 ) (S + 5 x 10^) (3 + 5 x 10 +3 > 

(s 2 + 2 x 10” 3 s +10 -4 ) (8 + 1) (s 2 + 1 x 10 +4 S+ 10 +8 ) 

Parallel form 

1,253955 X 10~% + 6.122844 x 10~ 6 + 1,221073 X 1(T 3 ( 9.975250 x Ifl^S + 5.024754 x 10+ 3 
S 2 + 2 x 10~ 3 s + 1 X 10” 4 s + 1 s 2 + 1 x 10 +4 s + 10 +8 

Frequency Range: 1 x 10” 4 to 1 x 10 +5 Hz. 

# Points/Decade ■ 30 
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DENOMINATOR COEFFICIENTS 


Term 



Exact Coefficient 


1.000000 


10001.001 2 


100010012.0021 


100200211.0001 


210001.0 

100004.0 


Additive Components 
1.0 

1.002 + 1.0 x let 4 

2.1 x 10 -3 + 1.002 x 10 +4 + 10 +8 
1 x 10' 4 + 2.1 x 10 +1 + 1.002 x 10 +8 
1.0 + 2.1 x 10 +5 
1.0 x 10 +4 


= Single Precision Variable Representation 
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LINEARIZED APPROACH 

Initial Error Function: 

P0“ k ) 

E = FQ'o) ) - v 

k k Q (ja) k ) 

where: FQo> k ) = measured frequency response at w k 
PG w k ) = estimated numerator polynomial at w k 
QG w k ) = estimated denominator polynomial at<»> k 

Weighted Error Function: 

E k =E k Q(i “k> = F (i“ k ) Q0“ k ) -P0“ k ) 

Iterative Error Function: 



[A] [X] = [B] 

Solve for parameter vector [x] 

2 

Iterations converge to minimization of | e k I 



5th Order example: Polynomial Results 


Exact 

Transfer 

Function 

Linear 

Results 


Cost 

Function 


(St5x 10~ 3 ) (s + 5 x IQ' 1 ) (S + 5 x 10 +1 ) (S + 5 x 10 +3 ) 

(s 2 + 2 x 10 -3 s + 10 - 4)(s + 1) (s 2 + 1 x 10 +4 s + 10 +8 ) 

(S + 5 x 10~ 3 ) (s + 5 x IQ" 1 ) (S + 4.24 x IQ 4-1 ) (s - 7.17 x 10 +3 ) 

(S 2 + 2 x 10" 3 S + 10 _4 ) (S + 9.98 x 10 _1 ) (S - 7.17 x 10 +3 ) (s + 1.72s 10 +4 ) 

2.37 x 10" 10 



270.0 n 


180.0-4 



- 1800 [ ■ i 1 rr !ff ^ 1 1 TH Tff n 1 rrm w ^TTm^ r rrm w [ ' mrn^ r vrvm[ n n rw; ' 


1 X 10 


1 X 10 


-2 


6A 


1X10 W 1 X 10 4 

FREQ-RAD/SEC 


1 X 10 


1 X 10 
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FACTORED FORM APPROACH 


Nonlinear method of solution - strong initial guess required 

Initial Guess - Examine bode plots to approximate frequency and damping 

of modes 

Problem - Error function relatively insensitive to perturbations in para- 
meters of high frequency modes 

- Gradient expressions small compared to those of the lower 
frequency parameters 

Solution - Normalizing window 

Scale Data Such That: 

- Frequency is centered around 1 rad. 

Amplitude is centered around 1 (abs. units) 

Benefits - Need for strong initial guess eliminated 

effectiveness of optimization algorithm maximized 
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STEPWISE FACTORED FORM TECHNIQUE 
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Exact 

Transfer 

Function 

Factored 

Form 

Results 


Cost 

Function 
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i 
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OUT/IN PSI-DEG 




OUT/IN AMP-ABS 


16th Order Transfer Function Estimation 
Pcg/Vg: Roll Rate Measured at C.G. of Aircraft 

vs. 


1 X 10 2 -a 

1 X to 1 - d 

1X10°-' 

- 

1 X 10 -1 
1 X 10 " 2 

1 X 10 “ 3 -| 

1 X 10" 4 __, 
1 X 10 


Unit Gust Along Y-Body Axis 


+ + + = measured 



FREQ-RAD/SEC 



-270.0 4 


-360.0 


"rnrrmTn — *-r 


i * r r mn — 1 i 1 i , m , i , i , n — t 

1X10* 1 1 X 10° 

FREQ-RAD/SEC 


TTTTTITf > I 
1 X 10 1 


1 X 10 


-3 


1 X 10 


-2 


FTTriTJ 

1 X 10 2 


TTTTTr) 

1 X 10 2 


Cost Function: 4.5 x 10" 12 
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CONCLUSIONS 


Development of technique to estimate transfer functions directly in factored form 
Advantages: 

Ability to fix damping and frequency of a mode, if known, during the esti- 
mation process 

Avoidance of numerical difficulties associated with high order polynomials 
Ability to obtain Jordan-form time domain equations directly 

Progressive approach to transfer function estimation through use of a frequency 
and amplitude normalizing window 

Development of frequency and amplitude normalizing window that eliminates the 
initial guess problem and maximizes the effectiveness of the optimization algo- 
rithm 
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Amplitude 


SCOLE DAMPING 










In One Dimension 


Without hysteresis 

x(t ) + co 2 x(t) + 2(0 C,x(t) 

+ y x(t) 2 m |x(r)| a x(tf n + 1 WOI 
+ Bu(t) + FN(t) 

= 0 

0 < a, p < 1 . 


GSS- 






AMPLITUDE 


NONLINEAR DAMPING: £ = 0 



NUMBER OF CYCLES 
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Beam Model 


u(s, t) + Xu”" {s, t) - 2t^X ii"(s, t) 


f ~ L 


- Y 


f u\s,t)u\s,i ) ds I 
o 


2 (n + p)+ 1 


u"{s, t) 


= 0 , 


0 < s < L; 0 < t 


n : zero or positive integer 

0sp<i 

C , : Linear Damping Ratio 
Prime represents space derivative 
Dot represents time derivative 
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A ~ 


d. 

: clamped beam 


9k 

'l A (j)* = (tik$k 
~ (- 1 ) 


* 


ds 2 


Lj 

J u'(s, t) u'(s, t) ds 


0 


= J [«(J, t) u\s, f)] 

L 

- J u(s, t) u'\s, t ) ds 
0 
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x(t ) = u(‘, t) 


— > ~ [x(0, Va x(0] 


Dx) = y([x, Va x]) 


2(a + ( 3 )+ 1 


+ XAx(t) + DJc(t) 

+ F(x(t), Dx{t)) + Bu(t) 

= 0 

£> = 2t,^K<A 
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x(t ) = u(-, t) 


Mx(t ) + XAx(t) + Dx{t) 

+ F{x{t), Dx(t )) + Bu{t) = 0 


Energy 

E(t) = j {[x(r), i(f)] + X[Ax(t), x(0]} 

4 jE(t ) = [* + XAx, i:] 

= ~[Dx(t), x(t)] - [F(x(t),Dx(t)), m] 

[F(x, DJc), Jc(0] = ([x,^A x]) 2(n + p)+2 

> 0 

=> < o 

at 

GSI 



x = a k (t){ h 


F(x, Dx) 

= y (a k (t) (H k a k (t)) 2(n + ?)+1 (O k a k mk 

= Y a t «) 2( ' , * PK2 <o? < " P, * J ^W 2< " SK1 ^ 

a = P ; m = n 
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X(t) = a k (t) <$> k 


aicd ) + A,( 0 la k (t) + 2C,^X d) k a k (t ) 

+ y (a k (t)a k (t)d) k ) 2{n + ^ +l) co k a k (t) 

= 0 

Y a k (t) 2m \a k (t)\ a \a k (t)f a k (t) 2n+1 

~ m = n ; 


a = (3 
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Alternate Form 


u(s, t) + Xu""(s, t) - 2C^X ii"(s, t ) 


+ Y 




Lj 

J u(s, t) u"(s, t) ds 


1 2(« + P)+ 1 


0 


u"(s, t) 


= 0 , 


0 < s < L; 0 < t 


n : zero or positive integer 

0<p<i 

£ : Linear Damping Ratio 
Prime represents space derivative 
Dot represents time derivative 
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USE OF THE gUASIUNEARIZATION ALGORITHM 
FOR THE SIMULATION OF LSS SLEWING 

By 

Feiyue Li and P. M. Bainum 
Howard University 
Washington. DC 20059 

ABSTRACT 


The use of the Maximum Principle for the large angle dewing of LSS usually results In 
(he so-called two-point boundaiy-value problem. In which many requirements (e.g., 
minimum time, small amplitude, and limited control power, etc.) must be satisfied 
simultaneously. The successful solution of this problem depends largely on the use of 
an efficient numerical algorithm. There are many candidate algorithms available for 
this problem (e.g.. quasllinearizatlon, gradient, etc.). Here we discuss only the 
quaslllnearization method which has been used for several cases of large angle slewing 
of LSS. The basic idea of this algorithm is to make a series of successive 
approximations of the solution from a particular solvable case (linear or nonlinear) to 
a more general practical case. 

For the rigid spacecraft slewing problem with no constraints on the controls, the 
solution procedure can be found in the literature. This procedure needs to be modified if 
a minimum time for the slewing problem is desired with control limits given. Recently, 
an Indirect method for finding the minimum time was developed to meet all these 
requirements. 

For the general mixed (including both rigid and flexible parts) problem, an additional 
constraint of small vibrational amplitude on the flexible parts is Imposed. To solve 
this problem several steps in which the complexity increases gradually are needed, Le., 
from a linearized version to a final nonlinear problem, from a less constrained case for 
the control to a more constrained one. from a nonminimum-time level to a near- 
minimum-time slewing in which a trade-off needs to be made between minimum time 
and small flexural amplitude requirements. Some examples of these algorithms are 
presented for planar slewing maneuvers of the SCOLE configuration. 
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Use of the Quasilinearization Algorithm 
for the Simulation of LSS Slewing 
Feiyue Li 

Graduate Research Assistant , (202)636-7124 

and 

P. M. Bainum 

Professor of Aerospace Engineering, (202)636-6612 
Department of Mechanical Engineering 
Howard University, Washington, D.C. 20059 


Abstract 

The use of the Maximum Principle for the large angle 
slewing of LSS usually results in the so-called two-point 
boundary-value problem, in which many requirements (e.g., 
minimum time, small amplitude, and limited control power, etc) 
must be satisfied simultaneously. The successful solution of 
this problem depends largely on the use of an efficient 
numerical algorithm. There are many candidate algorithms 
available for this problem (e.g., quasilinearization, gradient, 
etc.). Here we discuss only the quasilinearization method which 
has been used for several cases of large angle slewing of LSS. 
The basic idea of this algorithm is to make a series of 
successive approximations of the solution from a particular 
solvable case (linear or nonlinear) to a more general practical 
case. 


For the rigid spacecraft slewing problem with no 
constraints on the controls, the solution procedure can be 
found in the literature. This procedure needs to be modified if 
a minimum time for the slewing problem is desired with control 
limits given. Recently, an indirect method for finding the 
minimum time is developed to meet all these requirements. 

For the general mixed (including both rigid and flexible 
parts) problem, an additional constraint of small vibrational 
amplitude on the flexible parts is imposed. To solve this 
problem several steps in which the complexity increases 
gradually are needed, i.e., from a linearized version to a 
final nonlinear problem, from a less constrained case for the 
control to a more constrained one, from a non-minimum-time 
level to a near-minimum-time slewing in which a trade-off needs 
to be made between minimum time and small flexural amplitude 
requirements. Some examples of these algorithms are presented 
for planar slewing maneuvers of the SCOLE configuration. 
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INTRODUCTION 


MAXIMUM PRINCIPLE IS APPLIED TO 

THE ATTITUDE MANEUVER AND VIBRATION CONTROL 

OF LARGE SPACE STRUCTURES 

(A) PERFORMANCE INDICES 

(B) BOUNDARY CONDITIONS 

(C) CONTROL REQUIREMENTS 

THIS LEADS TO THE TWO-POINT BOUND ARY- VALUE PROBLEM 
(TPBVP) 

ONE OF THE METHODS OF SOLVIN6 TPBVP IS THE 
QUASILINEARIZATION ALGORITHM 
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MAXIMUM PRINCIPLE 


STATE EQUATIONS 

X - f(x) ♦ B(x)u, x(0)-Xq, X(tf)-X( ( 1 ) 

PERFORMANCE INDICES 


Jj=(1/2) J (x T Qx ♦U T Ru)dt (2) 

0 

J 2 = f (Odt= t f |U|U u 1b , 1=1 ... n (3) 

0 

NECESSARY CONDITIONS 

H|=(1/2)(X T Qx ♦ U T RU) ♦ A T (f(x) ♦ Bu) (4) 

A = - (dHj/dx), A (0) unknown (5) 

(3Hj/du)=0, Ru=-B t A (6) 

H 2 = 1 ♦ A T (f(x) ♦ Bu) (7) 

A = - (dHj/ix), A (0) unknown (8) 

Uj= - Ujjj slgn(B T A ), i=1 ... n (9) 

TPBVP 

i = g(z) . z=[ x . A 1 T = [ z,. z 2 ] T (10) 

Z|(0). Zf(t f ) known; 
z 2 (0), z 2 (t f) unknown. 
z 2 (0) to be determined. 
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QUASILINEARIZATION ALGORITHM 


(A) LINEAR DIFFERENTIAL EQUATION: 

Nonhomogeneous: z - Az ♦ B. z=lzj. z 2 l T . (11) 

zj(0), Zj(t f ) known, z 2 (0) to be determined 
Homogeneous: z - Az (12) 

(a) n solns. of (12) ♦ 1 particular soln. of (11) 

(b) n ♦ 1 particular solns. of (11) 

(B) NONLINEAR CASE: 

Linearized equation of (10): 

z< k * * ) =(dg/dz) 1 ) ♦ h( z( k ) ) (13) 

where 

z'*' is the k** 1 approximate solution 

of the nonlinear equation (10), 
z (k+1) =z (k) +Az (k) 

z=tz,, z 2 1 T . 

Z|^ k+ ^)(0). Z|(* + *)(tf), known 
z 2 (k+1) (o) to be determined 
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(SCOLE) 
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ZL9 


TOWARDS 
THE CENTER 
OF THE EARTH 



Fig. 1b. Attitude of the SCOLE Showing Antenna Line of Sight 











ELANAR SLEWIN6 OF FLEXIBLE SCOLE 


LINEARIZED EQUATION OF MOTION: 

U S 
u l 
u 2 
[ u 3, 

where 

0 is the angle of rotation, 

7 nX | is the amplitude vector of the flexible modes, 

«> 

n is the number of mode used, 

1 is the moment of inertia about the axis of rotation 
■, M are the Inertia parameter vector, matrix. 

E is the stiffness matrix. 

<£(z) Is the mode shape function vector. 

<P j^(^) . Zj is the coordinate along z axis, 

L is the length of the beam, 
u s is the control torque on the Shuttle, 

U| are the control actuators on the beam and the 


I « T 


• • 1 

e 

♦ 

0 

fi T 

e 

s 

1 



£ 

■ II 


• « 
7 


0 

» 

t 

7 


0 



& 

1 


67S 


reflector. 



STATE EQUATIONS 


s - As ♦ Bu 



s l 


e 


s = 

s 2 

. Sl = 

'I 

. s2 = 

BOUNDARY CONDITIONS FOR s 




w 

© 


r i 

0 


s(0)= 

0 

, S(t f )= 

0 



0 


0 



0 


0 

2(n* 1 )x 1 


where n is the number of mode shapes used. 
PERFORMANCE INDEX 

l f 

J =( I /2) j (x T Ox *u T Ru)dt 

/ 

0 




TPBVP 


z = Cz. z=I s , A 1 T = lz t .z 2 ] T 

A is the costate vector, 

Z|(0), Zj(t f ) known; 
z 2 (0) to be determined. 
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FIGURE 1.1: DRAWING OF THE SCOLE CONFIGURATION 
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NUMERICAL RESULTS 


(A) SLEWING ABOUT X-AXIS (ROLL) 


OnUj ( -first skfxs ) IS used 

Of = 20 y 

J (s T QS + 2t r e u ) & s- f e,7, , e, ?, ] 



CASE / : V, /s used , 

Q, - L \] 


CA S£ 2 : 


U.J is KS*-^ 

<P*“ L 0 °- C ”° 0 ] 


CA SE 3 - : 



CA SE 3 : 


'll, Ar^. usee/ 
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see 
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Q*< 0.000 0.010 0.010 0.010 ) 


CONCLUDING REMARKS 


1) Solution has been obtained for nonlinear rigid spacecraft 
attitude maneuver (including the rigidized SCOLE). 

2) Use of the Maximum Principle can make the states 
satisfy the boundary conditions very well. 

3) Due the fact that the costates must be used in the method, 
the dimension of equations of the system is doubled, and 
higher computational ability is needed in this method. 

4) Further work on more complicated models (nonlinear 
differential equation) is needed. 

5) Need to consider different cost functions and perform 
parametric studies. 
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SESSION V - CONTROL SYNTHESIS AND SIMULATION 
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CONTROL LAW SYNTHESIS AND OPTIMIZATION SOFTWARE 
FOR LARGE ORDER AEROSERVOELASTIC SYSTEMS 

By 

V. Mukhopadhyay. A. Pototzky. and T. Noll 
NASA Langley Research Center 
Hampton, Virginia 


ABSTRACT 


Motivation: A flexible aircraft or space structure with active control is typically modeled 
by a large-order state space system of equations in order to accurately represent the rigid 
and flexible body modes, unsteady aerodynamic forces, actuator dynamics and gust 
spectra. The control law of this multi- input/multi-output (MIMO) system is expected to 
satisfy multiple design requirements on the dynamic loads, responses, actuator deflection 
and rate limitations, as well as maintain certain stability margins, yet should be simple 
enough to be implemented on an onboard digital microprocessor. This paper describes a 
software package for performing an analog or digital control law synthesis for such a 
system, using optimal control theory and constrained optimization techniques. 

Software Capabilities: The primary software capability is the optimization of the system 
by changing the control law design variables to improve stability and performance. A 
block diagram of the optimization scheme is shown in Fig. 1. 

1) The optimization module minimizes a linear quadratic Gaussian (LQG) type cost 
function, while trying to satisfy a set of constraints on the conflicting design 
requirements such as design loads, responses and stability margins. Analytical 
expressions for the gradients of the cost function and the constraints, with respect to the 
control law design variables, are used for computation. This facilitates rapid convergence 
of the numerical optimization process. The designer can choose the structure of the 
control law and the design variables. This enables optimization of a classical control law 
as well as an estimator-based full or reduced order control law. Selected design responses 
are incorporated as inequality constraints instead of lumping them into the cost function. 
This feature is used to modify a control law to meet individual root-mean-square (RMS) 
response limitations and design requirements. 

2) In order to improve the multiloop system stability robustness properties in the 
frequency domain, the minimum singular value of the return difference matrix at the 
plant input and output are as additional inequality constraints. 

3) Other supporting capabilities Include: (a) singular value analysis evaluation and 
plotting at the plant input and output: (b) linear quadratic optimal control law synthesis; 
(c) Kalman Filter design, LQG Loop transfer recoveiy; (d) pole-zero computation; (e) 
frequency response. Nyquist and Bode Plot; (1) root locus plot; (g) block dlagonalizatlon; (h) 
modal residuallzatlon and truncation; (1) transient response to deterministic and white 
noise input; Q) transfer of quadruple data to and from MATRIX-X and DIGIKON; (k) 
parameter search to stabilize an unstable control law, and (1) both interactive and batch 
mode execution using the Cyber NOS system. 
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2 


Applications: The software has been used In the past for the following applications: (l) 
flutter suppression control law for the AKW-I wind tunnel wing model; (2) gust load 
alleviation control law for the AKW-II drone: (3) flutter suppression control law synthesis 
for AKW-II drone and the DC- 10 Derivative wind tunnel wing model; (4) robust Digital gust 
load alleviation control law synthesis for AKW-II drone; and the (5) Active Flexible Wing 
(AFW) flutter suppression control law synthesis which Is presently being carried out. 
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FIG. 1 Optimization scheme block diagram 
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Abstract 


A flexible aircraft or space structure with active control is typically modeled by a large 
order state space system of equations in order to accurately represent the rigid body and 
flexible modes, unsteady aerodynamic forces, actuator dynamics and gust spectra. The 
control law of this multi-input multi-output (MIMO) system is expected to satisfy multiple 
design requirements on the dynamic loads, root mean square (RMS) responses, actuator 
deflection and rate limitations as well as maintain certain guaranteed stability margins, yet 
should be simple enough to be implementable on an onboard digital microprocessor. This 
paper describes an interactive software named DESIGN for analysis and synthesis of analog 
and digital control laws for such a system, using optimal control theory and constrained 
optimization techniques. 






Overview 


A multi-input multi -output aeroservoelastic system is typically represented by a large order 
state-space system of equations in order to accurately represent the rigid body and flexible modes, 
unsteady aerodynamic forces, actuator dynamics, gust spectra, antialiasing filters, computational 
delays etc. The active control law is expected to satisfy a set of conflicting design requirements on 
the performance and stability margins, yet should be simple enough to be implementable on an 
onboard digital microprocessors. This objective can be achieved using the synthesis software 
described in this paper. The methodology used are optimal control theory, order reduction 
techniques, unconstrained and constrained optimization with constraints on the design RMS 
responses and the minimum singular value of the return difference matrix at the plant input and 
output. Optimization can be performed for both continuous system and discrete systems. The 
methodology has been used to synthesize a) Analog and digital gust load alleviation control laws 
for a remotely controlled drone b) Analog and digital flutter suppression control laws for Active 
Flexible Wing (AFW) wind tunnel model . Other potential future applications include a) Rapid 
maneuver load control for AFW d) Vibration suppression for large space structure and control 
structure interaction study. 


OVERVIEW 


CONTROL LAW SYNTHESIS AND OPTIMIZATION 
SOFTWARE FOR FLEXIBLE STRUCTURE 


DESIGN OBJECTIVES 

LOW ORDER ROBUST CONTROL LAW FOR A 
HIGH ORDER AEROSERVOELASTIC SYSTEM 

METHODOLOGY 

OPTIMAL CONTROL THEORY 
CONTROL LAW ORDER REDUCTION 
NUMERICAL OPTIMIZATION 

COST FUNCTION LOG TYPE 
CONSTRAINTS RMS RESPONSES 

SINGULAR VALUES 

SYSTEMS 

CONTINUOUS DISCRETE 

APPLICATIONS 

GUST LOAD ALLEVIATION OF A DRONE 
FLUTTER SUPPRESSION OF AFW MODEL 
RAPID MANEUVER LOAD CONTROL 
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Optimization Block Diagram 


The optimization procedure minimizes a linear quadratic Gaussian (LQG) type cost function, 
while trying to satisfy a set of constraints on the conflicting design requirements such as dynamic 
loads, design RMS responses and singular value based stability margins at the plant input and 
output. The analytical expressions for the gradients of the cost function and the constraints, with 
respect to the control law design variables are used for computation. This facilitates rapid 
convergence of the optimization process. The designer can choose the structure of the control law 
and the design variables. This enables optimization of classical control law as well as an estimator 
based full or reduced order control law. Selected design responses are incorporated as inequality 
constraints instead of lumping them into the cost function. This feature is used to modify a control 
law to meet individual RMS response limitations and design requirements. 


DISCRETE SYSTEM 
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| CONTINUOUS SYSTEM CONTROL LAW DESIGN 
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CONTINUOUS PLANT + CONTROL LAW 


REDESIGN 
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METHOD OF 
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CONSTRAINTS 


DIRECTION 
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Software Organization 


The interactive software DESIGN is organized to interact with several well used 
softwares such as 1) IS AC (Interaction of Structure, Aerodynamics and Control) for receiving 
state-space quadruple data, 2) DIGKON for discretization, interconnection, model 
generation, digital design, verification and graphics and 3) MATRIX-X for matrix 
manipulation, interconnection, quadruple data transfer, graphics and design 
verification. DESIGN can also be run in batch mode on the CYBER/NOS system for large 
order problems involving systems with more than 120 states with large number of design 
variables and constraints. This batch version was previously known as PADLOCS 
(Program for Analysis and Design of Linear Optimal Control Systems). 


isac \ 


l 


DIGIKON 


DESIGN 


MATRIX 


► 

INTERACTIVE 


X 



PADLOCS 




BATCH 
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Basic Command Summary 


The quadruple data is generated and stored in a sequencial binary file called QDATA. The 
design starts with the file command 

GET, QDATA. 

GET, DESIGN. 

DESIGN. 

The random access files DBASE, and sequencial file PLDATA are used to transport quadruple 
data to and from DIGIKON and MATRIX-X, while random access file TAPE7 is ued to transfer 
data from IS AC.using the UTILITY commands. The system parameter and quadruple data are 
read by the SYSTEM INPUT commands as shown in the figure above. The primary capability of 
this software is the optimization of the system by changing die control law design variables to 
improve the stability robustness and performance requirements. The supporting capabilities include 
a) Linear quadratic optimal control law synthesis; b) Kalman filter design, linear quadratic 
Gaussian design (LQG) and loop transfer recovery (LTR); c) Singular value analysis , evaluation 
and plotting at the plant input and output; d) Pole-zero computation; e) Open and closed loop 
frequency response, Nyquist and Bode plot, and loop breaking test; f) Root locus plot g) Block 
diagonal transformation; h) Modal residualization and truncation; i) Transient response to 
deterministic and white noise input; etc. 
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Basic Design Commands 


The basic design commands are shown in the figure above. For systems with known stable 
control laws the optimization procedure can be executed directly using the command EXC OPTM 
for continuous systems and EXC OPTD for discrete systems. For MIMO systems with no known 
initial stabilizing control laws, first an Linear Quadratic Gaussian (LQG/LTR) design is performed 
to obtain a full order robust control law using a set of LQG design commands. The order of the 
control law is then reduced by truncation, residualization or balanced realization method using 
DESIGN, DIGIKON or MATRIX-X. The singular value analysis and block diagonal 
transformation procedure is very helpful in the reduction process. Since this reduced order control 
law is not optimal and may not satisfy the design requirements, constrained optimization procedure 
is used to update the reduced order control law. Constraints can be imposed on the design RMS 
responses and minimum singular values at the plant input and output. 


BASIC DESIGN COMMANDS 
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Gust Load Alleviation of A Flexible Drone 


The synthesis procedure was applied to the gust load alleviation problem of a flexible drone. 
The basic control scheme is shown in the figure. In longitudinal motion, the symmetric elevator 
and outboard aileron deflections are used as the two control inputs. The accelerometer sensors at 
the outboard aileron and on the fuselage near the center of gravity are used as two measurement 
outputs. The output signals are filtered through first order antialiasing filters 50/(s+50) before 
digitization at 100 Hz. The two input two output system was modeled by a 32nd order system 
flying symmetrically through a Dryden gust. 


GUST LOAD ALLEVIATION OF 
FLEXIBLE DRONE 



703 




Gust Load Alleviation Design Requirements 

The objective is to obtain a low order robust digital GLA control law which would reduce the 
open loop root-mean- square values of the wing root bending moment and shear by 50% without 
increasing the wing outboard bending moment and torsion The control law should maintain certain 
guaranteed stability margins based on minimum singular value of 0.6 at both the plant input and 
output. The control surface deflections and rates should be within the allowable limits. First a full 
order LQG control law is synthesized to satisfy the design requirements. This 32nd order control 
law is then reduced to a second order control law and then discretized. This control law does not 
satisfy the design requirements. After unconstrained optimization most of the requirements are 
satisfied except the wing outboard bending moment and the singular values. Using constraints on 
the RMS wing loads and on the minimum singular values of the return difference matrix at the 
plant input and output, the control law parameters are reoptimized (Ref. 1,2). 


GUST LOAD ALLEVIATION 
DESIGN REQUIREMENTS 



32nd order 
airplane eqns. 



control law 


Physical quantities 

Design objectives 

How we do it 

Root bending moment 

50% reduction 

1. LQG design 

2. Control law 
order reduction 

3. Discretize 

4. Optimization 

5. Apply constraints 

a) on rms loads 

b) on singular val. 

Root shear 

50% reduction 

Outboard bending mom 

No increase 

Outboard torsion 

No increase 

Elevator deflection 

Within max limit 

Elevator rate 

Within max limit 

Aileron deflection 

Within max limit 

Aileron rate 

Within max limit 


7 0 + 






















Symmetric Flutter Suppression System 

The software has been used in the past for the following applications: a) Robust flutter 
suppression control law synthesis for ARW-I wind tunnel wing model; b) Flutter suppression 
control law synthesis for ARW-II drone and DC- 10 derivative wind tunnel wing model and c) s 
plane summation of forces load model (Ref. 3). A brief survey of the research activities is 
presented in RefADigital robust control law synthesis for the Active Flexible Wing (AFW) wind 
tunnel model is presently being carried out in collaboration with Rockwell International. The basic 
block diagram for a two input two output symmetric flutter suppression system is shown in the 
figure for a sting mounted model using leading edge outboard (LEO) and trailing edge outboard 
(TEO) symmetric actuators and colocated accelerometer sensors.The sampling rate is 200Hz.The 
design takes into account the effects of actuator dynamics, 4th order 100Hz Butterworth filters and 
one cycle computational delay at each channel. Full order and reduced order analog and discrete 
robust control laws were synthesized based on an approximate 38th order system at 300 psf design 
dynamic pressure. The discrete 8th order control law was able to stabilize the system over the 
range 300 to 150 psf.The more detailed 80th order model was also stable at 300 and 200 psf. 
Starting with these preliminary control laws detailed analysis will be carried out using the discrete 
system optimization procedure EXC OPTD. 


SYMMETRIC 


FLUTTER SUPPRESSION 
SYSTEM 


IVA adtutof anil alasing 
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CONCLUDING REMARKS 


• Software Improvement 

1. Direct time response constraints 

2. New reduction techniques 

3. H-infinity design 

4. Additional derivative/sensitivity capabilities 

5. User help and online documentation 

• Portability improvement 

1 . Microvax 

2. NOS VE 


• Future Applications 

1 . AFW rapid roll maneuver 

2. Large space structure-control interaction 


AFW RAPID ROLL CONTROL 


Active Flexible Wing Rapid Roil Mechanization 
Wind Tunnel Model Digital Control 



LARGE SPACE STRUCTURE 
ACTIVE CONTROL 
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FLEXIBLE AIRCRAFT DYNAMIC MODELING 
FOR DYNAMIC ANALYSIS AND CONTROL SYNTHESIS 

By 

David K. Schmidt 
Purdue University 
West Lafayette. Indiana 


ABSTRACT 


The linearization and simplification of a nonlinear, literal model for flexible aircraft 
is highlighted. Areas of model fidelity that are critical if the model is to be used for 
control system synthesis are developed and several simplification techniques that can 
deliver the necessary model fidelity are discussed. These techniques include both 
numerical and analytical approaches. An analytical approach, based on first-order 
sensitivity theory is shown to lead not only to excellent numerical results, but also to 
closed-form analytical expressions for key system dynamic properties such as the 
pole/zero factors of the vehicle transfer-function matrix. The analytical results are 
expressed in terms of vehicle mass properties, vibrational characteristics, and rigid- 
body and aeroelastic stability derivatives, thus leading to the underlying causes for 
critical dynamic characteristics. 
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TOPICAL OUTLINE 


What Constitutes a Valid Model? 

How Validity Will Be Measured 
What’s Important in a Feedback System 

Some Approaches to Obtain (Simple) Valid Models 
Numerical 
Literal 

Physical Causes of Critical Dynamic Characteristics 



WHAT VEHICLE DYNAMIC CHARACTERISTICS 
ARE CRITICAL IN A FEEDBACK SYSTEM? 

They Are Important If They Can: 

— Induce Critical Pole/Zero Interactions 
— Significantly Affect Frequency Response Near Crossover 
— Significantly Affect Time Response 
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Table 6. Frequency Weighted Internally Balanced Reduction 
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Step 2: Find T and I where XY = TE 2 r‘,T = (T t T M ], r T = [U„ U M ] 
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FOR QUASI-STEADY AERODYNAMICS 
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Tabfe 4. Elastic Alttrtft Linear Longitudinal Equations In Polynomial Matrix Form 
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Table 10. Transfer Functions For the True Model 
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Table 12. Transfer Functions ForFWIB Reduction Model 
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Table 14. Lower-Order Transfer Function Forms 
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Table 8. Polynomial Coefficient and Factor Relationships 
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Table 9. Polynomial Coefficient and Factor Difference Relationships 
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Table 14. Transfer Functions from Literal Approximations 
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CLOSED-FORM EXPRESSIONS 
FOR THESE TERMS 


Specifically, it has been shown that for a slender (low aspect ratio) vehicle 
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FURTHERMORE, IT CAN BE SHOWN 


For Example, That 
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f c wdy w + s c c 


-b/2 C W ' dx 


where 


c i = wing section lift-curve slope 

a w 

(-— ) r .\ = mode i slope at location (•) 
dx w 
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SUMMARY 


• Key Issues In Feedback Systems Reviewed 

• Physics Of Model Uncertainty Addressed 

• Mode Frequencies Near Crossover 

• Mode Shapes And Dipoles 

• With Aeroelastic Coupling, Both Affected 

• Closed-Form Expressions Developed To Show Sources Of Interactions 

• Sensitivity To Uncertainty In These Parameter May Be Further Explored 


( Purely 
r Elastic 
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INTRODUCTION 


• Flexibility is important for high speed, 
high precision operation of lightweight 
manipulators. 

• Accurate dynamic modeling of flexible 
robot arms is needed. Previous work has mostly 
been based on linear elasticity with prescribed 
rigid body motions (i.e., no effect of flexible 
motion on rigid body motion). 

• Little or no experimental validation of 
dynamic models for flexible arms is available. 
Experimental results are also limited for 
flexible arm control. 

• We include the effects of prismatic as 
well as revolute joints. 

• We investigate the effect of full coupling 
between the rigid and flexible motions, and of 
axial shortening. 

• We also consider the control of flexible 
arms using only additional sensors. 





BACKGROUND 


Research since 1970's 

(e.g., [Book, Maizzo-Neto, Whitney 75]) 

Modeling of flexible mechanisms and structures 
(e.g., Elasto-Kineto Dynamics, Floating 
Frames, 70's) 

Approaches to control 

Trajectory planning [Meckl, Seering 83,85] 

Open loop (none) 

Closed loop with micromanipulator 
[Cannon et al, Book et al] 

Closed loop with additional sensors only (none) 

Experimental work 

[Zalucky and Hardt 84] 

[Cannon et al 83, 84] 

Theoretical control studies 

[Book et al, Cannon et al, etc, early 1980's] 
Various control strategies proposed typically 
assuming all states available, no spillover, simple 
models, no implementation considerations. 
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MODELING AND SIMULATION OF FLEXIBLE 
ROBOTS WITH PRISMATIC JOINTS 


• Robots with both rigid and flexible links 
attached with revolute and/or prismatic joints 
can be modeled and analyzed. 

• The equations of motion are derived using 
Lagrange's equations. The prescribed motion, 
and prescribed torque/force cases can both he 
handled. 

• Flexible elements are represented as 
Euler-Bernoulli beams, and the axial shortening 
effect is also included. 

• Finite element analysis is used for the 
discretization of the resulting hybrid equations 
of motion. 

• Constraints are handled using Lagrange 
multipliers. 

• The resulting algebraic-differential 
equations are solved numerically using 
constraint stabilization methods. 



■ Nominal configuration 
B Actual configuration 



Schematic of a two-link robot. 
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X P* 




Axial shortening of a beam under plane transverse deflection. 



Schematic of revolute joint i. 
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-H K- 



link i- 1 
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Time (sec.) 


Tip displacement in "fast pull" case with 
Ci = 0.025 m, C 2 = -0.7 m and T = 0.7 sec. 
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LABORATORY ROBOT 

Small table top spherical coordinate robot with 3 DOF 

Designed and built at UM 

Interfaced to an IBM PC/XT 

Convienent test bed experimental research work 

r and 0 axes are dc motor driven through leadscrews 

+ axis is dc motor driven driven through a gear train 

all axes have tachometers and optical incremental 
encoders with counter circuits 

last link is intentionally designed to be flexible 

accelerometers (in two orthogonal directions) 
measure end of arm accelerations which are 
integrated to get velocities and positions 
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Lead Screw 



Schematic of the experimental setup. 
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Elastic lip displacement (m.) 



The elastic tip displacement obtained from the equations of motion 
with prescribed motions. 


.03 
.025 

~ .02 
•§ .015 
| .01 
8 .005 

Q. 0 

1L -.005 

o -.01 

M -.015 

UJ 

-.02 
-.025 
-.03 

0 .5 1 1.5 2 2.5 3 3.5 4 4.5 5 

Time (sec.) 

Elastic tip displacement obtained from the equations of motion 
with prescribed torque/force. 
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Translational velocity (m./s.) Angular velocity (dcg./s.) 



Angular velocity of the first joint of the rigid manipulator 
with controller. 



Translational velocity of the third joint of the rigid manipulator 
with controller. 
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l ' I » I ' T » I 1 1 * I * T l T 

0 .5 1 1.5 2 2.5 3 3.5 4 

Time (sec.) 


Angular velocity of the first joint of the flexible manipulator 
with controller. 



Time (sec.) 

Translational velocity of the third joint of the flexible manipulator 
with controller. 
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Elastic tip di 



:cments in the case with T = 2 sec. 



icement in the case with T = 1 sec. 
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SUMMARY AND CONCLUSIONS 


• A general modeling procedure for robot 
arms consis ting of rigid and flexible links 
connected by revolute and/or prismatic joints 
has been developed and experimentally 
validated. 

• The significance of full coupling (effect 
of flexible motion on rigid body motion) has 
been demonstrated. 

• The axial shortening effect is shown to be 
significant for high speed operation of 
lightweight manipulators. 



CONTROL OF A LEADSCREW DRIVEN FLEXIBLE 

ROBOT ARM 


• The laboratory robot is used to compare 
the performance of a rigid body motion 
controller with that of a rigid and flexible 
motion controller. 

• The rigid body motion controller uses 
only the joint motion measurements and joint 
actuators. The rigid and flexible motion 
controller also uses the end of arm motion 
measurements, but no additional actuators. 

• The leadscrew transmission 
characteristics as veil as observation and 
control spillover are considered. 

• The numerical and experimental results 
show good agreement, and indicate that 
significant reductions in arm vibration are 
possible through use of the rigid and flexible 
motion controller. 
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o 

K k 


Arm geometry and coordinates. 
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PHYSICAL CONSTRAINTS 


The physical constraints that are considered in this 
work are the ones imposed by the leadscrews only. 

. Condition for self locking assumption to be valid 
is: 

H > tan( Vl ) 

where 

\l is the thread coefficient of friction. 

V, is the thread helix angle. 

. Effect of the self locking condition. 

. Effect of coulomb friction. 
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CONTROLLER DESIGN 


Equations of motion: 

M(*) r £ (* ,r) = f(T) 

*= L r ,0> 4 > f/i , %/z > ‘fra > ? 2 * 7 

I* r = L~ r t /* , T 5 ] 

Linearized equations: 

^ ^ + & U 

y T = £«JV r «Ttf T ] 5 <d=£l 

<J* T r £ir, <fe, S4> , $<£„ j Sj tl ] 

Integral plus state feedback controller: 

ijtl- ; fcy 2 -R z )dt ; y lz 
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Block diagram of the integral plus state feedback controller. 
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fluUt imaart 



Homj kttd 

Figure 3. Flexible motion coordinate, qj 1 (t), In response to 

the rigid and flexible motion controller In the 
reduced order model case. 


a^») Insert 



Figure 4. 


Day imnf 

Flexible motion coordinate, qgiO)* In response to 

the rigid and flexible motion controller In the 
reduced order model case. 
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fluULtofiiad 



2kDfcJ iAMUt 

Figures. Flexible motion coordinate, q^t), In response to 
the rigid and flexible motion controller In the 
control spillover case. 



Figure 6. 
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Flexible motion coordinate, q 12 (t), In response to 
the rigid and flexible motion controller In the 
control spillover case. 
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Dmj l i cait 

Figure 7. Flexible motion coordinate, q-f 2 <t), In response to 

the rigid and flexible motion controller In the 
control and observation spillover case. 
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Figure 8. 


Flexible motion coordinate, q-^O)* |n response to 

the rigid and flexible motion controller In the 
control and observation spillover with structural 
damping Included. 
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Standard Set of Physical System 
Parameters 


VALUE 


Mass of the first beam (m,) 

0.454 Kg 

Mass of the second beam (m s ) 

0.816 Kg 

Mass of the Payload (m, ) 

0.07 Kg 

Cross sectional area of the second 


beam (A 2 ) 

0.000151 m 2 

Length of the first beam ( L ,) 

0.233 m 

Length of the second beam {L t ) 

2 m 

Gravitational acceleration (g ) 

0.81 m/te e : 

Aluminum density (p) 

2707 Kg/m * 

Flexural rigidity (El) 

770.87 Pa 

Reference position for r 

1.85 m 

Reference position for 0 \ 

0 rad 

Reference position for ^ 

0 rad 

Desired reference position for r 

2m 

Desired reference position for 6 

0.5 rad 

Desired reference position for d 

0.5 rad 

Servo nataral frequency for r (u>„ ) 

. 4 rad/sec 

Servo naMral frequency for # (w, ,) 

■ 4 rad/»ec 

Servo nataral frequency for d (w.*) 

8 rad/sec 

Flexible motion gain, K £ 

•0.000178 

Flexible motion gain, K £» 

•0.084 

Flexible motion gain, K 

i 1.568 


TABLE 1 



settling time 
(seconds) 

maximum deflection 
(peak to peak) 

rigid body controller 

11.0 

7.5mm 

rigid and flexible 
motion controller 

3.0 

2.7mm 


Table 2. 
110 
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TW I jMOXSf 

Figure 10. e response obtained from the rigid body controller In 
the experimental work. 
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{mm) 



24.00 


TVMLt iMODOl 

Figure 11. Total vertical deflection In response to the rigid 
body controller In the experimental work. 
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T>na_ t iMOBBtf 

Figure 12. • response obtained from the rigid and flexible 

motion controller In the experimental work. 


CaMSgn* 

t«*»l 



QttU iMBDOt 

Figure 13. Control signal for the second Joint obtained from the 
rigid and flexible motion controller In the 
experimental work. 
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Bnaj ImobA 

Figure 14 . Total vertical deflection In response to the rigid and 
flexible motion controller In the experimental work. 


77 3 



RIGID BODY CONTROLLER VERSUS 
RIGID AND FLEXIBLE MOTION CONTROLLER 


Simulation results: 

* Control spillover effect can be observed, but 
does not cause significant deterioration. 

* Control and observation spillover can 
destabilize the residual mode. However, a 
small amount of damping (0.0145) eliminates 
the problem. 

* Settling time is reduced from 3.5 to 1 .07 
seconds, and maximum vibration amplitude is 
reduced by 50%. 


Experimental results: 

* With low pass filtering and light structural 
damping, no detrimental spillover effects were 
observed. 

* Settling time is reduced from 1 1 to 3 
seconds, and maximum vibration amplitude is 
reduced by 75%. 
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SUMMARY AND CONCLUSIONS 


A dynamic model of a spherical coordinate robot arm, 
whose last link is flexible, is developed. The 
constraints imposed by the leadscrew transmission 
mechanisms are also considered. 

The interrelationships between the robot arm 
structural flexibility and the controller design 
are investigated using a rigid body controller. 

The rigid and flexible motion controller, which employs 
additional sensors only, has led to an approximate 50% 
reduction in the magnitude of the flexible motion even 
in the presence of the observation and control 
spillover. 

The experimental results of the rigid and flexible 
motion controller show good agreement with those 
of the digital simulation. 


SUMMARY AND CONCLUSIONS 

• A general modeling method for robot arms 
with flexible and rigid links connected by 
prosmatic and revolute joints has been 
presented and experimentally validated. 

• A flexible arm controller which uses end 
of arm motion meaurements, but only joint 
actuators has been numerically and 
experimentally studied and found to give good 
rigid body control with significant reduction in 
end of arm vibrations. 
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CONTROLLING FLEXIBLE STRUCTURES - A SURVEY OF METHODS 


By 

A. A. Anlsslpour, R A Benson, and E. E. Coleman 
The Boeing Company 
Seattle, Washington 


ABSTRACT 


In response to the demand for higher fuel efficiency, new Jet transport designs are using 
composites to reduce weight. Since these composites have less Inherent damping than the 
materials they are replacing, flexible body dynamics are becoming more significant. Low 
frequency structural modes that in the past had sufficient open loop damping to avoid 
stability and ride performance problems now require compensation. Control law design 
can no longer be concerned solely with rigid-body dynamics. 

Control systems design techniques which have proven effective for rigid dynamics may no 
longer be applicable when Including flexible modes. This paper studies the different 
control strategies and their application to systems Including flexible dynamics. The 
control law design goal with respect to flexible modes Is twofold. First, the modal damping 
must be Increased to some minimum acceptable level. Second, the rigid body modes must 
be controlled while minimizing excitation of the flexible modes. A necessary, and often 
nontrivial first step Is identification of the flexible modes (both mode frequency and 
shape). The following addresses these Issues with respect to a particular control law design 
technique. 

1) Off-line Modeling/ OfT- line Controller Design 

The off-line modeling controller design technique has been used historically as the 
standard for airplane control system design. This technique consists of building up a 
model on mathematical prediction, wind tunnel analysis, and flight test data. Controllers 
are then designed off-line based on the model. Controller implementation may Include 
scheduling based on flight condition, but with no on-line modifications made to the 
controller. 

2) Off-line Modellng/On-llne Controller Design 

The same type of model used for the previous technique Is developed using the off-line 
modeling/on-line controller design with the primary difference being that the controller 
structure Is defined off-line, while the controller gains are adapted on-line to minimize 
certain performance criteria. With this technique no attempt is made to model the open 
loop system on-line, and all controller gain changes are made solely in response to 
performance criteria. 
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3) On-line Modcllng/Off-lIne Design 

On-line modeling/off-line design technique defines the controller In terms of a general 
open loop system model for a specified structure. With this technique each controller gain 
Is defined as a function of model parameter values. Implementation Is executed by 
building an on-line adaptive open loop estimator. The controller gains are then set based 
on the model estimate values per the definitions developed off-line. 

4J Off-line Modellng/Off-llne Conlroller/On-llne Adjustment 

The control law Is designed using the off-line modeling/ off-line controller to stabilize the 
nominal plant. With this technique an adaptive loop is placed around the system, creating 
an adaptive system that tunes the control input to improve the ofT-nominal performance. 
This adaptive loop is designed to never destabilize the nominal plant and when the plant 
is at nominal, the adaptive control signal Is zero. 
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CONTROLLING FLEXIBLE STRUCTURES - A SURVEY OF METHODS 
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EDWARD E. COLEMAN 

The Boeing Company 
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P. O. Box 3707 

Seattle, Washington 98124-2207 


ABS T RACT 

Most of the presently available control system design techniques applicable to 
flexible structure problems were developed to design controllers for rigid body 
systems. Although many of these design methods can be applied to flexible 
dynamics problems, recently developed techniques may be more suitable for flexible 
structure controller design. The purpose of this presentation is to examine briefly 
the peculiarities of the dynamics of flexible structures and to stimulate discussion 
about top level controller design approaches when designing controllers for flexible 
structures. 

This presentation contains a suggestion of a set of categories of design methods for 
designing controllers for flexible structures as well as a discussion of the advantages 
and disadvantages of each category. No attempt has been made herein to select one 
category of design techniques as the best for flexible structure controller design. 
Instead, it is hoped that the structure suggested by these categories will facilitate 
further discussion on the merrits of particular methods that will eventually point to 
those design techniques suitable for further development. 


CHARACTERISTICS OF FLEXIBLE STRUCTURE DYNAMICS 

Flexible structure dynamics tend to differ from rigid body dynamics in several 
important ways. First, flexible dynamics are higher order than rigid body dynamics. 
By definition, rigid body dynamics involve six degrees of freedom. Since each 
degree of freedom results in two states, the full set of dynamics for a rigid body 
system will involve only twelve states. (Additionally, servo, actuator, sensor, and 
controller compensation states must be added.) By comparison, a flexible structure 
model may have 100 or more states. This increase in the number of states derives 
an increase in the complexity of the control problem. Hence, those design 
techniques which work well for tenth order systems may have difficulties handling 
systems with ten times that many states. 

A second important difference between rigid body and flexible dynamics is flexible 
dynamics tend to be more difficult to predict than rigid body dynamics. It is the 

78 \ 


structure of a system which derives its flexible dynamics. Parameters such as mass 
distribution, material stiffness and damping, and unsteady aerodynamics become 
influential. Often mathematical models developed to predict the flexible dynamics 
differ with the physical system representation. As a result, the controller design 
based on these prediction models must be made robust to withstand the 
discrepancies between the model and physical system. 

A third difference between these dynamics is that rigid body dynamics can often be 
treated as decoupled, whereas flexible body dynamics are most often highly coupled. 
As a result, control problems that can often be treated as a single input / single 
output (SISO) or as a series of SISO problems when dealing with rigid body systems, 
become multiple input / multiple output (MIMO) problems when dealing with 
flexible systems. SISO methods appropriate for decoupled rigid body system 
controller design may be unsuitable for flexible system controller design. 

A fourth difference between these dynamics is the goal of the systems designed to 
control them. Rigid body control usually involves commanding the rigid degrees of 
freedom to follow desired trajectories. For an airplane these might be altitude, 
heading, and airspeed. By contrast, the goal of most flexible structure controllers is 
either to perform the desired rigid body control without exciting flexible modes, or 
to provide active damping for structural modes that are excited. In almost all 
applications the objective is to keep flexible structure dynamic responses at the 
lowest possible levels. The difference between the goals for rigid body control and 
flexible structure control may require different controller design approachs. 


CONTROLLER DESIGN METHOD CATEGORIES 

Four categories of controllers are detailed in the following paragraphs with a 
discussion on the advantages and disadvantages of each category with respect to the 
design of controllers for flexible structure systems. It is not the author's intent to 
favor any category of controller design techniques over another. The divisions 
herein are made simply to facilitate comparison of different top level strategies for 
the design of controllers for flexible structures. The categories are deliniated by the 
types of models each use for controller synthesis, and whether the controller is 
designed off-line, on-line, or both. 

1) Off-line Modeling / Off-line Controller Design 

Off-line modeling / off-line controller design techniques have been historically used 
as the standard for control system design. This technique consists of building up a 
model using mathematical prediction, wind tunnel analysis, and flight test data. 
Controllers are then designed off-line based on the model. Although separate 
controllers may be designed for different flight conditions (i.e., requiring gain 
scheduling based on flight condition), only these on-line modifications defined 
previous to flight are made to the controller. 
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This category of controller design techniques has the advantage of all controller 
synthesis work being completed ahead of time off-line. As a result, the on-line 
computational load is kept to a minimum. In addition, since the controller is well 
defined for each flight condition, rigorous analysis is possible for predicting 
performance and robustness characteristics. 

A disadvantage of this controller design approach is that discrepancies between the 
model and the physical system itself must be handled solely by controller 
robustness. The controller is unable to tune itself to account for modeling errors or 
changes in the dynamics as a result of different flight conditions or weight 
distributions. Hence, this design approach requires development of an accurate 
system model. Whenever possible, the off-line model is updated to concur with the 
obtained test data using the physical system to be controlled. In those cases where 
test data is not available, analysis must be done to show that the controller will 
function in an acceptable manner for the set of anticipated dicrepandes between the 
model and the physical system. 

2) Off-line Modeling / On-line Controller Design 

The same type of model used for the previous category of controller is used for off- 
line modeling / on-line controller design techniques. However, the primary 
difference with this category is that while the controller structure is defined off-line, 
the controller gains are adapted on-line to minimize certain performance criteria. 
With this technique no attempt is made to model the open-loop system on-line, and 
all controller gain changes are made solely in response to performance criteria. 

This approach has the advantage of the controller being able to tune itself to account 
for parameter variations between the model and the physical system. The off-line 
system model does not have to be as precise for this type of controller as for those 
described in category 1. The off-line model is used to note the structure of the 
system and the general trend of the dynamics. Furthermore, a related advantage is 
that since the controller is able to tune itself, gain scheduling does not have to be as 
detailed as for a controller that is designed completely off-line. 

A disadvantage of this approach is that the controller must be tuned on-line, thus 
requiring more computation power. Another disadvantage is that while the 
controller is able to tune itself to account for parameter variations, its structure is 
fixed. If the structure of the flexible system changes or there are wide swings in its 
general dynamics, the controller may not be able to tune itself sufficiently to provide 
the necessary control. A third disadvantage is that while the controller is 
continuously tuning itself, it is impossible to predict the gains for any given flight 
condition. As a result, it is impossible to obtain the level of performance and 
robustness analysis possible with each flight condition assigned a fixed set of 
controller gains. 
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3) On-line Modeling / Off-line Design 

On-line modeling / off-line design techniques define the controller in terms of a 
general open-loop system model of a specified structure. With these techniques 
each controller gain is defined as a function of model parameter values. 
Implementation consists of building an on-line adaptive estimator to estimate the 
model parameters. The controller gains are then set based on the model estimate 
values per the definitions developed off-line. 

An advantage of this approach is that the controller gain definitions can be chosen 
to give the desired performance and robustness properties regardless of the model 
parameters. (It is assumed that sufficient controllability and observability exist for 
all variations of the model parameters.) Gain scheduling is not an issue since the 
controller gain definitions automatically give appropriate gains for any given 
operating condition. 

A disadvantage of this approach is that the controller gains must be defined 
symbolically rather than numerically. Controller synthesis requires solution of 
symbolic rather than numeric equations. Fortunately this task is done off-line and 
will not require real-time computing resources, but it is still a formidable task 
nonetheless. Another disadvantage is that on line estimation of the model 
parameters is required along with computation of the controller gains. Evaluation 
of the equations defining the controller gains may be quite expensive to compute. 
An additional disadvantage is that the controller gain definitions may include 
singular points within the region of possible model parameter sets. A method to 
avoid singularity is needed. 


4) Off-line Modeling / Off-line Controller / On-line Adjustment 

The final controller design approach category involves off-line modeling, off-line 
nominal controller design, and on-line controller adjustment. The nominal 
controller is designed to stabilize the nominal system. With this technique an 
adaptive loop is placed around the system to tune the control input for improving 
the off-nominal performance. This adaptive loop is designed to consistently 
stabilize the nominal system. Furthermore, when the system is at nominal, the 
adaptive control signal is zero. 

An advantage of this system is that while the controller is able to tune itself to 
account for modeling errors, the tuning is restricted so that the closed-loop system 
remains stable. This method can be thought of as a compromise between a fixed 
gain controller and fully adaptive controller. The nominal controller gains are fixed 
while the adaptive algorithm is free to vary the gains within a range about nominal. 

A disadvantage of this system is that if the physical system varies greatly from the 
nominal given by the off-line model, there is no longer a guarantee of stability. 
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Hence on the one hand, the adaptive tuning is restricted to keep from destabilizing 
the nominal system. On the other hand, this restriction may lead to a situation 
where the actual system is driven unstable and the limitations on the adaptive 
tuning are such that the controller is unable to tune itself sufficiently to provide 
closed-loop stability. 


SUM M ARY 

Rather than point to a single design approach as the best for designing controllers 
for flexible structures, the goal of this presentation is to simulate thought and 
discussion. Most likely a single approach is not well suited for all problems. The 
key is to realize that there are fundemental differences between rigid body and 
flexible structure dynamics and that these differences may require different 
approaches to controller design. 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


Goal of Presentation: 


Stimulate discussion of merits and 
maturity of different top level approaches 
to controller design for flexible structures 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 

Outline of Presentation: 

o Characteristics of Flexible Structure Dynamics 
o Four Categories of Controller Design Approaches 


o Open Discussion on Merits and Maturity of Design Approaches 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


Characteristics of Flexible Structure Dynamics 

o High order (100 or more states) 

o Difficult to predict 


o Coupled dynamics 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


Four Categories of Controller Design Approaches: 

1) Off-line Modeling / Off-line Controller Design 

2) Off-line Modeling / On-line Controller Design 

3) On-line Modeling / Off-line Design 

4) Off-line Modeling / Off-line Controller Design / 

On-line Adjustment 



BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


1) Off-line Modeling / Off-line Controller Design 

o Model developed off-line using mathematical predictions, 
wind tunnel data, and/or flight test data 

o Controller designed based on model 

o Controller structure and gains may vary from condition 
to condition but controller is fixed for a given operating 
point 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


Advantages: 

o Low on-line computation requirement 
o Well defined controller allows for rigorous analysis 

Disadvantages: 

o Controller cannot react to modeling errors 
o Model must match physical system closely 
o Controller must have significant robustness 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


2) Off-line Modeling / On-line Controller Design 

o Model developed off-line using mathematical predictions, 
wind tunnel data, and/or flight test data. 

o Controller structure chosen based on off-line model 

o Criteria outputs specified to drive controller adaptation 

o Controller gains adapted on-line in effort to minimize 
criteria outputs 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


Advantages: 

o Controller can tune itself to overcome modeling errors 

o Model must represent structure and general trend of 
the dynamics, but need not be exact 

Disadvantages: 

o On-line gain tuning requires greater computational power 

o Gains may diverge resulting in an instability 

o Because gains are not fixed for a given flight condition, 
rigorous closed loop analysis is not possible 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


3) On-line Modeling / Off-line Controller Design 

o System model adapted on-line to fit a specified structure 

o Controller gains defined as functions of model parameters 

o Controller gain equations developed off-line as functions 
of generalized model parameters 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


Advantages: 

o Controller gain equations give explicitly the desired 
performance and robustness characteristics 

o Controller gain equations developed off-line and must 
be computed only once 

Disadvantages: 

o Controller gain equations must be solved symbolically 
rather than numerically 

o On-line model estimation is computationally expensive 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


4) Off-line Modeling / Off-line Controller Design 
/ On-line Controller Adjustment 


o Model developed off-line using mathematical predictions, 
wind tunnel data, and/or flight test data 

o Controller structure defined based on off-line model 

o Nominal controller gains chosen based on off-line model 

o Controller gains allowed to adapt on line in such a way 
that the nominal plant (as defined by the model) is 
never driven unstable 
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BOEING COMMERCIAL AIRPLANES 

GUIDANCE AND CONTROL RESEARCH 


Advantages: 

o Controller can tune itself to adjust to modeling errors 

o Controller self tuning is restricted to maintain stability 
of nominal plant as defined be the model 


Disadvantages: 

o Controller may not be able to tune itself sufficiently 

o Analytical analysis difficult because controller not fixed 
for a given flight condition 
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Discussion of Merits and Maturity 
of Top Level Approaches to Controller 
Design for Flexible Structures 
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AIRCRAFT MODAL SUPRESSION SYSTEM: EXISTING DESIGN APPROACH 

AND ITS SHORTCOMINGS 

By 

J. K. Ho, T. J. Goslin.and C. B. Tran 
The Boeing Company 
Seattle, Washington 


ABSTRACT 


The bending of flexible body aircraft may degrade the ride comfort of passengers. This 
Is especially noticeable towards the aft end of the aircraft (due to the relatively large tail 
surfaces) which may easily be excited when flying through turbulence. In addition, 
some aircraft may experience a front body bending mode which can be annoying to the 
cabin crew and first class passengers. Normally, this dominant body bending mode 
falls between 1-5 Hz. This range is easily perceived by the human body. Also, in some 
situations, the rigid body control law may be out of phase with the mode and aggravate 
the vibration. Hence, an active modal suppression system is desirable for improving 
the ride quality of the airplane. 

The size of the mathematical model, which has both the airplane rigid body and 
flexible characteristics, could easily exceed 100 states. This paper addresses the 
computational burden and fidelity of this large structural model. Later, the design 
methodology of the control law, which could be categorized into three steps:-(l) sensor 
selection. (2) modal phase determination and (3) modal suppression filter design-will 
be discussed. Each of these steps will be discussed in detail. Then we will present the 
theoretical results and compare them with flight test results. Here we will highlight the 
shortcomings of this design approach and briefly discuss what can be done in light of 
these deficiencies. Finally, we will include a brief description of the software tools. 
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AIRCRAFT NODAL SUPPRESSION SYSTEM 
EXISTING DESIGN APPROACH AND ITS SHORTCOMINGS 

BY 

J. HO, T. GOSLIN, AND C. TRAN 
BOEING COMMERCIAL AIRPLANES 
SEATTLE, WASHINGTON 
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OBJECTIVES OF NODAL SUPPRESSION VAN DANPER 


0 CONTROL DUTCH ROLL RESPONSE 
o PROVIDE GOOD TURN COORDINATION 

o IMPROVE LATERAL RIDE COMFORT BY SUPPRESSING FLEXIBLE BODY MODES 
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DESIGN PROCESS 


o MODEL GENERATION 
o BASIC YAW DAMPER DESIGN 
o MODAL SUPPRESSION SYSTEM DESIGN 

o SELECTION OF SENSOR LOCATION, CONTROL SURFACE 
o ESTABLISHMENT OF MODAL PHASE RELATIONSHIP 
o MODAL SUPPRESSION FILTER DESIGN 


o FLIGHT TEST 
0 ITERATE IF NEEDED 



soe 


STRUCTURAL MODEL 


O ORDER OF STRUCTURAL MODEL EXCEEDS 60 STATES 

o WITH ANTI-ALIASING FILTER, TRANSPORT DELAY, SAMPLE AND HOLD, CONTROL 
LAM, RUDDER PCU, ORDER OF MODEL EXCEEDS 100 STATES 

0 MAIN FRAME COMPUTER REQUIRED FOR REASONABLE TURNAROUND TIME 
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Design Requirements 


The damping ratio for the closed loop airplane should be 0.4 or greater for 
the Dutch roll mode. Phase and gain margin requirements are: 


Mode frequency (fM) 

Gain Margin 

Phase Margin 

Hz 

dB . 

degs 

fM <= 0.06 

> ±3 

> ±20 

.06 <= fM <= 1st 
aeroelastic mode 

> ±4.5 

> ±30 

fM > 1st 
aeroelastic mode 

> ±6 

> ±45 



807 


.00016 t 


PSD 


.00014 

.00012 

.00010 

.00008 

.00006 

.00004 

.00002 


0.00000 k 


0.0 


Peak at 0.007 


M 


V'* 

T 


y " 


2.5 


5!o 


Frequency In Hz 


7.5 


PSD Level of the Lateral Acceleration at the Aft Galley 



Q08 


PSD 



PSD Level of the Lateral Acceleration at the Pilot Station 
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PSD 
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.00032 

.00028 

.00024 

.00020 

.00016 

.00012 

.00008 

.00004 
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YAW DAMPER CLOSED LOOP AIRPLANE 
PSD OF THE LATERAL ACCELERATION AT THE PILOT STATION 


WITH BASIC YAW DAMPER LOOP CLOSED, USE MPAC TO COMPUTE OBSERVABILITY 
AND CONTROLLABILITY 


f LATERAL ACCEL AT PILOT STATION ^ 

PICK { ) AS SENSORS 

1 LATERAL ACCEL AT AFT GALLEY ' 

PICK RUDDER AS CONTROL SURFACE 


DESIGN AFT FILTER FIRST BY DETERMINING ITS MODAL PHASE RELATIONSHIP 
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Plant Model lor Aft Filter Design 
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CD 
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-60 
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Bode Plot of tho Transfer Function From Rudder Command 
to Aft Galley Sensor 



Phase In 
Degrees 


-200 



Bode Plot of the Aft Filter 





Plant 


Plant 


Yaw Damper 




Front Sensor 


el for Pilot Station Filter Design 
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Gain In 
dB 



Bode Plot of the Transfer Function From Rudder Command to 
Pilot Station Sensor 



Phase In 
Degrees -200 



Bode Plot of Pilot Station Filter 



I Modal Suppression I 
I System I 


Closed-Loop System 
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PSD of the Lateral Acceleration at the Aft Galley 
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FLIGHT TEST RESULTS 
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MODELLING DIFFICULTIES 


o LACK OF ACCURACY IN STRUCTURAL MODEL 

o RUDDER ACTUATION SYSTEM INCLUDES NON-LINEAR ELEMENTS: 

RATE SATURATION, POSITION LIMITER, HYSTERESIS, DEADZONE, VARIABLE 
GAIN 

o USE OF TRANSFER FUNCTIONS RESULT IN LOSS OF ACCURACY IN FREQUENCY 
DOMAIN 
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STRUCTURAL STABILITY AUGMENTATION 8Y8TEM DESIGN USING BODEDIRECT: 

A QUICK AND ACCURATE APPROACH 

By 

T. J. Goslln and J. K. Ho 
Boeing Commercial Airplane Company 
Seattle. Washington 


ABSTRACT 


Both all and forward sections of long slender aircraft suffer from lateral accelerations 
in turbulence. Some of these accelerations can be attributed to flexible body bending of 
the airplane in the 1-6 Hz region. Given an accurate flexible body model of the aircraft, 
a control law can be employed using lateral acceleration feedback loop augmented with 
the basic dutch roll yaw damper system to actively damp out these modes. 

In this paper a methodology will be presented for a modal suppression control law 
design using flight test data instead of mathematical models to obtain the required gain 
and phase Information about the flexible airplane. This approach will be referred to as 
BODEDIRECT. The purpose of the BODEDIRECT program is to provide a method of 
analyzing the modal phase relationships measured directly from the airplane. These 
measurements can be achieved with a frequency sweep at the control surface Input 
while measuring the outputs of interest. The measured "Bode-models" can be used 
directly for analysis in the frequency domain, and for control law design. Besides 
providing a more accurate representation for the system inputs and outputs of interest, 
(his method is quick and relatively inexpensive. 

To dale, the BODEDIRECT program has been tested and verlfled for computational 
Integrity. Its capabilities include calculation of series, parallel and loop closure 
connections between Bode-model representations. System PSD, together with gain and 
phase maiglns of stability may be calculated for successive loop closures of multi- 
inpul/multi-output systems. Current plans include extensive flight testing to obtain a 
Bode model representation of a commercial aircraft for design of a structural stability 
augmentation system. 

In addition to the BODEDIRECT approach, an indirect approach using flight test data to 
derive a mathematical mode for analysis using a Transfer Function Matching Routine 
will be presented along with the strengths and weaknesses of each approach. 
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STRUCTURAL STABILITY AUGMENTATION SYSTEM 
DESIGN using BODEDIRECT; A QUICK and ACCURATE APPROACH 

By 

T. J. GOSLIN and J. K. HO 


BOEING COMMERCIAL AIRPLANES 


LZQ 


THE RIDE QUALITY PROBLEM 


1 Rigid Body Airplane 

a) Roll 

b) Pitch 

c) * Yaw 

d) * Lateral Motion 

e) Vertical Motion 

f) Longitudinal Motion 


2 Structural Airplane (Long Slender Bodies) 

a) * Body Bending Modes 

b) Tortional Modes 

c) Wing and Empennage Modes 


THE MODES WHICH MUST BE ACTIVELY CONTROLLED INCLUDE 


BASIC YAW DAMPER : 

The Dutch Roll Mode 

Spiral Mode for turn coordination 
MODAL SUPPRESSION SYSTEM ; 

1st and 2nd Body Bending Modes 

All other modes must be passively controlled or have 
sufficient gain and phase margins. 


Note : This discussion is limited to the sysnthesis of a modal 
suppression system using the BODEDIRECT program. It is assumed 
that a good basic yaw damper already exists. 
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TRADITIONAL METHOD 


THEORETICAL 

ANALYSIS 


WIND TUNNEL 
AERODYNAMIC 
DATA 


FLIGHT TEST 
DATA 


LARGE STATE - SPACE 
STRUCTURAL MODELS 


NUMBER NJ BODE PLOTS 

CRUNCHING X| AND PSD 


CONTROL LAW 
DESIGN 


FIGURE 1 A 


BODEDIRECT 



FIGURE IB 
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Mathematical Models 

Ad vgnteqes -i 

1) No prototype is necessary 

2) Analysis In both time and frequency 

domain 

3) Observability/controliability 

directly available 

4) Eigenvalues & damping ratios 

directly available 


Disadvantages : 

1) Require large computing 

budgets (main frame computer) 

2) Lack required precision 


BODEDIRECT 

Disadvantages : 

1) A prototype must exist 

2) Analysis in frequency domain only 

3) Observability/controllability 
not directly available 

4) Eigenvalues/damping ratios 

not directly available 


Advantages ; 

1) Quick and inexpensive (work 

station environment) 

2) Accurate 

nasa workshop/v3*t jgosXin 





BODE -MODEL GENERATION 


FREQUENCY SWEEP 


<b 

U) 



Outputs 

(Used in BYD and 
MSS control laws) 


30-60 SECONDS 


Bode-models are measured for all inputs and 
outputs of interest for the free-airplane 





MODAL SUPPRESSION SYSTEM 
CONTROL LAW DESIGN 



For control law design, we determine the 
modal phase relationship in the feedback 
loop such that a compensator can be synthesized 
to bring the feedback signal into phase with 
the desired mode(s). 
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DATA PRE-PROCESSING 


Each bode-model requires magnitude and phase at each frequency point. 
With the current program version, 

Up to 500 frequency points may be stored for the frequency range. The 
frequency vector may be loaded from a file or generated. 


Example: 


‘LOAD FREQUENCY POINTS 
201 _ 


Number of frequency points 


(from file) 


‘GENERATE FREQUENCY POINTS 
LOG 


Starting 


0.1,10.0,201 



point Ending point 


Number of frequency points 


Frequency points may be loaded or generated in log or linear scale. 



Similarly bode-models may either be loaded from a file 
in terms of magnitude and phase or generated from a 
transfer function. 

Example: 

*LOAD 'modelname' 

‘GENERATE ’modelname’ 

transfer 
function \ 
coefficients! 

In addition, a bode-model may be extrapolated from a previous 
bode-model or a set of lab test data using the command: 


‘CREATE BODE-MODEL ’modelname’ 

data f 

points to < 
extrapolate 1 


This allows quick bode-model generation from a few test data points. 


£34 


B0DE-M0DEL1 


B0DE-M0DEL2 



‘CONNECT IN SERIES : 

MAG - GAIN * MAGI * MAG2 
PHASE - PHASE1 + PHASE2 
IF( GAIN .LT. 0 ) PHASE - PHASE + 180° 

*CONNECT IN PARALELL 

MAG-V<COM 2 + YCOM ^ 

PHASE - TAN" 1 (YCOM/XCOM) 

where 

XCOM - GAIN1 * MAGI * COS (PHASE1) 

+ GAIN2 # MAG2 * COS (PHASE2) 
YCOM - GAIN1 * MAGI * SIN (PHASE1) 

+ GAIN2 * MAG2 * SIN (PHASE2) 

S3S 









‘RESOLVE FEEDBACK JUNCTION 


MAG = l/V XCOM 2 + YCOM 2 
FBAC 

PHASE - -TAN ' lYCOM/XCOM) 
FBAC 


OUT 


IN 



where 

XCOM - 1.0 - MAG (G*H) * COS (PHAS (G*H)) * GAIN 
YCOM - -MAG (G*H) * SIN (PHAS (G*H))* GAIN 


‘RESOLVE CLOSED-LOOP SYSTEM 
CLOSED-LOOP - FBAC * G (SERIES CONNECTION) 


‘COMPUTE BODE-LOCUS 

(This routine computes a locus of bode-plots for a range 
of specified closed-loop gains. The user provides G,H, and 
the range of gain values.) 


‘COMPUTE GAIN AND PHASE MARGINS 

(The gain and phase margins are determined from the broken loop 
system bode plot. For good stability the gain should be -6 db or q 
less at the ±180° crossings. Additionally, the phase should be ±45 
away from ±180° at the 0 db crossing. For clarity, all phase 
values are scaled between 0 and -360 in all bode-model plots. 
Uncertainty values for gain(db) and phase(deg) may be declared by 
the user. The program will warn of near 0 db or 180 deg crossings 
based on these uncertainties.) 
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EXAMPLE USING SIMPLE (YAW RATE ONLY ) BASIC YAW DAMPER 



■ , foH ' 

NYAFTB 

where 

G = act model 

MPR - " 1 + GH* 

NYPSR 

and H 

= ROR * K 


12 



•RESOLVE OPEN-LOOP SYSTEM 


THIS ROUTINE WILL RESOLVE THE OPEN LOOP BODE-MODELS 
(G1 , G2, G3...) FOR A MULTILOOP SYSTEM IF THE CLOSED-LOOP AND 
FEEDBACK BODE-MODELS ARE KNOWN. A GAUSS-JORDON 
SOLUTION IS REPEATED FOR EACH FREQUENCY POINT TO SOLVE 
THE SET OF SIMULTANEOUS EQUATIONS. 


U 



THIS ROUTINE MAY BE USEFUL WHEN THE FLIGHT TEST DATA IS 
OBTAINED WITH THE BASIC YAW DAMPER ON AND THE EFFECTS OF 
THE YAW DAMPER ON STRUCTURAL MODES ARE THE DESIRED 
INFORMATION. 
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BODEDIRECT INCLUDES NONLINEAR CAPABILITIES 


which incorporates a describing function (DF) analysis 


for the following : 

a) hysteresis 

b) rate or position limits 

c) deadzones 

d) variable gain 



common nonlinearities 
in control actuation 
(rudder) systems 


Using this analysis, different bode-models may be 
generated as a function of input amplitude. Using 
polar plots for the describing function a prediction 
can be made if limit cycling will occur and at what 
frequency and amplitude. Notice that BODEDIRECT 
may use actual lab test data for the rudder system 
if desired over the theoretical model. 
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CRUISE CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT PILOT STATION 
BODEDIRECT PERFORMANCE PREDICTION 



CRUISE CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT AFT STATION 
BODEDIRECT PERFORMANCE PREDICTION 



CRUISE CONDITION 
SINE WAVE SWEEP DATA 
LATERAL ACCELERATION AT AFT STATION 
BODEDIRECT PERFORMANCE PREDICTION 




POWER 

SPECTRUM 


CRUISE CONDITION 
SINE WAVE SWEEP DATA 
LATERAL ACCELERATION AT PILOT STATION 
BODEDIRECT PERFORMANCE PREDICTION 



FLIGHT TEST RESULTS - CRUISE CONDITION 


POWER 

SPECTRUM 

(G 2 ) 


CD 


POWER 

SPECTRUM 

(G 2 ) 



0.18 

0.135 


0.09 


0.045 



0 










Mathematical Models 



1) No prototype is necessary 


2) Analysis in both time and frequency 
domain 


3) Observability/controllablllty 
directly available 

GO 

4^ 4) Eigenvalues & damping ratios 

directly available 



1) Require large computing 

budgets (main frame computer) 


2) Lack required precision 


BODEDIRECT 



1) A prototype must exist 

2) Analysis In frequency domain only 

3) Observability/controllablllty 
not directly available 


4) Eigenvalues/damping ratios 


not directly available 



1) Quick and Inexpensive (work 
station environment) 


2) Accurate 


nasa workshop/v3-t jgoslin 




PSD 


DESCENT CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT PILOT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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PSD 


DESCENT CONDITION 
SQUARE WAVE SWEEP DATA 
LATERAL ACCELERATION AT AFT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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PSD 


DESCENT CONDITION 
SINE WAVE SWEEP DATA 
LATERAL ACCELERATION AT PILOT STATION 
BODEDIRECT PERFORMANCE PREDICTION 
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LATERAL ACCELERATION AT AFT STATION 
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OPTIMAL Q-MARKOV COVER FOR FINITE PRESICION IMPLEMENTATION 

By 

Darrell Williamson 

Australian National University Canberra 
ACT 2601. Australia 

and 

Robert E. Skelton 
Purdue University 
West Lafayette. Indiana 


ABSTRACT 


The existing q-Markov COVER realization theory does not take into account the 
problems of arithmetic errors due to both the quantization of states and coefficients of 
the reduced order model. All q-Markov COVERS allow sonde freedom In the choice of 
parameters. In this talk, we exploit this freedom In the existing theory to optimize the 
models with respect to these finite wordlength effects. 
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Optimal q-Markov Cover for 
Finite Precision Implementation 

Darrell Williamson* and Robert E. Skelton** 


Abstract 

The existing q-Markov COVER realization theory does not take into account 
the problems of arithmetic errors due to both the quantization of states and 
coefficients of the reduced order model. All q-Markov COVERS allow some 
freedom in the choice of parameters. In this paper we exploit this freedom in the 
existing theory to optimize the models with respect to these finite wordlength 
effects. 


♦Dept. of Systems Engineering, Research School of Physical Sciences, Australian National 
University Canberra, ACT 2601, Australia 

♦♦School of Aeronautics and Astronautics, Purdue University, West Lafayette, IN 47907, U.S.A 
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Introduction 

An asymptotically stable system can be characterized in terms of its impulse 
response sequence (Markov parameters) and its output covariance sequence 
(covariance parameters) due to a zero mean white noise input process. A general 
approach has been developed [3] for realizing a system which matches q Markov 
parameters and q covariance parameters. Such a system is referred to as a q- 
Markov COVER, and q-Markov COVERS may be generated from output data 
[3,4] or from higher order models [5,6]. The Markov and covariance parameters 
are not independent and consequently the q-Markov COVER is not unique. In 
particular, all q-Markov COVERS are not related by state space similarity 
transformations [4]. In this paper we shall exploit the remaining degrees of free- 
dom to optimize the q-Markov COVER realization with respect to an aspect of its 
finite wordlength realization. 

Specifically, when digital controllers are to be implemented, both the con- 
troller coefficients and the controller states must be represented in finite 
wordlength precision. This finite wordlength (FWL) representation (or quantiza- 
tion) causes inaccuracies in the response when compared to the ideal (i.e. infinite 
precision) behaviour. Effects of quantization on the controller are increased noise 
at the output due to internal state quantization, and errors in time and frequency 
response characteristics due to coefficient errors. 

In digital filter design, the FWL effects are known to be most significant 
when the poles of the filter are very close to the unit circle [12]. In particular, 
narrow band filters have all these poles near z = l±jo. For digital control, the 
zero-order-hold equivalent of a continuous time model (or controller) with a pole 
at X will have a discrete pole at exp (XT). Hence for fast sampling and/or low 
damping of the continuous models, the discrete model will behave like a narrow 
band filter. The synthesis of optimal digital controllers with respect to 
arithemetic quantization noise is an important consideration in design especially 
for continuous time systems operating under a fast sampling rate [9,10]. The 
effects of quantization depend highly on the structure of the controller. This 
paper seeks to reduce these errors in the synthesis of q-Markov COVERs. 

1. Discrete q-Markov COVER 

Consider the asymptotically stable nominal discrete system 


856 



-3- 


x(k+l) = Ax(k) + Bu(k) ; x^R*. u(k)eR"- 

y(k) = Cx(k) ; y(k)eR ,v (L1) 

where {u(k)J is a zero mean process with unit intensity E |u(k)u*(j)} = 15^ and 
E (x(k)u*(j)} = 0 for k £ j. The Markov parameters Mj and covariance parameters 
Rj of (1.1) are defined by 

= CA'B ; Rj = CA j XC*, j £ 0, Rj = CXA* j C*. j £ 0 (1.2) 

where the state covariance matrix X satisfies the Lyapunov Equation 

X = AXA* + BB* . (1.3) 


These parameters Mj and Rj appear as coefficients in the expansion of the transfer 
function H(z) and power spectral density H(z)H*(z -1 ); that is 

H(z) = C(zI-A)- 1 B = £M i z-< i+1) ; H(z)H*(z) = £ Rjz" j 

i=0 

We suppose that as data we are given the first q-Markov and first q-covariance 
parameters {Mj, Rj; i = 0, 1, .... q-1) of an asymptotically stable system from 
which we construct the two data matrices 

£ R W 


D q ±R q -M q A/*eR n » qxn ’ q 

where R q , M q and M q are the Toeplitz matrices of the data as defined by 

Ro Ri* ... K-i 
Rj Ro ... R«j-2 


R n = 


R<t-2 

[Rq-1 Rq-2 ••• Ro J 


" q = A 


[ 0 0 ... 0 0 


Mo 0 ... o' 

O 

O 

• 

O 

£ 

,M q = 

M, Mj ... 0 

Mo ••• • 

• • • • 

• # • • 

• • * 

• ♦ • 

« « • ••• 


Mq_ 2 0 

Mq-j M<j_3 ... M,, 0 


Mq_i Mq_2 ••• Mo 


(1.4a) 


(1.4b) 


The first data matrix D„ in (1.4a) is Hermitian and it is shown in [3-4] to be 
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positive semidefinite. Hence we can obtain a (nonunique) full rank factorization 


D q = P q P *; P,eRW 

(1.5a) 

where 


r q = rank (D q ) = rank(P q ) £ n y q 

(1.5b) 

If we partition P q according to 


P q = [E q F q ]; E q eR n ^, F^R^ 1 ^”' 

(1.6) 

then it follows that the second data matrix D q can be factored as 


D q = P q P q ; P q eR n ^ 

0.7) 

where 


P q = [F* G*]; GqeR^ 

(1.8) 


for some G q (to be determined). The following result has been established. 
Theorem 1.1 [3] 

Given the q Markov parameters {M^ i = 0,1, ..., q-1 } and the q covariance 
parameters {Rj; i = 0,l, .... q-1} and a matrix G q in (1.8) such that (1.7) is 
satisfied, then the realization { A q , B q , C q } of order r q defined by 

A q = P q P q ; B q = P q [M q ••• M^]*; C q = E q (1.9) 

where P q denotes the Moore-Penrose inverse of P is a q-Markov COVER. The 
corresponding controllability grammian is given by 

X q = I (1.10) 

Furthermore 

P q = [C q A*C* • • • (A q -1 )*C q ]* (1.11) 


aa □ 

This theorem describes a large but not complete class C q of q-Markov COVERS 
parameterized by {G q } such that for some E q ,F q the data matrices D q , D q satisfy 
(1.5)-(1.8). Each matrix G q will (generally) result in a q-Markov COVER having 
a different transfer function. In order to compute the set of all such G q , observe 
in (1.5)-(1.8) that 
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D ,= 


[Eq Fq] . 


(1.12a) 


Then 



Dq-1 



' F ql 

°q = 


dqq 

" 



[F q * G q *J 


(1.12b) 


implies 


dqqe/?^ 

EqEq* = Rq , FqFq = 5^, F q G q * = d,, G q G q * = dqq (1.13) 


Now expand D_ in terms of its singular value decomposition 




Zi 0 


ur 

0 0 


u 2 * 

t - 




Then from (1.12a) 


(E,*Fq*) = Z 1 W U; 


(1.14) 


(U5) 


so that E q = C q is defined by the first n y rows and F q by the last (q-l)n y rows of 
Ui£i 1/2 . Define 

p q = rank (F q ) . (1.16a) 

Then from (1.15) 

p q £ min (r q , (q-l)n y ) . (1.16b) 

Next, expand F q in (1.13) in terms of its singular value decomposition. If strict 
inequality occurs in (1.16b) we have 

- r .*1 


Fq=[U 0 Upl 


0 0 


Vu 

V P 




pi* p« 


The Moore-Penrose inverse F q of F q is then given by 

F, + =v„2 q - i u; 


(1.17) 


(1.18) 


Corollary 1.1 
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\ 


Define 

(i) (1.19) 

(ii) G q2 e such that G q2 G q2 = d^-^D^dq 

where 

s q = rank [d,, - d^D^Lj d,] (1.20) 

and 

(iii) G q3 = Vp e . (1.21) 

Then if strict inequality occurs in (1.16b) the set of all G q which satisfy (1.13) are 
given by 

G q = G q | + G q2 U q Gq 3 (1 .22a) 

where 

U q e RV (r ^ ; S q <; r q - p q £ ny (1.22b) 

is an arbitrary row unitary matrix (i.e. U q U* = I). Furthermore, if the Moore- 
Penrose P q of 

Pq = [E q *F q T (1.23) 

is expressed as 

Pq = [L„ L 12 ] ; L n e R r ^- l)n r , L 12 e« w (1.24) 

then the corresponding state space representation { A q , B q , C q } of the q-Markov 
COVER is given by 

Aq = Ln + L 12 G q ; L|j = LuF q € R 

B q = P q + [MoX * • ‘ 5 C q = E q ‘ d-25) 

If r q = p q , then G q = G ql is unique. 

Proof: The expression forF q G q in (1.13) implies G q is of the form 

Gq* = F q + dq + G*,M*; 

for some M. Then expanding G q G q using (1.13) we have 
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d qq = dqVqVFjd, + d^(F^)*G q3 M* + MG^d, - MG q3 G q ‘ 3 M* 

Also from (1.13) and (1.21) 

(F+J’F^D*,, G, 3 g£ = I; (F£'g‘, = 0 (1.26) 

so that 

MM* = d qq -d q *(F q + )*F q + d q 

Since MM* has rank s q , 

s q = rank(G q2 G q2 ) £ r q - p q 


2. Optimal Finite Wordlength q-Markov COVER 

A fixed point finite wordlength realization of the idea) (i.e. infinite precision) 
q-Markov COVER (1.1) shall be referred to as a q-FWL Markov COVER and is 
described by 

x(k+l) = AQ[x(k)] + Bu(k) 

y(k) = CQ[x(k)] (2.1) 

Q[x(k)] = x(k)-e(k) 

A A 

where e(k) is the error in computing x(k). The components of the matrices A, B, 
£ are assumed to have a W 0 bit fractional representation obtained by quantization 
of the components of A, B, C in (1.1). The components of x(k) have a W+W 0 bit 
fractional part while components of Q[x(k)] and u(k) all have a W bit fractional 
part The components of the state residue vector e(k) has a W+W 0 bit fractional 
representation in which the most significant W bits are zero. The LHS and RHS 
of (2.1) are therefore consistent with respect to their fractional wordlength 

* A 

representation. The number of bits required to represent the integer parts of A, B 
and C depend on the dynamic range of the coefficients. State space structures in 
which all coefficients are less than unity are therefore advaintageous in this regard. 
The required integer representation of Q[x(k)] will depend on the dynamic range 
of the input signal u(k). Inadequate dynamic range will result in arithmetic 
overflow. The accuracy in the computation of x(k) is determined by its fractional 
wordlength W. 

Define the state error vector e,(k) and output error vector £y(k) by 
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e,(k)^*(k)-x(k); e,(k) £ ?(k) - y(k) (2.2) 

Then from (1.1), (2.1) and (2.2) 

e,(k+l) = Ae,(k) - Ae(k) + AAQ[x(k)] + ABu(k) + BAu(k) (2.3) 
£y(k) = Ce*( k) - Ce(k) + ACQ[x(k)] 

where 

AA = A - A; AB = B - B; AC = C - C 
Au(k) = u(k) -u(k) 

There are five terms which contribute to the output error (i) internal arithmetic 
errors e(k) due to state quantization (ii) coefficient errors due to errors AA in A 
(iii) AB in B (iv) AC in C, and (v) input quantization errors Au(k). Under weak 
’sufficiently exciting’ conditions on the input {u(k)J it can be shown [6] that if 
Q[»] in (2.1) denotes ’roundoff quantization, then (e(k)} is a zero mean uniform 
white process with covariance 

B(e(k)e‘(k))=Y 2 I; Y 2 =^- 2 - 2W . (2.4) 

Similarly { Au(k)} is assumed to be a zero mean white uniform process with 

E (Au(k)A*u(k)} =y*I (2.5) 

We assume that the quantized coefficients A, B, C are obtained by rounding A, B, 
C to W 0 bit fractions. Consequently, all components Ap of the error matrices AA, 
AB, AC satisfy 

IApl <Y„; Y„=y2~ W '. (2.6) 

For simplicity we normalize the error matrices and define 8A, 8B, 8C by 

8A= — AA; 8B = — AB; 8C= — AC (2.7) 

Y 0 Y 0 Yo 

so that all components 8p of SA, SB and SC satisfy 

I8pl<l. (2.8) 

The steady state output error covariance Y of (£y(k)) is then given by (we assume 

independence of e(k), e(k) and x(k)). 
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Y = CPC* + fee' + y 0 2 (5C)(XVi)(5C)* + y 0 f[C$Q* + (8C)C*] , (2.9) 

where 

P = £(e x (k)e x *(k)} 

= APA* + ^AA* + y 0 2 (8A)(XVl)(5A)* + y 2 <8B)(8B)* + ^BB* 
and 

X = E (x(k)x*(k)} = AX(A)* + t 2 A(A)* + (lV)BB* 

For the remainder of this section we assume no coefficient errors (i.e. Yo = 0 in 
(2.9)) and consider only the effects due to finite state wordlength (FSWL). The 
issue of coefficient error shall be resumed in Section 4. 


Theorem 2.1 

Define the output noise measure 

J- tr[Y] . 

Then for y Q = 0 

J = Y 2 {tr[K] + tr[B*KB]} 

where 

K = A*KA + C*C . 


Proof: From (2.9) 


where 


Now 


and 


so that 


Y = CPC*; P=APA* + y 2 Z = P + y 3 I 
Z = I + BB*; 


P = Y 2 £ A k Z(A k )* 

k=0 


K= £ (A k )*C*CA k 

k=0 


(2.10) 

( 2 . 11 ) 
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tr[CPC*]= 7 2 tr(ZK) . 


□□□ 

A fixed point q-FSWL Markov COVER corresponding to the (ideal) q- 
Markov COVER (1.1) is therefore described by 

x(k+l) = AQ[x(k)] + Bu(k) 


y(k) = CQ[x(k)] 


( 2 . 12 ) 


Q[x(k)] = x(k)-e(k) 

The output noise gain (q x ) due to state quantization and the output noise gain 
(rj u ) due to input quantization are defined by 

T\ x = tr[K]; TJ U - tr[B*KB] (2.13) 

The noise gain rj x generally varies with state space representation whereas r| u is 
independent of the coordinate basis. Specifically, consider the q-FSWL Markov 
COVER 


z(k+l) = A Q[z(k)] + B u(k) 


y(k) = C Q[z(k)] (2.14a) 

Q[z(k)] = z(k)-f(k) 

where 

A = T -1 AT, B = C = CT (2.14b) 

and Q[z(k)] has a W bit fractional representation. Assuming ’sufficient excita- 
tion’ by u(k), the state residue sequence (f(k)} in (2.14a) due to roundoff quanti- 
zation will again be a zero mean white uniform process with covariance ^1 as in 
(2.5). The corresponding output quantization noise gains T| z and rj u due respec- 
tively to state and input quantization are given by 

ri z = tr[KJ; (2.15) 

where B is given by (2.14b) and 

K Z = AK Z A* + C*C . 

But from (2. 11), = T*KT, so that 
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Ti z = tr[T*KT]; tj u = tr[B*KB] (2.17) 

Notice from (2.13) that the noise gain rj u due to input quantization errors is unaf- 
fected by a similarity transformation. Conversely the noise gain due to state 
quantization generally changes with co-ordinate bases. There is no change if T is 
unitary. The q-FSWL Markov COVER (2.14) is superior to the q-FSWL Markov 
COVER (2.12) if 

Tlz<Tl x . (2.18) 

However the comparison in (2.18) must be made under the assumption of identi- 
cal scaling of the states x(k) and z(k). Specifically, equal ^- scaling of gain a 
from a zero mean unit intensity white noise input u(k) to the state components 
Xj(k) of x(k) requires 

Xjj = aforallj (2.19) 

where Xjj denotes the jth diagonal component of the state covariance matrix X 
given by (1.3). Equal l 2 -scaling of gain a of components of z(k) in (2.14) 
requires 

Zj^oc; Z=AZA*+BB* (2.20) 

Equality in l 2 -scaling of representations (2.12) and (2.14) is equi valent to equality 
in the state dynamic range (i.e. number of bits in the integer representation of 
states) for a given probability of overflow. We now state a resu It which is impor- 
tant for establishing l 2 -scaling. 


Lemma 2.1 [8,9] Suppose M = M* > 0 is an nxn matrix. Then a necessary and 
sufficient condition for the existence of a unitary matrix V such that 

VMVjj = a for all j 

is 


tr[M] = not 


cm 

We have shown in Lemma 1.1 that different similarity transformations of an 
ideal q-Markov COVER corresponds to different factorization of the first data 
matrix D q in (1.5a). Our aim is to optimize this factorization. 

Definition 2.1 
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The Optimal q-FSWL Markov COVER minimizes the output quantization 
noise gain T) due to state quantization errors; that is 

TW = min tr[T*K q T]; T*T = A' 1 (2.21) 

^ T.G, M 

subject to the l 2 -scaling constraint: 

Ajy — a for all j (2.22) 

where the observability grammian Kq satisfies 

Kq = AqKqAq + C*C q (2.23) 

with ( Aq, B q , C q ) defined by (1.22)-(1.25). 

□□□ 

In corollary 1.1 we have shown that all the degrees of freedom available to 
select G q are confined to an arbitrary row unitary matrix U q . We now show how 
to optimize U q . 

Theorem 2.1 

a. The optimal q-FSWL Markov COVER (1.25) is defined by 

n opt = r q 1 min(tr[K^]) 2 (2.24) 

where U q e R s ‘ ,x ^ r,_p,) is an arbitrary row unitary matrix and Kq satisfies 
(2.23). 

b. The transfer function of the optimal q-FSWL Markov COVER has Hankel 
singular values given by the eigenvalues of Kq defined by the minimizing 
U q . 

c. Suppose Uq = U qo is the minimizing solution corresponding to the optimal 
G q = G qo in (1.22a). Let {Aq 0 , B qo , C q0 } be the corresponding state space 
realization in (1.24). Then the optimal q-FSWL Markov COVER has a 
(nonunique) state space representation {T^AqoT^ T^Bqo, C qo T} where 

T 0 = U 0 7t 0 V 0 * (2.25) 

such that 

(i) the unitary matrix U 0 is defined by 
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U 0 *K q0 U 0 = X 0 2 (2.26a) 

where 

Kqo = AqoK^A * 0 + C^Cqo ; 5^ = diag{a 1 2 0 ,o4...,o r ^ ) ) (2.26b) 

in which {o^} are the optimal Hankel singular values (eigenvalues of K^). 

(ii) 

"o = 4- <IOko)£,'‘ (2.27) 

k=l 

and (iii) V 0 is unitary such that 

r < 

(VoZoV^ = — for all j (2.28) 

r q 

Tlopt ~ Tl„ (optimal) = ^-(£(g 2 (2.29) 

a r q k=l 

Proof: By corollary 1.1 we have for G q defined by (1.22) for any row unitary 
matrix U q (of appropriately specified dimensions) that 0 q defines a q-Maikov 
COVER. The corresponding realization { A q ,B q ,C q ] for each such U q has identity 
controllability grammian and observability grammian defined by (2.23). Now 
given a particular U q , apply a similarity transformation 

T = U 0 Jt 0 V * 

to the given q-Markov COVER. Then 

tr(T*K q T) = tr(x 0 2 U 0 *KqU 0 ) 
and 

(T*T) -1 = V 0 Jto 2 Vq 

By lemma 2.1, the l 2 -scaling constant can be satisfied for some V 0 provided 
tr(7t" 2 ) = not. Following Williamson [1, Theorem 4.1] (with a minor modification 
of the 1 2 - scaling constraint), the optimal performance is given by 
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r ' -> 

„ _ cc=l 

nw< " a% 

where {a 2 } are the eigenvalues of Kq. That is, 

M 

The optimal q-FSWL Markov COVER therefore achieves the minimum in (2.24). 
The structure of U 0 , 7 t 0 , V 0 in (2.25)-(2.29) follow directly from Williamson [1] 
(see proof of Theorem 4.1 with U = I). 


3. Computation of the Optimal FSWL Markov COVER 

Necessary conditions for the optimal solution in Theorem 2.1 can be 
obtained using the method of Lagrange multipliers. Specifically, let 

J = (tr[K^]) 2 + tr[A(-K q +A*K q A q + C*C q )] + tr[fi(I-U q U*)] (3.1a) 

where 

K, = K* K?‘ A = A* e /l w ; Cl = Cl* 6 R ***' (3.1b) 

are symmetric Lagrange multipliers. After taking derivatives of J using (1.22) 
and (1.25) 

= -K q + + ^qCq 


= 21 - 2AK* + 2A q AA q *K* (3.2) 

OKq 

Jj- = 20iL,* 2 K^ 1 AGi - 2£1U, 

By setting these derivatives to zero we obtain the following result 


Lemma 3.1 Necessary conditions for the derivation of the optimal q-FSWL Mar- 
kov COVER are given by 
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K q = A*K q A q + C*C q 
A = A q AA q * + K q "*; A = A*e/?W 
U q U q * = I ;Uq6 ^V (r ,-P^ (3.3) 

QU q -P q U q Q q = R q ; Q = Q*eR s ^ 

where 

P q = P q * = GqaL^L^Gqa e R 

Qq = Qq* = G q3 AG q3 € R <VP,)*fr,-P,> (3.4) 

Rq = G q2 L 1 * 2 K q (L 11 +L 12 G ql )AG < ;3 € R ****** 
and A q , Gqj, Ly are defined by (1.20)-(1.24) 

□□□ 

These necessary conditions cannot be solved explicitly lor the optimal row 
unitary matrix U q and so an iterative procedure is required. One possible algo- 
rithm is now described. 

Recursive Algorithm for Optimal q-FSWL Markov COVE R: 

(0) Set j = 0 and choose any row unitary U q (0) in (1.21a) 

(1) Form Aq(j) from 

AqO) = (Li i + Li 2 G ql ) + L 12 G q2 U q (j)G q 3 (3.5a) 

(2) Compute K^): Kj(j) = + C q C q (3.5b) 

(3) Compute A(j):A(j)=A q (j)A(j)A q (j)+K^^(j); A(j) = A*(j) (3.5c) 

(4) Compute P q (j), QqCj) Rq(j): 

P q (j) = G q 2L* 2 Kq(j)L 12 G q2 ; 0,0) = G q3 A(j)G q3 ; 

RqO) = GiL^KqOCLn+LjjGqjJAO^ (3.5d) 


i 
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(5) Update U q (j) by solving the nonlinear algebra problem: 
fl(j)U q (j+D - P q (j)U q 0+l)Qq(j) = RqO); fl(j) = n*0) (3.5e) 

U q (j+l)U q *a+D = I 

The most difficult step at each stage of the algorithm is to solve (3.5e) for a row 
unitary U q (j+1) and symmetric Q(j). There is generally no explicit solution 
except for the following special cases. 

Lemma 3.2 Consider the equation 

QU q -P q U q Q q = R q ; CleR s ^ (3.6) 

where 

P q = P q *6/?^; Qq = Q q eR q— pq) q (3.7) 

are given. Then there exists an analytical solution (Q,U q ) with Q symmetric and 
U q row unitary when s q = 1 or Qq = PL (p scalar) 

a. When s q = 1, 12 and P q are scalars and Rq is a row vector. Then U q is arbi- 
trary for Rq = 0 while for Rq * 0 

U q = RqCQI-PqQq)- 1 ; IIU q ll = 1 (3.8) 

b. When Qq = PL let RqR^ have the singular value decomposition 

P m 

fol °1 v r 
R q < = (v,v 2 )[ 0 0 J v . 

to * 

where is invertible. Then 

u, = O = K * VilqW (3.9) 

In particular, when RqR^ has full rank, 

U q = (RqRqYX (3.10) 

Proof: For case (a) 

Uq(m-P q Qq) = Rq ; Q scalar 

so that (3.8) follows by the row unitary property U q U q = I. In case (b) 
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(G - pP q )U q = Rq 

and using the row unitary property of U q 

(Q-pP q ) 2 = R q Rq 

Hence using the SVD of RqR^ 

VlZq^VrUq-Rq 

But Vj Vi = I and V* has full row rank which gives (3.9 >. 

□CO 

Strictly speaking, (3.8) is not an analytical solution since the scalar Q must still 
be chosen so that IIU q ll = l. Note that by Corollary 1.1, G q3 G q 3 = I so that 
Qq(j) = I in (3.5b) if A(j) = I. The necessary condition (3.5e) is equivalent to 
assuming Kq(j), A(j), P q (j), Qq(j) and Rq(j) are known and optimizing over row 
unitary U q (j+1). That is, after dropping the index j and j+1 in (3.5e) we have the 
following result. 


Lemma 3.3 Suppose P q , Qq and Rq in (3.7) are known. Then a necessary condi- 
tion for a row unitary matrix U q to achieve optimally for the problem: 

min tr[QqUqP q U q + 2R q U q ]; U q e (3.11) 

is that there exists a symmetric matrix Q such that (3.6) is satisfied. 

Furthermore, the optimization in (3.1 1) is equivalent to 

min J(U); Ue R (r<rP< ‘ )x(r<rP < ) (3. 12a) 

where 


J(U) = tr[QU*PU + 2RU] 

over unitary matrices U* = [Uq V*] where Q = Qq and 


(3.12b) 


P = P = 


Pq 0 




0 OJ 

R = [R 0 0] e /? (r< i _p ‘» )x(r< t _ P« ) 


nm 

The advantage of the point of view (3.12) is that U ctin be treated as a square 
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matrix. The solution to (3.12) when U is a 2x2 unitary matrix is provided in the 
following lemma. The result can be derived by directly substituting into (3.12). 


Lemma 3.4 Suppose P = P* = [py], Q = Q* = [qy] and R = [ry] are 2x2 matrices. 
Then the minimum in (3.12) over 2x2 unitary matrices U is achieved by either 

(i) U = diag{u],u 2 ) where u 2 = 1, u 2 = 1 minimize 

J i = r ll u l + r 22 u 2 + 2qi2Pi2UiU 2 (3.13) 


or(ii) 



where I x I £ 1 minimizes 

J 2 (x) = ax 2 + 2bx + 2(cx+d)Vl— x 2 


(3.14) 


a ~ (Pll-P22)(qil-<l22). b — r ll+ r 22 


c - qi 2 (Pn-P 22 ) + Pl 2 (q 22 _< lll)» d - r 2 j - r 12 


□□□ 


Note that we must optimize over the disjoint sets of 2x2 unitary matrices consist- 
ing of signature matrices (as in (3.13)) and rotations (as in (3.14)). The optimal 
solution of (3.13) can be obtained by inspection of the magnitudes of the 
coefficients in Uj. For example, suppose 

lr u l > lq 12 Pi2> ^ lr 22' 


Then 


ui = -sgn(r u ); uju 2 = -sgn(q 12 p 12 ) 

However the optimization in (3.14) requires numerical solution. 

A general nxn unitary matrix U is either a signature matrix (i.e. a diagonal 
matrix £ such that £ 2 = I) or a product of 1/2 n(n-l) rotations Uy where the 
components of Uy(k,l) Uy are defined by 

U ij (i,i) = U ij (jj) = cos0 ij (3.15a) 
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Uy(i,j) = -Uy(j,i) =sin 0y 
Uy(k,k) = 1 for k * i, k # j 

Uy(k,l) = 0 otherwise (3.15b) 

A particular signature matrix is also defined by (3.15b) where 

Uy(k,k) = ±l fork = ij 

Ujj(k,l) = 0 for k * 1 (3.16) 

By letting 

u = nu y 

The optimization in (3.12) can be reduced to a sequence of one dimensional 
optimizations over the angles 0y. To be complete, J(U) should also be evaluated 
separately for all 2 n (n = r q -p q ) signature matrices. A compromise during the 
iterative procedure is to include the possibility of compo nents Uy being defined 
by (3.16) as well as (3.15a). Rather than present the general result we only illus- 
trate by means of an example. 

Specifically, suppose we express a 3x3 unitary matrix U as 

U = U 12 U 13 U 2 3 (3.17) 

Then by invoking the trace property, J in (3.12b) can equivalently be expressed as 

J(Uy) = trfQyUyPyUy + 2RyUy] (3.18a) 

where 

Q 12 = U12U23QU23U13! Pj 2 = P; R 12 = b 23U13R 

Q 13 = U 23 QUi; P 13 = U 12 PU 12 ; R 13 = U^RU^ (3.18b) 

Q23=Q; p 23 =ur 3 u; 2 pu 12 u 1 3; R23 = RU 12 U 13 

With i = i 0 , and j = j 0 fixed in (3.18a), J can be optimization over Uy o . The pro- 
cedure is recursive. That is, first assume i = 1, j = 2 with U 13 and both initial- 
ized to (say) the identity. After optimizing over U 12 , fix U 12 and U 13 and optim- 
ize over U^j, etc. Many cycles may be necessary for convergence. 

In order to explicitly demonstrate the formulation for each of the 2x2 optim- 
izations consider the case i = 1, j = 2, and express 
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Ql2“ 

Ql2 Ql2 

2 3 

P 12 = 

1 2 

P 12 P 12 
2 3 

Rl2- 

1 2 
Rl2 R 12 

4 3 


Q 12 Q 12 


P 12 P 22 


P 12 P 22 


iii 

where Qi2» P i2* p i2 e p • Then from (3.15), (3.16) the optimal 0 12 which 
minimizes Ji 2 (U 12 ) also minimizes 

Ji 2 (0i2) = tr[Q 12 UeP 12 U e + 2(R 11 -K3 12 P 1 * 2 )U e ] 

where components of the 2x2 unitary matrix U 0 is defined by (3.15a) or (3.16) for 
i, j, e { 1,2) The 2x2 optimization of J 12 (0i2) over 0 12 is partially solved in lemma 
3.4. 

Before concluding this section it is important to reiterate that the dimension 
of the problem for optimizing over the row unitary matrices U q is generally low. 
In particular from (1.21b) both the number of rows and columns of U q is not 
greater than the number of outputs. For a single output systems, U q is a scalar 
and so there are at most two possibilities, and no optimization is necessary. That 
is, for p q < r q we merely evaluate the cost in (2.24) for two values of G q in (1.21a) 
corresponding to U q = ±1, while if p q = r q , then G q = G ql is unique. 

4. Coefficient Errors 

Recall that Y in (2.9) is the error in the covariance of the output {y(k)J due 
to finite precision implementation of both states and coefficients of the q- Markov 
COVER. The optimal q-FSWL Markov COVER minimizes the trace of Y when 
there are no coefficient errors (corresponding to y 0 = 0). Furthermore, when there 
are no coefficient errors, there are no errors in either the Markov parameters Mj or 
covariance parameters Rj in (1.2). Once coefficient errors are introduced and all 
finite wordlength (FWL) errors are considered, there is no longer a clear interpre- 
tation of what should constitute the optimal q-FSWL Markov COVER. One pos- 
sibility is to again attempt to minimize the trace of Y. Alternative performance 
criteria could be based on the errors AM t and ARj in the Markov and covariance 
parameters as given by 

Mj + AMj = (C+AC)(A+AA)'(B+AB); 

R + ARj = (C+AC)(A+AAyX(C+AC)* (4.1) 

where X satisfies X = AXA* + BB*. For example, one could attempt to minimize 


874 



- 21 - 


Cm~ £ tr( AMjfAMj)*] or Cr - £ tr|ARJ (4.2) 

i=0 i=0 

However there are no results which directly connect or Cr with errors in time 
or frequency response of the q-Markov COVERS. Furhermore, the analytical 
and computational aspects involved in the resulting optimization would be very 
difficult if not practically impossible. 

A convenient approach to parameter optimization is to assume a statistical 
model for parameter errors. A statistical design can be justified along the follow- 
ing lines. Suppose (as is the case in practice) that both the Markov parameters M* 
and covariance parameters Rj are known only to be accurate up to a specified 
wordlength, and any higher precisional representation is regarded as uncorrelated 
random noise. Then the calculation of all q-Markov COVERS (for a particular 
row unitary matrix U q ) will also only be accurate to a finite precision beyond 
which the parameter representation contains uncorrelated random noise. 

Lemma 4.1 Suppose M = M* > 0 and K = K* > 0 are given nxn matrices. Let 
vj e R n be a zero mean random variable uniformly distributed between ±1 with 
uncorrelated components which are also uncorrelated with components of Vj. 
Then we have 

E (vj*Mvj) = -i- tr[M] . (4.3) 

Furthermore 

E {tr[V*MVK] } = j tr[MK] (4.4) 

where 

V = [V|V 2 • • • v n ] eR™ . 

□□□ 

Unfortunately these results cannot be applied directly to (2.9) since X itself is a 

A. 

random variable. However if we approximate X by X we can deduce the follow- 
ing result 

Theorem 4.1 

Suppose the components of 8A, 8B and 8C are zero mean uncorrelated ran- 
dom variables uniformly distributed between ±1. Then E{J} where J = tr[Y] is 
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approximated byE{J) where 

2 2 

E { J} = ^tr[B*KB] + (V>+ y )tr[K] + ^-(tr[XK]+tr[X]) (4.5) 

where K, X are defined by (2.11) and (1.3). 

□CD 


Proof: From (2.9) ignoring the linear term in 8C 

J- 7 2 {tr[K] + tr[B*KB]) + 


Yo 2 {tr[(SA)*X(8A)K]+tr[(8B)*K(6B)] + tr[(8C)*X8C]} 


The result then follows using Theorem 2.1. 


□□□ 


Under a similarity transformation T, the performance measure (4.5) becomes 

2 2 

E {J T } ^Y 2 tr[B*KB] + (Y 2 +y)tr[T*KTl + y(tr[XK]+trrr 1 X(T- 1 )*] (4.6) 

Note that both tr[B*KB] and tr[XK] are invariant In fact the invariant eigen- 
values {o^} of XK are the squares of the Hankel singular values of the system 
defined by { A,B,C). Consequently we need only consider the minimization of 

2 2 

(T 2 + y)tr[T*KT] + y tr['T 1 X(T- 1 )*] (4.7) 

over similarity transformations T. We make use of an earlier result [8] to provide 
the minimum in (4.7). 


Theorem 4.2 [8] 

Consider a minimal asymptotically stable order system {A,B»C} with con- 
trollability grammian X and observability grammian K. Let X and K be the 
transformed grammians as a result of applying a similarity transformation T; that 
is 

X = T^XCT 1 )*; K = T*KT (4.8) 

Then 
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tr[a 2 X+K] £ 2a £ a k (4.9) 

k=l 


where {a k } are the Hankel singular values. Moreover equ ality is achieved in 
(4.9) if and only if 


K = a 2 X 

(4.10) 

In particular, in (4.7) 


min E (1 T ) = y’trlB'KB] + ^-(£ o}+2a^\) 

T 3 k=l k=l 

(4.11a) 

where 


a=Vl + 3(7^ 

(4.11b) 


The minimum value is achieved in (4.11a) when K, X satisfy (4.10) with a given 
by (4.11b) 


□□□ 

One optimal realization (4.10) is a scaled internally balanced structure; that is 

Xj = a^diagfoj, cr 2 a Iq ); Kj = a diag{o lf a 2 , • • • , a r<i ) (4.12) 

From the point of view of l 2 -scaling, equal diagonal components of X guarantee 
equal dynamic range of the state components. It is evident from (4.10) that any 
unitary transformation U applied to the coordinate basis having X and K as the 
respective controllability and observability grammians will not alter the optimal 
performance. Consequently an optimal realization in which all diagonal com- 
ponents of the controllability grammian are equal exists with controllability 
grammian U XjU and observability grammian U K^U such that 

tfXjUjj = — £ <J k for all j (4.13) 

OTq fc=l 

where Xj, Kj are defined by (4.12) and U unitary. The existence of U is 
guaranteed by lemma 2.1 and an explicit algorithm for constructing a (nonunique) 
U is available in [9, Appendix A]. 

Corollary 4.1 

The optimal q-FSWL COVER which minimizes (2.21) subject to the 1 2 - 
scaling constraint 
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A ii = ^£ 0 k f<* aUi (4.14) 

OT q k=l 

A 

also minimizes E { J T ) in (4.6) 

□□□ 

This result provides a connection between the optimal q-FSWL COVER structure 
which minimizes only the effects due to state quantization noise, and the subop- 
timal q-FWL Markov COVER structure which minimizes E { J T ) subject to the 
assumed random parameter error model stated in Theorem 4.1. Once again we 
note that the result is suboptimal in the sense that X and X in (2.9) and (4.5) are 
only approximately equal. The result of Corollary 4.1 is also only of academic 
value since the l 2 *constraint (4.14) is not known until the design is complete since 
the Hankel singular values {oj} depend on the optimal row unitary matrix U p as 
provided in Theorem 2.1. However a more explicit result can be stated. 

Corollary 4.2 

The optimal q-FSWL cover subject to the ^-scaling constraint (2.22) also 
minimizes E (J x ) in (4.12) subject to (2.22). 

□□□ 


5. An Example 

Consider a 5 mode simply supported beam of length n having 2 inputs u lf u 2 
and 2 outputs y lt y 2 

Ml = F(0.27i,t) , u 2 = T(tt,t) 
yi = 0(0, t) , y 2 = |i(0.6x,t) 

where F(0.27t,t) denotes a force applied at .2 tc units from the left end of the beam, 
T( 7 t,t) denotes a torque at the right end of the beam, 0(O,t) denotes angular 
deflection at the left end, and |i(0.6jc,t) denotes rectilinear deflection at 0.6 jc from 
the left end of the beam. The equations of motion are assumed to be described by 

% + 2 ^ 1 -COiJk + CD^Tlk = [sin(0.27tk) kcos(ick)] 


ui 

u 2 
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yi 

5 

kcos(Ortk) 

. y 2 

= Z 

k=l 

sin(0.6ttk) 


(5.1) 


where eo^ = k 2 rads/sec. and £ = 0.005. A continuous time 10th order state space 
model is defined by 


x = Fx + Gu , y = Cx 


where 


x=(Tii'ni ••• iisns)* 


(5.2) 


A zero order hold equivalent 10th order discrete model (1.1) is defined by 

T 

A = e Fr ; B = Je Fo daG 


For the numerical work, a sampling period T = 0.025 sec. was selected which 
corresponded to approximately 10 samples in the shortest period. The eigen- 
values of A are at 


0.996±j0.0250, 0.9985±j0.0500, 0.9968±j0.0750, 0.9945±j0.0998, 0.99 16+jO.l 246 . 


Using the algorithm described in Corollary 1.1 the following results were 
obtained. 


q 

S q 

r q 

Pq 


2 

2 

4 

2 

) 

3 

2 

6 

4 

[ U q is 2x2 

4 

2 

8 

6 

) 

5 

2 

8 

8 

) 

6 

2 

9 

9 

>no freedom 

£7 

2 

10 

10 

) 


Hence for q = 2, 3, 4, U q in (1.22b) can be an arbitrary 2x2 unitary matrix, while 
for q > 5 there is no remaining freedom in the q-COVER. 

Optimal q-FSWL COVER designs: 



COS0 q 

-sin0 q 


sin0 q 

cos0 q 


0 2 ~ 40° 
0 3 - 0 ° 
0 4 - 65° 


(other cases 


1 0 
0 -1 


and 


-1 0 
0 1 


were also checked and neither was optimal). 
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The cost ranges from (2.29) for a = 1 were 

Ti^ = 0.3143X10 6 £ T|2 £ 0.8478X10 6 

Tlsopt = 0.2570x1 0 6 ^ tl 3 ^ 0.4764xl0 6 
Tl4opi = 0.0019x10 s £ n 4 £ 0.1308x10 s 
The actual FWL output roundoff noise is given by 

where W bits are assigned to the fractional wordlength of the state. Hence a fac- 
tor of 4 improvement in ri q corresponds to a wordlength saving of 1 bit. There is 
little savings in this example when q = 2,3. However for q = 4 we have a saving 
of 4 bits. In practice, for fast sampling and low structural damping, the savings 
would increase as the dimension of the model increases (e.g. a simply supported 
beam of 50 modes with q = 8). 
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INPUT-OUTPUT ORIENTED COMPUTATION ALGORITHMS 
FOR THE CONTROL OF LARGE FLEXIBLE STRUCTURES 

By 

K. D. MInto 

GE Corporate Research and Development Center 
Schenectady, New York 


ABSTRACT 


Tills presentation will overview work in progress aimed at developing computational 
algorithms addressing two important aspects In the control of large flexible space 
structures; namely, the selection and placement of sensors and actuators, and the 
resulting multivariable control law design problem. 

'Hie issue of sensor/actuator set selection is particularly crucial to obtaining a 
satisfactory control design, as clearly a poor choice will inherently limit the degree to 
which "good" control can be achieved. Moreover, it is becoming increasingly clear that 
systematic methods are required for determining prior to the control law design phase 
whether a particular candidate sensor/actuator set will yield acceptable closed-loop 
performance, irrespective of the particular control system design methodology used. 

With regard to control law design we are driven by concerns stemming from the 
practical issues associated with eventual implementation of multivariable control 
laws, such as reliability, limit protection, multimode operation, sampling rate 
selection, processor throughput, etc. Naturally, the burden imposed by dealing with 
these aspects of the problem can be reduced by ensuring that the complexity of the 
compensator is minimized. 

Our approach to these problems is based on extensions to input/outiput oriented 
techniques that have proven useful in the design of multivariable control systems for 
aircraft engines. In particular, we are exploring the use of relative gain analysis and 
the condition number as a means of quantifying the process of senslor/actuator 
selection and placement for shape control of a large space platform. Complementing 
this activity is the development of a new multivariable design approach that allows the 
designer to precisely control the complexity of the resulting compensator. The 
technique Incorporates input-output performance criteria such as the popular singular- 
value loop-shaping approach, yet without resorting to high-order compensators 
inherent to observer-based design approaches. 
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OUTLINE 


Motivation for research 

Approach 

design philosophy 

focus: shape control 
Results 

preliminary experiences with 
of design problem: S/A set 
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Motivation 


• control of flexible structures recognized as a 
key emerging technology for GE 


• leverage considerable design experience with 
MIMO design process for medium complexity 
problems (aircraft engines) to high complexity 
systems (LFSS, IFPSC) 


• a particularly important unresolved issue: 
decoupling the process of sensor/actuator 
selection from control law design phase 


GOALS: 


- Quantitative, systematic approach to problem 
of decoupling S/A selection from control 
law design process 

- Complexity reduction/management 

- Design process 

- Final product 
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Complexity clearly a major issue ... 


(1) countably infinite number of S/A sets 

number 
placement 
- types, etc. 

(2) large dynamic-order models (many flex, modes) 


order reduction a critical step, due 
primarily to limitations in traditional 
control law design approaches 
( observer-based compensators imply 
high-order) 


model accuracy/fidelity often sacrificied to 
accomodate these inherent computational 
limitations (spillover effects, etc.) 


conflict with S/A selection process, where 
numerical behaviour improves with model 
dynamic order 


(3) Shape Control - very large I/O dimensionality 
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Approach 

traditional F.E.M. /MIMO design approach, based 
on the following cycle 

Focus: 

(1) Control Structure Design 

- selection, pairing 

(2) Control Law Design 

- MIMO design w/ fixed 

order compensators 

(3) Uncertainty Modeling 

(4) Computations 

- S/A placement 

- Frequency-domain control-law design 

- Stable Factorization (balancing, order reduction) 

- demonstrate via shape control problem for LFSS 
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Control Structure Design: 


" That portion of the control system design 
process which deals with the selection and 
pairing of measurement and manipulation 
variables " 

- S.O.A.: typically Ad - hoc, often arbitrarily 
chosen 


much iteration, involving control law 
design phase 

- probably the most critical step in entire design 
process (certainly true for shape control via 
MIMO techniques ... ) 

(1) Sensor/Actuator Selection 

- how many? 

- locations? 

- types* (* ignor for the present ... ) 


colocated vs. non-colocated 



- systematic search for candidates that 
guarantee "good" closed-loop control 
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(2) Pairings and Decentralization 


given a S/A set, how do we interconnect 
for minimal closed-loop interaction? 


assume standard unity-feedback configuration 


*o 


e 


u 


Decentralization: Choose C as follows ... 



x 


vs. 



vs. 


X X X X X X 

X X X X X X 

X X X X X X 

X X X X X X 

X X X X X X 

X X X X X X 


I/O Pairing: 



Example: 

Colocated, fully 
decentralized: 

ui <--> ei, for all i. 
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(1) S/A Selection 


High level algorithm ... 


- determine a large number of candidate 

sensor/actuator sets 

- reduce to a manageable number of acceptable 
control structures (design - by - analysis) 


Specifics ... 

- develop a necessary condition for assessing 
a candidate control structure based on 
stability robustness 


uncertainty characterization: modified additive 

perturbations ... 


G = G + AG 


5(G - G) 

3(G) “ 5 


- connect control structure design process with 
stability robustness via following accepted fact: 
" plants with low condition number are 'easy' 
to control ... " 


base selection process on condition number 
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we have the following result 


• •• 


Theorem: If C stabilizes P, then necessarily 


K(P) < 2 r- 

o( PC(I + PC) ) 


5*\ for 


0) < o) c 


assumes "perfect" control at DC ... 


Selection Process: 

- compute condition number of candidate structures 
at DC 

• discard those with large condition number ... 


Computational Aspects: 


RGA (relative gain array) yields lower bounds 
on condition number, hence a necessary 
condition for viable control structures 


computational burden of RGA calculation 
small, but problem with exponential growth 
in complexity required to examine all possible 
combinations 


- Example: > 3 x 10D10 ways to choose a 12 x 12 
control structure from a set of 20 
possible I/O pairs 
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Heuristic Solution: 


direct selection of inputs/outputs to minimize 
condition number 

based on SVD of plant DC-gain. 

/n 3(G) _ max || G u || 

K 5 (G) " min || G u || 

reduce condition number <-> (i) reduce max. sv 

(ii) increase min. sv 


- DC gain: G(0) = D + C(-A) _1 B 

- SVD: G(0) = U S V' 

- introduce notion of input/output coupling 
operators ... 



m = s V' 


N' 1 = s 2 u* 
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Note: 


z = M u, 


= N -1 


express u, y as sums of standard basis 
vectors, i.e. 

u ■ X °i®i ’ y s ZPi®i 


then we have 


z = ^ a.He. = ^ a. M. 


z = xpi N '’ - 2>;' 


Design Heuristic: 


drop those inputs (outputs) corresponding to 
the maximum and minimum gains from 
u to z (y to z), i.e. 


Il*ll 2 « XKI HMjij 

i 


encouraging results for shape control 
application ... 
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Some Observations : 


systems characterized by large condition numbers 
(in general), and relatively large uncertainty 

- extremely difficult to control 
• severe performance limitations 

tendency towards partially colocated feedback 
structures (as a result of S/A selection process) 


lower condition numbers with non-colocated, 
and in particular, non-square control structures 
( non-square systems a challenge for MIMO 
control law design ) 


numerical conditioning improves (generically) 
with increasing information content 

RGA: G(0) = C(-A) -1 B 

■ D Da - □ 

Rank(G(0)) £ dim( [J) 

K(G(0)) i as dim(A) T 
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Tools: 

RGA, BRG (Bristol, Arkun, Maniousiouthakis) 


Usage: 

- Assess interactions between various feedback 
blocks 

- Account for cross-feed degradation due to 
use of decentralized control structure 

Preliminary result: 12 x 12 system of previous 

example ... 



□ actuators 
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• sensors 



Control Law Design 


- loop-shaping design philosophy 

• account for practical implementation constraints 

- fault-tolerance and reliability 

- limit-protection and multimode operation 

- digital control aspects 

All severely impacted by complexity of 

compensator! 

- 2 - stage design procedure allowing explicit 
constraints on compensator complexity 

1st Stage: "Ideal” Compensator computation 

model-matching performance specification 
(sv loop-shaping basis) 

controller parameterization (IMIC, SF) 

identify constraints on achievable performance 
( inner-outer factorization) 

- compute "ideal" compensator, i.e. K* such that 


2nd Stag e: "Low-order" Compensator Computation 


- analysis of "ideal" compensator (frequency 
response) to determine 

(i) approx, required complexity (dynamic order) 

(ii) values of compensator denominator terms, 

initial values of numerator terms 


parametric optimization to adjust compensator 
numerator terms (least-square approach) to 
minimize 

II H - H d || 2 = 0 


Features: 


- NOT an open-loop order-reduction procedure 

- Closed-loop low (fixed) order design procedure, 
with a flavour of order reduction (Ideal comp. 
=> lower-order parametric design) 

- Frequency-domain oriented design, hence 
complexity proportional to I/O dimension, NOT 
state dimension 

- Reduced emphasis/reliance on explicit order- 
reduction 
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Balancing & Order Reduction (SFPACIO 


balancing a popular method for model/controller 
order reduction 


- also useful for avoiding large coefficients in 
state-space manipulations 


Problems: Unstable or marginally stable systems 


- decomposition solution: G s G + + G. 

- balance, reduce components, recombine 

- Stable Factorization approach: 

G = N D* 1 , N, D stable, right-coprime 

Form a composite system F, with state-space 
realization ... 


r - 


r A 

B 

D 


c 

E 

N 

1 

1 



c. 

E„ 



2 

2. 


( as stable as we wish ... ) 


Order-reduction on F: 








6 

A 


6 = A 6 - a reduced-order model of G 
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THE ACTIVE FLEXIBLE WING AEROSEB VOELASTIC 
WIND-TUNNEL TEST PROGRJlM 

By 

Thomas Noll and Boyd Perry 
NASA Langley Research Center 
Hampton, Virginia 


ABSTRACT 


The evolution of advanced, high performance aircraft Is rediuiring that the engineering 
disciplines of aerodynamics, controls, and structures be Integrated Into a unified 
aeroservoelastic technology. To provide for technology maturation, sophisticated 
analysis and design methodologies must be developed and Verified through data 
correlation with experimental results. The most economical! means of obtaining test data 
that Includes the effects of these three disciplines without actually conducting full-scale 
night tests Is through the use of flexible wind-tunnel models scaled for aeroelastic 
phenomena. For a specific application of aeroservoelastic technology, Rockwell 
International Corporation developed a concept known as th e Active Flexible Wing (AFW). 
The concept incorporates multiple active leading- and traillmg-edge control surfaces with 
a very flexible wing such that wing shape Is varied in an optimum manner resulting in 
improved performance and reduced weight. As a result of a cooperative program between 
the AFWAL’s Flight Dynamics Laboratory, Rockwell, and NJKSA LaRC a scaled aeroelastic 
wind-tunnel model of an advanced fighter was designed, fabricated, and tested In the 
NASA LaRC Transonic Dynamics Tunnel (TDT) to validate the AFW concept. Besides 
conducting the wind-tunnel tests NASA provided a design of an Active Roll Control (ARC) 
System that was Implemented and evaluated during the tests. The ARC system used a 
concept referred to as Control Law Parameterization which Involves maintaining 
constant performance, robustness, and stability while using different combinations of 
multiple control surface displacements. Since the ARC system used measured control 
surface stability derivatives during the design, the predicted performance and stability 
results correlated very well with test measurements. 

The wind-tunnel model described above serves as the basis of a follow-on program to 
validate LaRC's and Rockwell's aeroservoelastic analysis methodology and multifunction 
digital control law design capability. This program provides an excellent opportunity for 
NASA and Rockwell to obtain an experimental database for the subsonic, transonic, and 
low supersonic speed regimes on an advanced aircraft configuration and to obtain 
experience with digital control systems and simulation methods. Significant activities to 
be conducted by NASA LaRC during the next 2 to 3 years to support the program include: (1) 
the design of multifunction digital control laws for flutter suppression and rolling 
maneuver load alleviation acting singularly and simultaneously: (2) the design and 
fabrication of a digital controller and the Implementation a|nd coding of advanced control 
laws: (3) a "hot bench" simulation of a flexible model with unsteady aerodynamics to 
verify the functionality of the digital controller; and (4) ground vibration, control system 
functional, and wind-tunnel tests on a model with violent flutter characteristics. Besides 
providing a multimilllon dollar aeroelastic model for the program, Rockwell will design 
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and fabricate a wing "tip missile" capable of either inducing flutter within the TDT 
performance envelope or preventing flutter through a decoupler mechanism, assist In the 
development of the advanced digital control laws, and participate during the testing and 
evaluation phases. 

Some results from the two previous wind-tunnel entries which describe the ARC system 
and the Control Law Parameterization concept will be presented during the workshop to 
establish the background for the more advanced studies now being pursued. In addition, a 
status report on the follow-on cooperative program will be discussed covering all facets of 
the effort. 
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THE ACTIVE FLEXIBLE WING 
AEROSERVOELASTIC WIND TUNNEL TEST PROGRAM 

The Active Flexible Wing (AFW) AeroServoElastic (ASE) Wind Tunnel Test Program is a 
recently-initiated cooperative effort between the NASA LaRC and the Rockwell International 
Corporation. The objective of this effort is to develop the analysis, design and test methodologies 
required to apply Active Controls Technology (ACT) for controlling and exploiting the aeroelastic 
characteristics of a flexible aircraft to improve performance. The approach selected to accomplish 
the program objectives includes the demonstration of various ACT concepts on a flexible full-span 
wind tunnel model, and the testing of the model to obtain an exj>erimental data base for validating 
the analysis and design methodologies associated with ACT. This effort is being directed by the 
Aeroservoelasticity Branch of the Structural Dynamics Division at LaRC. 
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OUTLINE 


I. Program Objectives 

II. Background 

III. Current Program/Description of Tasks 

Controller Development 
Design of Tip Missile Concept 
Development of EOM 
Synthesis of RMLA and FSS Control Laws 
"Hot Bench" Simulation 
Wind Tunnel Model Testing 
V. Concluding Remarks 
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OUTLINE 


This chart shows an outline of the topics to be discussed. Since the effort has only recently begun, 
today's presentation is a status report on where we are today. I will begin by reviewing the 
objectives of the program to scope the tasks involved. I will then describe some of the NASA 
analysis and test results obtained during a previous AFW wind tunnel test program between the Air 
Force Flight Dynamics Laboratory, the Rockwell International Corporation and the LaRC. This 
portion of the presentation will demonstrate the requirement to go beyond what had been 
accomplished, pushing the state-of-the-art into more challenging arid rewarding areas for ACT 
application. Next, a few charts describing each of the major tasks associated with this program 
will be discussed along with the progress and milestones recently accomplished. Finally, some 
concluding remarks and projections will end the presentation. 
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DEMONSTRATION AND VALIDATION 
OF ACTIVE AEROELASTIC CONTROL TECHNOLOGY 



ACTIVE FLEXIBLE WING MODEL MULTIFUNCTION DIGITAL CONTROL SYSTEM 


Objectives : 

• Demonstrate Multifunction Control Laws and Validate Analysis Methodologies 

• Obtain Experimental Data Base on Advanced Aircraft Configuration 

• Obtain Experience with Digital Control Systems and Simulation Methods 

NASA Role : 

• Design and Fabricate Digital Controller 

• Design Advanced Multifunctional Digital Laws and Code Controller 

• Perform Real Time Simulation for a Flexible Vehicle Model 

• Conduct Ground Vibration, System Functional, and Wind Tunnel Tests 

Rockwell Role : 

• Provide Flutter Critical AFW Model ($M's) for Technology Demonstration 

• Assist in Development of Advanced Control Laws 

• Provide Consultation and Participate During Test and Evaluation Phases 
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DEMONSTRATION AND VALIDATION 
OF ACTIVE AEROELASTIC CONTROL TECHNOLOGY 


The objective of this effort is to demonstrate the potential of using multifunction active control laws 
for controlling or exploiting aeroelastic response to improve aircraft performance. In addition, it 
gives NASA an opportunity to obtain an experimental data base on a flexible high performance 
advanced fighter configuration for validating analysis and design codes, to develop simulation 
techniques that include structural flexibility and unsteady aerodynamics, and to gain experience 
with digital control law implementation procedures. The NASA LaRC team consists of about 
twelve researchers from three different directorates (Structures, Flighi; Systems and Electronics). 
The team as a whole has the multidisciplinary experience required to perform the tasks identified 
for the AFW Aeroservoelasticity Program. The team will be required to design and fabricate the 
digital controller, design multifunctional control laws and code the controller, perform simulation 
studies to verify controller operation and conduct all model ground and wind tunnel tests. 

Rockwell through a Memorandum of Agreement and a separate support contract will provide the 
wind tunnel model for use during the program, assist in the development of the active control laws 
and participate during the wind tunnel tests. This cooperative effort provides an excellent 
opportunity to directly transfer technology to the aerospace industry. 
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ACTIVE FLEXIBLE WING WIND TUNNEL MODEL MOUNTED IN THE 
16-FOOT TRANSONIC DYNAMICS TUNNEL 


In mid 1985 the Air Force Flight Dynamics Laboratory in cooperation, with the NASA LaRC 
awarded a research contract to Rockwell International to test advanced control concepts on an 
aeroelastically scaled full-span wind tunnel model representative of an advanced fighter 
configuration. The model, shown in the photo mounted in the LaRC 16-Foot Transonic Dynamics 
Tunnel, was designed and fabricated by Rockwell using company funds. To give some 
perspective to the photo, the wing span from tip-to-tip is about 9 feet. The model consists of a 
rigid fuselage with scaled inertia characteristics and flexible wings. The wing box contains an 
aluminum honeycomb core and aeroelastically tailored plies of graphite epoxy. Each wing has two 
leading edge and two trailing edge control surfaces powered by rotary vane electrohydraulic 
actuators. The control surfaces have a chord and span of 25 percent of the local chord and 28 
percent of the wing semispan, respectively. They can receive constant signals remotely or time 
varying signals from a computer for active control investigations. Deflection limits are imposed on 
the various control surfaces to avoid exceeding hinge moment and wing load limitations. The 
model has the capability to roll about the sting axis or can be held fixed at any roll angle using a roll 
brake assembly mounted in the sting. In addition, the model can be tested at various pitch angles 
remotely controlled by an actuator located in the sting. All actuators <ire powered by an onboard 
hydraulic system. 
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INTERNAL DETAIL OF THE AFW MODEL 


This chart shows the wind tunnel model with some of the fuselage and wing panels removed to 
expose the complex internal detail required for ASE investigations. The outboard trailing edge 
control surface is driven by one actuator while the other three are driven by two each. Therefore, 
the control surfaces are powered by 14 actuators all supplied by onboard hydraulics. Eleven 
accelerometers ( five on each wing and one on the fuselage) can be used as sensors for active 
control systems or for monitoring model response during testing. In addition, the model has 
sixteen strain gages, nine rotary variable differential transformers (R\DT) to indicate control 
surface and pitch actuator position, a roll rate gyroscope and 141 static pressure taps on the upper 
and lower surfaces of the left hand wing along five spanwise stations. A six-degree of freedom 
force and moment balance is also present. 
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ACTIVE FLEXIBLE WING WIND TUNNEL TEST PROGRAM 
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ACTIVE FLEXIBLE WING WIND TUNNEL TEST PROGRAM 


The AFW concept integrate# active control# technology with a flexible #tructurc by using leading 
and trailing edge control surface# to deform the wing m an optimum manner to enhance 
aerodynamic performance and control . Two wind tunnel test# were previously conducted to 
validate the AFW concept The purpose of the First test (March and April, 1986) was to measure 
static aeroelastic and flexiblized stability derivative data as model angle-of-attack and control 
surface deflection# were varied. Some typical comparisons between an experimentally determined 
stability derivative and the predicted (unconected analysis) value as a function of dynamic pressure 
is shown in the upper right hand comer of the figure. The corrected analysis results were 
determined by using two separate "correction factors" in the analysis The first factor was used to 
match the expected rigid value of the stability derivative (extrapolation of the experimental data to 
the zero dynamic pressure value). The second factor, which varied with dynamic pressure 
(flexibility effect), was used to match the experimental values of the stability derivative with - 
dynamic pressure and to match the reversal conditions for each appropriate control surface. These 
factors were then employed during the design of the active roll control law which was evaluated 
during the 2nd test period. The predicted performance for the roll control law design is shown in 
the lower right portion of the figure. For both flight conditions analyzed, the predicted 
performance exceeded the goals established. In addition to these tests, the model was flutter tested 
for safety considerations across the planned flight envelop even though the model was designed to 
be flutter free. 
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AEROSERVOELASTIC ANALYSIS VALIDATED BY 
WIND TUNNEL TESTS 

Mach » 0.90 Dynamic pressure « 250 psf 


Active flexible wing model 


Roll rate 
(rad/sec) 


Roll rate performance 

I — 

Command 
/"P" — Response 


Block diagram of active roll 
control system 


R , ol l r ate 

WTM Re sponse 


0 0.3 12 

Time, sec 

Control law parameterization 

o lAnalyJ E) 

Or x 



Analy. 

Expt. 

6, 


o 

6, 

— 

• 


e> V 1 

(degrees) „ K _ = n 7fi 





AEROSERVOELASTIC ANALYSES VALIDATED BY WIND TUNNEL TESTS 


The 2nd wind tunnel test period was conducted during February and March, 1987 to evaluate 
Active Roll Control (ARC), Maneuver Load Control and Structural Mode Control Systems 
developed by Rockwell under Air Force sponsorship. In addition, an /VRC system designed by 
NASA using the aerodynamic corrections factors discussed previously and a "parameterization" 
procedure was evaluated in the tunnel. This "parameterization" procedure allowed the designer the 
flexibility of maintaining a fixed closed loop stability and a fixed closed loop roll performance 
while using different commanded control surface deflections for the active surfaces involved. The 
consequence of this concept is that the deflections of one pair of contro l surfaces can be traded-off 
against the deflections another pair with no loss or gain in aircraft stability or system performance. 
This idea becomes very important when control surfaces are required to undertake multiple active 
control functions simultaneously. The chart shows some of the NASA ARC system test results 
obtained. The lower right hand figure presents a sampling of data to demonstrate the principle of 
the concept and illustrates the excellent correlation obtained between the test and the calculated data. 
The figure shows that by changing one control law parameter, stability and performance are 
maintained while different amounts of leading edge and trailing edge control surface deflections are 
used. 
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AFW SCHEDULE FOR FIRST WIND-TUNNEL TEST 
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AFW SCHEDULE FOR FIRST WIND-TUNNEL TEST 


Two wind tunnel test entries are planned for the current program; the 1st entry is scheduled for 
April, 1989 and the 2nd test, one year later. During the 1st entry active flutter suppression and 
rolling maneuver load alleviation systems will be demonstrated separately. A schedule showing 
the major activities prior to the 1st wind tunnel entry is presented on the chart These 
multidiscipline activities include: 1) the model digital controller design, acquisition/fabrication, 
checkout and software coding; 2) the design and fabrication of a wing tp missile device to cause 
flutter within the flight envelope of the model and to act as a flutter-stopper for safety purposes; 3) 
the development of the aeroelastic equations of motion for six different model structural conditions; 
4) the synthesis of the RMLA and the FSS active control laws; 5) the "hot bench" simulation of the 
digital controller and associated software; and 6) the appropriate ground testing of the model to 
define its structural and dynamic zero-airspeed characteristics. Each of these activities will be 
discussed separately in the following figures. Although the details involving the 2nd wind tunnel 
entry will not be discussed here, the goal is to demonstrate active flutter suppression while the 
model is undergoing rolling maneuvers and alleviating wing loads. 
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AFW CONTROL LAW DEVELOPMENT SYSTEM 


The control law development system for the AFW wind tunnel model includes a Sun 3/160M and 
a Sun 3/50M workstation, a 141 megabyte hard disk with a 60 megabyte tape backup and a Apple 
LaserWriter for printer output. The workstations, driven by the Unix operating system, are 
connected through a Ethernet line. This network provides an excellent environment for several 
people to develop software and implement control laws independently. The Sun 3/50M is only 
used during the control law development, implementation and coding phases of the program. To 
execute the control laws the Sun 3/160M workstation requires a SKY Challenger processor board 
and a SKY Warrior processor board to be attached to the VME bus. The SKY Challenger is a 
VME digital signal processor board that is required to perform the scheduling and interfacing of the 
control law to the AFW model during simulation and testing phases. The SKY Warrior is a VME 
array processor board which can be used by either the Sun or the Challenger for performing high 
speed floating point arithmetic. In addition, two DT 1401 (VME Data Translation Cards) each of 
which has 32 analog-to-digital converters (ADC) and 2 digital-to-analog converters (DAC), and 
two DT 1406 each of which has 8 DAC are required to interface the incoming and outgoing model 
signals. These four boards provide 64 ADC and 20 DAC for use by the control program. A 282 
megabyte SCSI is also attached to the bus for the storage of data. 
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AFW MODEL DIGITAL CONTROLLER 






AFW MODEL DIGITAL CONTROLLER 


This chart contains a schematic drawing of the "interface box" and the APiV Control and 
Development System the comprise the model digital controller. The interf ace box processes 
the signals coming from or going to either the wind tunnel model or the "hot bench" simulation 
through low-pass filters, anti-aliasing filters and electrical isolation networks. The purpose of the 
low-pass filters is to reduce the high frequency noise and to limit voltage spikes that might appear 
on any of the 64 analog input lines. Currently, a 4th order Butterworth filter with a cutoff 
frequency of 100 Hz is planned to be used during the wind tunnel tests for the anti-aliasing filters. 
To be compatible with the "hot bench" simulation computers, the cutoff frequency of the anti- 
aliasing filters will require time scaling. The 16 analog signals returning to the model or to the 
simulation computer will also be filtered to prevent sharp edge transitions from being sent to the 
actuators. The Development System consists of several components linked to the Sun 
Workstation. The SKY Challenger is required to command the SKY Warrior, control the 
management of the data acquisition system (reads the ADC and writes to the DAC), monitor and 
update the User Control Panel, check limits and act as the system timekeeper. As described on the 
previous chart the SKY Warrior performs the required high speed floating point arithmetic. 
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AFW FLUTTER BOUNDARY MODIFIED 
BY THE ADDITION OF A WING TIP MISSILE 


Goal: 

Problem: 


Demonstrate Significant Increase in Flutter Dynamic Pressure Using Digital Active Flutter 
Suppression System 

Basic AFW Flutter Boundary is Beyond the TDT Tunnel Limits 


Solution: Add Wing Tip Missile to Lower Flutter Boundary 


Benefit: Tip Missile Designed to be a Flutter Stopper for Safety 
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AFW FLUTTER BOUNDARY MODIFIED 
BY THE ADDITION OF A WING TIP MISSILE 


Since active flutter suppression is one of the concepts being investigated during the present 
program, it is necessary to modify the model so that it will have a flutter instability within the 
operational capabilities of the TDT. In addition, this flutter instability must occur at sufficiently 
low dynamic pressures such that flutter suppression may be demonstrated experimentally. Several 
options were considered for lowering the flutter speed of the wind tunnel model; the option most 
attractive was to add a wing tip missile. The tip missile significantly increases the wing pitch 
inertia while only slightly changing the wing total mass. This in effect decreases the zero-airspeed 
1st wing torsion and 1st bending mode frequencies, and brings the two frequencies closer 
together. Because of the aerodynamic/structural/inertia interaction, the two modes will coalesce 
and cause flutter at a significantly lower dynamic pressure than without the tip missile present The 
lower left figure on the chart shows typical flutter boundary calculations for the model with and 
without the tip missile present Because of the close proximity of a symmetric and an 
antisymmetric flutter boundary, the active flutter suppression system must be capable of preventing 
both flutter modes simultaneously if the concept is to be effective. An added benefit of using a tip 
missile for causing flutter is its ability, with a little ingenuity, of returning the model to a flutter-free 
and, thus safe condition. This is accomplished by decoupling the missile dynamics from the wing 
by the use of a soft pitch spring at the wing/missile interface. The decoupling mechanism could be 
two pins, one stiff and one soft, as shown in the figure to the right Wkh the two pins installed the 
model would be flutter critical; with the stiff pin retracted, the wing is decoupled from the missile 
and becomes flutter free. 
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EQUATIONS OF MOTION 


W (t) 

[M] {q(t)} + [D] {q(t)} + [K] {q(t)} + \ pv 2 [CKt)] (q(t)} + ^ pv 2 {Q G (t)}-j- = {0} 

FREQUENCY DOMAIN: 

• Fourier Transform EOM 

• Q(ito) and Q G (ioo) are transcendental tabular functions 
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NO 

00 LAPLACE DOMAIN: 

• Laplace Transform EOM 

b /b \ 2 6 sA m 

• Q(s) and Q G (s) approximated by: A 0 + A lV s + A 2 (. v s) + X , v„ \ 

m - 3 V, b H m - 2/ 

STATE SPACE EOM: 

• Inverse Transform Laplace Domain EOM 

• Obtain the first order representation: x = Ax + B c u c + B Q u G 
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EQUATIONS OF MOTION 


One approach to formulating the equations of motion of an elastic aircraft is based on a chosen set 
of vehicle vibration modes and the Lagrange energy equation. Considering only small 
perturbations from a level equilibrium flight condition, the aircraft can be represented by a set of 
linear equations expressed in terms of the generalized coordinates, q(t). An example of such an 
equation is provided at the top of the slide. This equation represents a summation of forces and 
includes the inertial, the dissipation and the internal restoring forces, and the aerodynamic forces 
caused by the aircraft's rigid body, control surface and flexible motions and caused by gusts. To 
determine the aeroelastic characteristics of the vehicle, these equations are classically transformed 
into the frequency domain so that state-of-the-art unsteady aerodynamic theories based on simple 
harmonic motion can be used. These unsteady aerodynamic generalized force coefficients are 
transcendental tabular functions of several parameters, including freo[uency. Analyses in the 
frequency domain are straight forward using common methods. To ; perform aeroelastic analyses 
and design studies that include the effects of active feedback control systems, the equations of 
motion are transformed into the Laplace domain. The transformation of these equations into the 
Laplace domain is complicated by the transcendental functions of the generalized forces. The use 
of rational functions to approximate the generalized forces provides one solution to this problem. 
Several procedures are available for determining the rational function approximations. The 
equation shown on the slide is one of the more common forms. Here, a least square fit of the 
aerodynamic data is performed to determine the coefficients of the polynomial for each element of 
the frequency dependent generalized force matrices. Once the transformed generalized forces are 
obtained, the equations of motion are then placed into state-space foim for design investigations. 
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ANALYSIS FLOW 
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ANALYSIS FLOW 


This chart illustrates the major analyses performed and the flow of data and information between 
analyses. Circles represent both input to and output from the various analyses; rectikngular boxes 
represent analyses. The starting point, at the upper left, is a circle containing lumped-mass 
matrices and either stiffness or structural-influence-coefficient matrices. These matrices have come 
from a structural analysis code (not shown) and go into an eigenvalue/eigenvector analysis yielding 
in-vacuum frequencies, mode shapes and generalized masses. These quantities then go to three 
other boxes, the first of which is labelled Aeroservoelastic Analysis. Within this box the open- 
loop (and, when control laws are available, the closed-loop) aeroelastic equations of motion are 
generated, various analyses are performed, and intermediate results are passed "downstream" to a 
Control Law Synthesis box and two Simulation boxes. When generated, control laws are passed 
back up to the Aeroservoelastic Analysis box for computation of closed -loop frequency responses, 
closed-loop time responses, closed-loop flutter, etc. Control laws and other data are also passed to 
the two Simulation boxes which ultimately provide a functionality check of the digital control 
computers and a "best" pre-test estimate of the stability and performanc; of the closed-loop wind 
tunnel model. 
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SYMMETRIC FLUTTER ROOT LOCUS 
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SYMMETRIC FLUTTER ROOT LOCUS 


Various analysis procedures can be used to obtain the aeroelastic charac teristics of the model. This 
chart shows typical stability results using a root locus approach. To adequately define the flutter 
stability for the AFW wind tunnel model it was necessary to develop eq uations of motion for six 
different model representations. These included the model undergoing symmetric motion with the 
tip missile attached to the wing with either the stiff spring or the soft spring as discussed 
previously, and the model undergoing antisymmetric motion with the tip missile attached to the 
wing with either the stiff spring or the soft spring with the roll brake or and off. The data shown 
on the plot represents the AFW model undergoing symmetric motion with the tip missile attached 
to the wing with a stiff spring. Mach 0.9 doublet lattice unsteady aerodynamics were used for 
these calculations. Velocity was held constant and the air density was varied so that a matched 
point solution was obtained For this analysis, the first ten vehicle elastic modes were used to 
define the generalized coordinates. The predicted flutter mode involves the coalescence of the 2nd 
and 3rd elastic modes at a dynamic pressure of 213 psf at a frequency cf 1 1.1 Hz as can be seen 
when the 2nd elastic mode root moves into the right half plane. The objective of the FSS is to 
move the unstable root back into the left hand plane of the plot. 
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SYNTHESIS METHODOLOGIES 


1 . Classical Analog Design ■ Digital Realization ► Constrained Optimization 

1 

Final Digital Control Law 


2 Analog LQG ► Order Reduction ► Digital Realization 

l 

Final Digital Control Law ^ Constrained Optimization 


3. Direct Digital Design 


4. Eigensystem Design Techniques 
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SYNTHESIS METHODOLOGIES 


Several approaches that include classical analog and modem/optimal control techniques are being 
evaluated for use in designing the digital active flutter suppression system. The classical 
techniques being considered are based on root locus, Bode (transfer functions) or Nyquist plots 
and are useful for single-input, single-output systems. Once an analog control law which provides 
at least minimum stability for the aircraft at a certain design point is found, it is transformed into the 
z-domain and then optimized based on constraints such as design loads, actuator deflection and rate 
limits, and stability margins. The optimization task results in improved stability margins and 
robustness characteristics. The Linear-Quadratic-Gaussian (LQG) metiod is a systematic approach 
for designing multi-input, multi-output control laws. The LQG methoc. is based on minimizing a 
cost performance index consisting of quantities such as control deflection, design loads, 
accelerations, etc. The control law developed using this technique, however, is the same order as 
the aircraft being modeled. For flexible aircraft with unsteady aerodynamic forces, the number of 
states required to represent the vehicle is usually quite large. This order problem is solved through 
the truncation or residulation of the Kalman Filter. As described above, the LQG reduced-order 
control law is transformed into the digital domain and further optimized to improve performance 
and robustness. A third approach being considered involves the direct digital design of the control 
law. The methodology for the direct synthesis (determination of the coefficients for the z terms) of 
the digital FSS uses constrained optimization, and will meet multiple design requirements if 
necessary while maintaining reasonable stability requirements. The las : method being evaluated is 
an eigensystem design technique. The method involves the placement of the closed loop roots to 
obtain a control law with satisfactory stability, performance and robust less characteristics. 
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SYSTEM OVERVIEW OF HOT BENCH SIMULATION 


Wind tunnel model - B648 
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SYSTEM OVERVIEW OF HOT BENCH SIMULATION 


The purpose of the "hot bench" simulation is to provide a comprehen sive evaluation of the 
functionality of the Sun digital controller and the user control software, and to provide a low order, 
linear check of the flexible/dynamic system coupled with the active control laws. It is planned that 
this activity will be accomplished by attaching the Sun digital controller to a Cyber 175. The Cyber 
represents the AFW aeroelastic equations of motion modeled to include a sufficient number of 
elastic modes and the unsteady aerodynamic forces needed to accurately predict the static and 
dynamic characteristics of the test article across its expected test envelop. The Cyber will send 
sensor and other model or test condition information to the Sun for pocessing by the digital 
controller and will receive control actuator displacements from the Sun. Issues which can be 
investigated during the "hot bench" simulation besides the user control software, and control law 
stability and performance evaluations include: 

1) the operation of the flutter stopper, 

2) actuator transfer function differences between left and rig] it wings which could cause 
coupling between symmetric and antisymmetric model ch<iracteristics, 

3) failed actuators and sensors, 

4) control surface displacement and rate limits. 

A schematic that demonstrates the procedure to provide the interface between the Sun digital 
controller with the Cyber 175 during the "hot bench" simulation or to the AFW model during the 
wind tunnel tests is shown on this chart. On a previous chart, the NASA/Rockwell Interface box 
was discussed. Recall that this box takes discrete and analog signals from the wind tunnel model 
or from the Cyber during the "hot bench" simulation studies. Becaus e of the clock step and time 
step differences between the Sun and the Cyber, it will be necessary to conduct the simulation in 
synchronized slow time. 
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PRELIMINARY AFW MODEL TEST PLANS 


o Ground Vibration Tests 

o Control System Functional Tests 

Open Loop Transfer Function Measurements 

o Servoelastic Coupling tests 

Closed Loop limit Cycle and Ground Resonance 

o Wind Tunnel Tests 

Measure Control Surface Stability Derivatives 

Obtain Static Pressure Distributions at Selected Wing Locations 

Perform Roll Maneuvers with and without Active Load Alleviation 

Define Flutter Free Test Envelop for Decoupled Tip Missile 

Define Passive Flutter Boundary 

Demonstrate Active Flutter Suppression 
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PRELIMINARY AFW MODEL TEST FLANS 


The wind tunnel model will arrive at LaRC during the summer of 1988 for extensive ground 
testing prior to installing the model into the Transonic Dynamics Tunnel . Initially, the model will 
be ground vibration tested with the tip missile attached to the wing using both the stiff and soft 
spring representations separately. Symmetric and antisymmetric elastic mode frequencies, mode 
shapes and structural damping coefficients will be obtained with the model roll brake on and off. 

In addition, all sensor signals expected to be used by the active control laws will be measured. 
During these tests all actuators will be hydraulically powered. Open loop end-to-end tests will be 
accomplished to obtain transfer functions over a broad frequency range for all control 
surface/sensor combinations using several different amplitude signals to evaluate the nonlinear 
effects. The transfer function of selected components, such as the actuators and sensors will also 
be measured separately. Closed loop tests will also be accomplished for each active control law to 
be evaluated in the wind tunnel. These tests will include limit cycle tests to measure gain and phase 
margins at zero airspeed, model stability evaluations following an impulse excitation, and ground 
vibration tests. The intent of these tests is to obtain measured data for validating math models at 
zero airspeed (without aerodynamics). The various math models will be corrected and the control 
system designs updated as appropriate prior to the wind tunnel tests. Finally, the chart shows the 
expected wind tunnel tests and the order of conducting these tests. Routine force, moment and 
static pressure data will be measured first. Next the performance of t ie Rolling Maneuver Load 
Alleviation System will be evaluated. The higher risk tests which include the passive flutter and 
the active flutter suppression tests for preventing a high frequency violent flutter mode will be 
accomplished last 


VALIDATION OF TOOLS 
FOR MULTIDISCIPLINARY TECHNOLOGY 


MULTIDISCIPLINARY AREAS 



Elasticity 


1. Aerodynamics 

2. Controls 

3. Structures 


5. Closed-Loop Stability and Control 

6. Servoelasticity 

7. Aeroservoelasticity 


PLANNED TESTS/VALIDATION 

Static Pressure Distributions 
Analog and Digital Transfer Functions 
Frequencies and Mode Shapes 
Control Effectiveness and Flutter 
None Planned 
Closed-Loop Ground Tests 
"Hot Bench" Simulation/RMLA and FSS 


4. Aeroelasticity 


VALIDATION OF TOOLS FOR MULTIDISCIPLINARY TECHNOLOGY 


In summary, aeroservoelasticity is a multidisciplinary technology tl at involves unsteady 
aerodynamics, active control systems and flexible structures. This chart illustrates the potential 
interactions of these three technologies with the aid of three intersecting circles to represent 
individual technical disciplines. ASE represents that area common :o all three circles. To 
adequately develop analysis and design tools for application to ASE , it is important and critical to 
validate software within each area of interacting technologies. Rec;ill, that one of the objectives of 
this program is to obtain test data for evaluating the usefulness and accuracy of our codes involved 
in the design of flexible vehicles. The approach being followed during this program is to obtain 
experimental data to validate each of the primary technical disciplines prior to proceeding to levels 
involving two interacting technologies or three (ASE for this case). In conclusion, the 
NASA/Rockwell AFW program began in October, 1987 and will continue for about three years. 
This presentation has been a status report that addressed: 

1) why the program is being pursued, 

2) where we are today, and 

3) what to look forward to in the coming months. 

For those interested in pursuing the progress of the program, additional status reports will be 
presented at various conferences and workshops during the program. 
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MODELING AND STABILIZATION OF LARGE FLEXI3LE SPACE STATIONS 


By 

S. S. Um and N. U. Ahmed 
University of Ottawa 
Ottawa, CANADA 


ABSTRACT 


In this paper we present a preliminary formulation of a large space structure. The 
system consists of a (rigid) massive body, which may play t he role of experimental 
modules located at the center of the space station and a fledble configuration, 
consisting of several beams, which is rigidly attached to the main body. The equations 
that govern the motion of the complete system consist of several partial differential 
equations with boundary conditions describing the vibration of flexible components 
coupled with six ordinary differential equations that describe the rotational and 
translational motion of the central body. 

In our investigation we consider the problem of (feedback) stabilization of the system 
mentioned above. 

nils study Is expected to provide an insight Into the complexity of design and 
stabilization of actual space stations. Some numerical results will be presented. 
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ACTIVE VIBRATION MITIGATION OF DISTRIBUTED PARAMETER, SMART-TYPE 
STRUCTURES USING PSEUDO-FEEDBACK OPTIMAL CONTROL (PFOC) 

By 

W. N. Patten, H. H. Robertshaw, D. Plerpont, R H. Wynn 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 


ABSTRACT 


A new, near-optimal feedback control technique Is introduced that is shown to provide 
excellent vibration attenuation for those distributed parameter systems that are often 
encountered In the areas of aeroservoelasticlty and laige space systems. The technique 
relies on a novel solution methodology for the classical optimal control problem. 
Specifically, the quadratic regulator control problem for a flexible vibrating structure 
Is first cast In a weak functional form that admits an approxima :e solution. The 
necessary conditions (first-order) are then solved via a time finite-element method. The 
procedure produces a low dimensional, algebraic parameterizatlan of the optimal 
conlrol problem that provides a rigorous basis for a discrete controller with a first- 
order "like" hold output. 

Simulation has shown that the algorithm can successfully control a wide variety of 
plant forms Including multi-input/multi-output systems and systems exhibiting 
significant nonllnearlties. In order to firmly establish the efficacy of the algorithm, a 
laboratory control experiment was Implemented to provide planar (bending) vibration 
attenuation of a highly flexible beam (with a first clamped-free mode of approximately 
0.5 Hz). Base actuation for the cantilever was accomplished using a three degree-of- 
freedom active bay (variable geometry truss) actuator. On-line processing was 
accomplished with a 14 mhz ’’AT’ type microcomputer with data acquisition capability. 
The results of the tests corroborate the utility of the method. 
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OVERVIEW OF ANALYTICAL DEVELOPMENT 
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CONCLUSIONS 


• ALGORITHIM VALIDATED 

• REAL TIME IMPLEMENTATION POSSIBLE 
ON LOW-LEVEL PC ARCHITECTURE 

(80286 WITH 27.5 KHZ D/A) 
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SHAPE CONTROL OF HIGH DEGREE-OF-FREEDOM VARIABIE GEOMETRY TRUSSES 


By 

R. J. Salerno, C. F. Relnholtz, and H. H. Robsrtshaw 
Virginia Polytechnic Institute and State University 
Blacksburg, Virginia 


ABSTRACT 


Common static trusses are constrained to permit no relative motion between truss 
elements. A Variable Geometry Truss (VGT1, however, is a tru:is which contains some 
number of variable length links. The extensible links allow the truss to change shape 
in a precise, controllable manner. These changes can also be used to control the 
vibrational response of a truss structure or to perform robotic tasks. 

Many geometric configurations, both planar and spatial, are possible candidates for 
VGT manipulators. For this presentation only two geometries will be discussed; the 
(hrce degree-of-freedom (DOF) spatial octahedral/octahedral tiuss and the three DOF 
planar tetrahedral truss. These truss geometries are used as the fundamental element 
in a repeating chain of trusses. This results in a highly dexterous manipulator with 
perhaps 30 to 60 degrees of freedom that retains the favorable stiffness properties of a 
conventional truss. From a fixed base, this type of manipulator could perform shape or 
vibration control while extending and "snaking" through com plex passageways oi 
moving around obstacles to perform robotic tasks. 

In order for this new technology to be useful in terms of robotJc applications the 
forward and Inverse kinematic solutions must be efficiently solved. The approach 
taken here is to first concentrate on fully understanding the fi irward and inverse 
kinematics of the fundamental elements and then utilizing the insight thus gained to 
solve the more complex problem of the kinematic chains. The Inverse solution of a 30 
DOF planar manipulator will be discussed. The discussion will focus on how to specify 
parameters for an underspecified system by using shape control algorithms. 
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The Position Control Problem 
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Spatial Truss Forward Kinematics 


Emp I oy 'the concept. of* 
\< i nemct. i ca I I y equ i va I ent 
dev i cee 




OOOf 


Lower Half of Spatial VGT 




( 00 / 


Constraint 


Li = 1-1,1-12 = 

L_g = I - ig n 1 3 = 

1 — 3 = n 3 n, = 



Equations 


f C S, , ©2 D 
fC0 2 , ©3 } 

f C©3 ,0! D 


i CO302 -Bs I n0j s i n0 2 +C 



ZOO! 


Continuing Research 



N90- 10123 


OPTIMAL INTEGRAL CONTROLLER WITH SENSOR FAILURE ACCOMMODATION 

B y 

Dr. T. Alberts 
Old Dominion University 
Norfolk, Virginia 

and 

Dr. T. Houlihan 
The Jonathan Corporation 
Norfolk, Virginia 


ABSTRACT 


An Optimal Integral Controller that readily accommodates Ser sor Failure - without 
resorting to (Kalman) Alter or observer generation - has been designed. The system is 
based on Navy-sponsored research for the control of high performance aircraft. 

In conjunction with a NASA developed Numerical Optimization Code, the Integral 
Feedback Controller will provide optimal system response even In the case of 
Incomplete state feedback. Hence, the need for costly replication of plant sensors Is 
avoided since failure accommodation Is effected by system sof.ware reconfiguration. 

The control design has been applied to a particularly lll-behavcd, third-order system. 
Dominant-root design In the classical sense produced an almost 100 percent overshoot 
for the third-order system response. An application of the newly-developed Optimal 
Integral Controller— assuming all state Information available- produces a response 
with NO overshoot. A further application of the controller design— assuming a one- 
thtrd sensor failure scenario— produced a slight overshoot response that still preserved 
the steady slate time-point of the full-state feedback response. 

The control design should have wide application In space systems. The design can be 
expanded to Include gain scheduling that enhances system response to large-scale 
transients. For this latter Instance, using the NASA optimization scheme, the 
guesswork normally required to determine feedback gains for large transients Is 
eliminated. 


1003 



Optimal Integral Contrc I 
With 

Sensor Failure Accommoda tions 


Dr. Thomas Alberts 
Old Dominion University 
Norfolk, Virginia 


Dr. Thomas Houlihan 
The Jonathan Corporation 
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Optimal Integral 
Control Design 

Introduction 

Optimal Regulator 

Augmented System 
(Rates of Change of Input Signals) 

Optimal Tracker 
Optimal Integral 
Control Design 

Sensor Failure 
Accommodation 

Preliminary Results: 

Third Order System 
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NASA WORKSHOP 


JULY 1988 


Introduction 


Optimal Control Designs Compromised By: 


Inaccessible States (Sensors) 
Noisy Feed back Signals 


OC Designs Resort To Use Of 
Filter / Estimating Techniques 
To Overcome These Obstacles; 


NAVY Research in 1970s 
Leads to Alternative Approach 
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Optimal Regulator - Classic Design 
Tradeoffs Between Accuracy of Control 
And Energy Expenditure Reflected 
In Weighting Matrices (Q and R) 

Of Performance Index (J) 

J = /( XTQX + UTRU) dt 

Performance Index Formulation 
Assumes Unconstrained Inputs. 

In Reality, Inputs are Limited. 
Futhermore, Rates of Change of 
Input Signals Are Limited. 
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Typical System Response. 
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FIGURE 3. OPTIMAL 
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Augmented System: 

Rates of Change of Input Signals 
Can be Considered 

New State Vector = Old State Vector 

+ 

Input Signals 

Optimal Regulator Solution of 
Agumented System: 

U* = -G*X* G* = IGi G2I 

Gain Matrix (G*) of Augmented 
System Carries Information on System 
States (X) and Inputs (U)! 
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Optimal Tracker: 

Add Gain Matrix (M) to 
Select Command Inputs 

NOTE: Tracker is NOT Integral Controller 
Since Control Commands are NOT 
Generated by Integral of Error 
Between Desired Signals (r) 

And Output Signals (z). 

NOTE: Solution to Tracker Control 

Configuration is KNOWN. It is 
Solution of Augmented System. 
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FIGURE 5. OPE I MAI. TRACKING SYSTEM 
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Optimal Integral Control Design 

Equality of Optimal Integral Control 
Design and Optimal Tracker Design 
Effected by Block Diagram Reduction 
Techniques (Laplace Domain). 


Results: 

I L H I = IGi 62 1 A B -1 

EC ED 

Knowns: 

A, B, C, D, E - Configuration Matrices 
Gi, G 2 - Augmented System Solution 
Thus: 

L and H Matrices are Determinable 
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Sensor Failure Accommodation 
Matrices: 

H = Error Gain Matrix 
L s State Gain Matrix 

If State Information Unavailable, 

Corresponding Column Elements of L 

Matrix Are Zeroed - Suboptimal Control! 
From Before 

IG 1 G 2 I = IL S HI A B 

EC ED 

Hence, New Gain Matrix IG S I = IGi s G2sl 
Can be Determined to Effect Control 

Preliminary Results are Encouraging 
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Fig. 9. Third Order 


(a) Conventional 
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(b) Optimal 
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Summary 

Optimal Integral Control Design 
Effected by a Combination of 
Multivariable Control Analyses 

Sensor Failure Accommodation 
Accomplished Without Resort to 
Supplemental Filter / Estimator Designs 

Suboptimal Control Response 
Effective for Ill-Behaved, 

Third-Order Test Case 
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Postscript to Computational Aspects... 


Lawrence W. Taylor, Jr. 

NASA Langley Research Center 

What started as an effort to transcend various project and 
reasearch activities has become an official program.. Computational 
Controls. The following charts describes that program at this early stage 
in its development. The next meeting on the subjects of the 
Computational Aspects Workshop will be the 3rd Annual Conference on 
.Aerospace Computational Control. The conference will be held August 
128-30, 1989 at Oxnard. 
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Computational Controls 


^OBJECTIVlTj 

To Develop the NEW GENERATION 

NIGH PEEEOEMANCE Aerospace 
Modeling, Control, and Simulation Tools 
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Contacts: 

Lee Holcomb - NASA HQ Code RC 

John Dibattista - NASA HQ Code RC 

Guy Man - JPL 

Larry Taylor - LaRC 

Harry Frisch - GSC 

Henry Waites - MSFC 

Ken Cox - JSC 
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justificationJ 


Current Practices in Formulating , 
Modeling and Simulating do not 

Meet today s needs. 

• Hypersonic Cruise Vehicles 

• Multi-Component Launch Vehicles 

• Aeroassisted Orbital Transfer 


Vehicles 
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^ORGANIZATION! 

Code RC 

Goddard J P L Marshall 

Johnson Langley 
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APPROACH 


C Reference Problems ]> 



CUser Expedient Software^ 


Advanced 
Modeling, 
Analysis & 
Control 
Capability 
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Reference Problems 

• Shuttle RMS 

• Earth Orbiting Satellite 

• Mini-MAST 

• Pinhole Occulter 

• Mariner Mark II 

• Optical Interferometer 

• Advanced Launch System 

• F - 18 Fighter 

• Trans-Atmospheric Vehicle 
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Advanced Formulations 


• Order(n) Algorithms 

• Distributed Parameter Modeling 

• Mass Referenced Modeling 

• Composite Modeling 


(LaRC] 

\ / 

(LaRC] 

v I / 


[LaRC] 

\ . J 


(LaRC) 

\ V 
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Parallel Proc 


• Multiple Processors 

• Array Processors 

• Benchmarking 




1037 


^OFTWAREJ 

• Macintosh-Like User Environment 

• Simultaneous Tasking 

• Real-Time and Off-Line 

• Modular (Particular Methods) (LaRC 

• Data Base Management 


interactive Graphics 



1038 


Related 



R & T Base 

Control Techniques 
Control Applications 
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PROPOSED LARC AERO TASKS 

1.1 

DYNAMICS INTEGRATION AND ADVANCED CONTROL THEORY AND 
MODELING 

• F-18 THRUST VECTORED HI- <X VERSION 

• TRANS-ATMOSPHERIC 


HIGH-ORDER, HIGH FIDELITY, NONLINEAR 
MATHEMATICAL MODELS OF HIGH PERFORMANCE AIRCRAFT 

ADVANCED MODEL ORDER REDUCTION METHODS 

ROBUST INTEGRATED CONTROL DESIGN METHODOLOGIES 



RAPID CONTROLLER IMPLEMENTATION METHODOLOGIES 
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MASS/INERTIA 


AERODYNAMIC 




REACTION CONTROL 
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Plant Order Reduction for Controller Synthesis 


Physical 

System 


Frequency Response 


Region of 
Interest 


Evaluation 

Model 


® High-Order 
• Uncertain 



Eval 

ROM 


Error 


Frequency 

Matches Evaluation Model 
in critical range 

Quantifiable errors 
outside region 


Benefits 


Enhanced physical 
insight 

Error structured for 
robust control 

Reduced-order 

controllers 
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ACTIVITIES 


• Advisory Committee /Quarterly 

• Workshops /Annually 

• Programmatic Status Repts /Quarterly 

• Technical Reports /As Available 



ANNOUNCEMENT & CALL FOR PAPERS 

3rd Annual Conference on Aerospace 
Computational Control 

Radlsson Suite Hotel August 28-30, 1989 Oxnard, California 


ANNOUNCEMENT of a CLASS on 

msamr" 






S3 A Control System Simulation for 
Flexible and Articulating Structures 


WIEIIliMs August 31,1 9 8 9( After Conference) 

Wimiz 3rd Annual Conference on 

Aerospace Computational Control 

Radisson Suite Hotel, Oxnard, CA 

(SMTIMz • Overview 

• Example Problems 

• Hands-On Experience 

• User s Manual 

No Charge for Class or Materials 
for Registered Conferees 

Larry Taylor 
NASA Langley 
804-864-4040 


/044 


NASA 

Na»cnj Ae<anajtcs and 
Scwce Aomnsjcaion 


Report Documentation Page 


1. Report No. 

NASA TM-101578, Part 2 


2. Government Accession No. 


4. Title and Subtitle 


Proceedings of the Workshop on Computational Aspects 
in the Control of Flexible Systems 


3. Recipient's Catalog No. 


5. Report Date 

May 1989 


6. Performing Organization Code 


7. Author(s) 

Lawrence W. Taylor, Jr. (Compiler) 


9. Performing Organization Name and Address 

NASA Langley Research Center 
Hampton, VA 23665-5225 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


8. Performing Organization Report No. 


10. Work Unit No. 


506-46-11-01 


1 1 . Contract or Grant No. 


13. Type of Report and Period Covered 
Technical Memorandum 


14. Sponsoring Agency Code 


15. Supplementary Notes 


This report is part 2 of two parts published as NASA TM-101578. 


16. Abstract 

This publication is a collection of papers presented at the Workshop on 
Computational Aspects in the Control of Flexible Systems held at the Royce 
Hotel, Williamsburg, Virginia, July 12-14, 1988* The papers address the 
formulation, modeling, computation, software and control for flexible spacecraft, 
aircraft and robotic manipulators. 


18. Distribution Statr oant 

Unclassified-Unlimited 
Subject Category - 18 


19. Security Classif. (of this report) 

1 

20. Security Classif. (of this page) 

21. No. of pages 

22. Price 

Unclassifed 

Unclassified 

505 

A22 


☆U.S. GOVERNMENT PRINTING OFFICE: 1 9 8 9X 29 .29 5/ 0 0 0 28 


17. Key Words (Suggested by Author(s)) 

Large Flexible Spacecraft 
Control, Strucutal Dynamics 


NASA FORM 1626 OCT 86 














